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Rate of Temperature Change in 
Short-Length Round Timbers 


By J. D. MacLEAN,! MADISON, WIS. 


This paper discusses factors affecting the rate of tem- 
perature change in wood and includes time-temperature 
curves for use in determining the temperature obtained at 
different points in short-length log sections when heated 
in steam or hot water. One group of curves shows the 
temperature distribution from the ends to the mid-length 
and at different points between the circumference and 
axis, of logs 12, 24, 36, 48, and 60 in. in diameter. Other 
charts show the time-temperature distribution at various 
points between the circumference and axis at the mid- 
length of log sections of any diameter between 12 and 60 
in. Examples are given illustrating the procedure for 
using the charts in finding the temperature at a par- 
cular point, under different heating conditions. 


INTRODUCTION 


REEN logs of various species are often heated in steam or 
(& in hot water to condition them for veneer-cutting. These 
logs are relatively short, and the ratio of diameter to 
length is large in comparison with round timbers in the form of 
posts, poles, and piling. In timbers in which the ratio of di- 
ameter to length is small, the effect of longitudinal heating can 
extend only a limited distance, while the portion between the end 
regions is raised to any desired temperature. Because of the 
comparatively short length of veneer logs, however, it is necessary 
to consider the effect of heat transfer in the longitudinal as well 
as in the transverse direction. The relative amount of heating 
that takes place in these two directions depends upon a number of 
variables, such as cross-sectional dimensions, length, heating 
period, and temperature of the heating medium. 

The purpose of this paper is to discuss methods of finding the 
rate of temperature change in timbers of this kind. No attempt 
is made to discuss the temperatures needed for satisfactory re- 
sults in veneer-cutting, for some woods apparently require less 
heating than others. When the optimum temperature is known, 
however, the data presented provide a guide for selecting the 
required heating period for logs of different species and of different 
dimensions. The data also furnish a useful aid in studying the 
effect of different heating conditions on veneer-cutting opera- 
tions. 

A large amount of experimental work has been done at the 
Forest Products Laboratory to study the rate of temperature 
change in wood when heated in different mediums and to in- 
vestigate the effect of variables, such as species, moisture content, 
density, or specific gravity, and variation in the rate of heating 
in the radial, tangential, and longitudinal directions. A number 
of papers have been published which discuss the results of these 
experiments (1 through 10).? 
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EFrect oF VARIABLES ON RATE OF TEMPERATURE CHANGE 


Heating Medium. Experiments with different heating medi- 
ums show that for the same heating period and the same tempera- 
ture conditions, steam heats wood faster than water, and water 
heats faster than such liquids as preservative oils (6). In 
general, water heated wood about 5 to 10 per cent more slowly 
than steam. All of the liquids, however, heated the wood some- 
what faster than hot plates (6, 8). The slowest rate of heating 
was obtained in air at a low humidity, but the rate was increased 
as the humidity of the air was increased (10). 

Variation in the rate at which wood heats in different mediums 
is apparently caused by such variables as type of surface contact, 
rate of circulation when liquids and gases are used, specific heat of 
the heating substance, and heat of vaporization as in the case of 
steam. 

Differences in Moisture Content. Wood seasoned well below 30 
per cent moisture content, which is about the fiber saturation 
point, was found to heat more slowly than green wood. Results 
obtained in heating green wood nevertheless indicate that 
differences in moisture content above the fiber saturation point 
have no significant effect on the rate of heating in the range of 
temperatures commonly employed (4). For practical purposes, 
therefore, it may be assumed that green timbers of any given 
species will heat at about the same rate at any moisture content 
above the fiber saturation point. 

Direction of Grain. Studies made of the relative rate of heating 
in the radial and tangential directions (at right angles and parallel 
to the rings, respectively) showed there was no important 
difference in the rate of temperature change for these two direc- 
tions, but the rate of heating in the longitudinal direction or 
along the fibers was found to be from 2!/, to 23/, times faster than 
in the transverse direction (2, 6, 7). For practical purposes it 
should be sufficiently accurate to assume that on the average 
the rate of temperature change longitudinally is about 2/2 times 
faster than it is across the grain or in a plane at right angles to 
the axis of the tree. This ratio has been used in the calculation 
of temperature data that will be discussed later. 

Density or Specific Gravity and Species. Experiments on speci- 
mens of the same species and on those of different species having 
different specific gravities showed that with any given heating 
medium, the rate of heating varied’ inversely with the specific 
gravity (6). The rate of temperature change is apparently in- 
fluenced by the density of the wood rather than by the species, 
which affects the results merely to the extent that different species 
may vary in density. 

The specific gravity of wood is usually based upon the volume 
at current moisture content and weight when oven-dry. 

Table 1 shows the average specific-gravity values for various 
species of hardwoods and softwoods, These values are based 
upon volume when green and weight when oven-dry and are 
convenient to use when employing the time-temperature curves 


2 Numbers in parentheses refer to the Bibliography at the end of 
the paper. 


g959l 
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AVERAGE SPECIFIC GRAVITY OF VARIOUS SPECIES 
OF SOFTWOODS AND HARDWOODS? 


Average specific gravity 


TABLE 1 


Species (based on volume when 
green, and weight when 
oven-dry) 
Hardwoods: ; 
Ash, white (Fraxinus americana)...........+.+5 0.55 
Aspen, bigtooth (Populus grandidentata)........ 0.35 


Basswood, American (Tilia americana).......... 
Beech, American (Fagus grandifolia)............ 
Birch, yellow (Betula lutea)...........04+.ee8- 
Chestnut, American (Castanea dentata)......... 
Elm, American (Ulmus americana)............- 
Elm, rock (Ulmus thomasii).................-. 
Hackberry (Celtis occidehtalis)..............005 
Hickory, mockernut (Hicoria alba)............- 
Maple, silver (Acer saccharinum)...........++4+ 
Maple, sugar (Acer saccharum)............++ 
Oak, red (commercial) (Quercus sp.)....... 0.08 
Oak, white (commercial) (Quercus sp.)...... 5.5 
Pecan (Hicoria pecan).....-....-+sssee esses es 
Sweetgum (Liquidambar styraciflua)............ 
Sycamore, American (Platanus occidentalis)..... 
Tupelo, black (Nyssa sylvatica) .........00..005 
Tupelo, water (Nyssa aquatica)..........0.000 
Walnut, black (Juglans nigra).........-.-s.s00- 
Yellow poplar (Liriodendron tulipifera)......... 
Softwoods: ; 

Baldcypress (Taxodium distichum)............. 
Douglas fir, coast (Pseudotsuga taxifolia)....... * 
Douelas fir, Rocky Mountain (Pseudotsuga taxi- 

HOLIA) wire aie picterninterlare re ereseaere rahe emir ad ene aratebeeaa de 
Bir, white (CAbies‘sp:) 2. hremiscse stat ate 
Hemlock, Eastern (Tsuga canadensis). . 
Hemlock, Western (Tsuga heterophylla) 
Larch, Western (Larix occidentalis) . 
Pine, Eastern white (Pinus strobus) 
Pine, jack (Pinus banksiana).. 
Pine, loblolly (Pinus taeda)... 
Pine, longleaf (Pinus palustris).... 
Pine, ponderosa (Pinus ponderosa 
Pine, red (Pinus resinosa)........ 
Pine, shore (Pinus contorta). 
Pine, shortleaf (Pinus echina 
Pine, slash (Pinus caribaea)... 
Pine, sugar (Pinus lambertiana). 
Pine, Western white (Pinus monticola)... oe 
Red cedar, Eastern (Juniperus virginiana)....... 
Red cedar, Western (Thuja plicata)............ 
Redwood (Sequoia sempervirens)............+4+ 
Spruce, Engelmann (Picea engelmanni)......... 
Spruce, Sitka (Picea sitchensis)................ 
Spruce, white (Picea glauca)................+.. 
Tamarack (Larix laricina)............... 
White-cedar, Northern (Thuja occidentalis) 


a 


* Data for other species can be found in “Strength and Related Properties 
of W code aor in the United States,’’ U. S. Dept. of Agriculture, technical 
etin 479. 


shown in the figures. The method of using the data on density, 
or specific gravity, is discussed in the following section. 


CoMPUTATION OF TEMPERATURE 


The rate of temperature change in any solid depends upon the 
diffusivity which, like conductivity, may be considered a constant 
for normal ranges of temperature. This factor is a measure of the 
rate of temperature change and may be defined as the change in 
temperature produced in a unit volume of the substance by the 
quantity of heat that passes in unit time through unit area of a 
layer of unit thickness and having unit difference of temperature 
between the faces. If the diffusivity is represented by the symbol 
h?, the specific heat by c, the density by p, and the conductivity 
by K, the relation of these variables is shown by the algebraic 

: K 

expression, h? = —. 

cp 

product of specific heat and density, and directly as the conduc- 
tivity. 

Fig. 1 shows the relation of diffusivity and the specific gravity 
of wood determined from tests made on green material heated in 
steam. As previously mentioned, water was found to heat a 
little more slowly than steam, and for the purpose of calculation 
it should be satisfactory to assume that the rate of heating in 
water is about 90 per cent as fast asin steam. On this basis the 
dotted line in Fig. 1 represents the relation of diffusivity and 
specific gravity for green wood heated in water. If, for example, 
green wood having a specific gravity of about 0.5 were heated in 
steam, the diffusivity would be determined from Fig. 1 as about 


Thus the diffusivity varies inversely as the 


JANUARY, 1946 
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Fig. 1 Revation or Dirrusiviry AND Spreciric GRAVITY 
0.000302, while the corresponding diffusivity for heating in water 
would be about 0.000272. These diffusivity factors are in inch- 
second units. 

Since the rate of temperature change depends to a certain ex- 
tent upon the method of heating, the diffusivity values deter- 
mined from experimental tests may be considered as values - 
which apply for the heating medium used. 

Wood is a substance that is naturally more or less variable in 
structure; hence temperature calculations cannot be made with 
the same precision that would be possible if it were a homogeneous 
material or a substance having the same physical properties in all 
directions. Nevertheless, the calculated temperature values, 
which are based upon average results from experimental tests, 
are sufficiently accurate for practical purposes, and whatever 
differences there may be between the computed temperatures 
and those actually obtained in a given timber will be small pro- 
vided care is taken to maintain the heating conditions assumed. 

In the early part of the heating period a variable amount of 
time is usually needed to bring the heating medium to the desired 
maximum temperature. A gradual change in the surface tem- 
perature is therefore obtained during this period, while the 
temperature calculations are necessarily based upon the assump- 
tion that the surfaces of the timbers are maintained at a constant 
temperature from the time the heating period is started. 

If the temperature rise is fairly uniform, a convenient way to 
make an adjustment for the lower average temperature to which 
the wood is exposed during the early period of heating is to 
assume that the full temperature is applied for one half of this 
period. For example, if 2 hr were taken to bring the heating: 
medium to the maximum temperature, the 2-hr period would be 
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counted as 1 hr at the maximum temperature. While this is an 
approximate method of adjusting for the period of variable 
temperature, it is generally satisfactory for moderate periods of 
changing surface temperature when the rate of temperature rise 
of the heating medium is reasonably uniform. A consideration 
of this lower average temperature at the beginning of the heating 
period naturally becomes of less importance as the total time of 
heating is increased. When the heating conditions are not con- 
trolled so that they are approximately uniform over most of the 


heating period, individual judgment must be exercised in deter- * 


mining the total time the timbers should be heated to compensate 
for the variable conditions. 


Factors INVOLVED IN TEMPERATURE CALCULATIONS 


Equation [4] in the Appendix was derived for computing the 
temperature data shown in Figs. 4 to 16, inclusive. Definitions 
of the symbols used and a brief outline of the development of this 
equation are also given in the Appendix. 

An examination shows that the following factors are included 
in the fundamental equation: 


Initial wood temperature = Uo 
Temperature of heating medium = U; 
Diffusivity factor for radial heating = h? 
Diffusivity factor for longitudinal heating = q? 
Heating period = T (sec) 
Radius of timber = a (in.) 
Length of timber = L (in.) 
Distance from nearest end surface, of point at which 
temperature is to be computed = Z (in.) 
9 Distance from axis, of point at which temperature is to be 
computed = r (in.) 
10 Temperature at the point (7, Z) under consideration = U 


ONOohwWNH 


Fig. 2 shows the distance from the end surface and the distance 
from the circumference at which temperatures were computed 
for plotting the time-temperature curves. Fig. 3 is a chart for 
finding the corresponding wood temperatures when the heating 
conditions are different from those assumed in computing the 
temperature data plotted. ‘ 

In the preparation of the time-temperature curves shown in 
Figs. 4 to 16, inclusive, the initial wood temperature U) was 
taken as 70 F, the heating-medium temperature U; as 212 F, the 
diffusivity factor for radial heating h? as 0.000270 (in-sec), and 
the diffusivity factor q? for longitudinal heating as about 2.5 
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(Dotted lines illustrate examples.) 


times h?, in round figures 0.000680, making the ratio slightly over 
2.5, or about 2.52. It may be noted in Fig. 1 that the diffusivity 
factor of 0.000270 corresponds to about that for green wood hay- 
ing a specific gravity of 0.5, when heated in water. 

Since the diffusivity for the longitudinal di- 
rection q? is taken as 2.52h?, the exponential 


TEMPERATURE CURVES IN GROUP | ARE TEMPERATURES FROM END SURFACE 
TO MID-SECTION A-B AT DIFFERENT DISTANCES FROM THE CIRCUMFERENCE, 


NUMBERED I, 2, 3,4, ANDS 


TEMPERATURE CURVES /N GROUP 2 SHOW TEMPERATURE DISTRIBUTION FROM 
CIRCUMFERENCE TO CENTER OF MIO-SECTION A-8, AND TEMPERATURE CURVES /N 
GROUP 3 SHOW TEMPERATURES IN MID-SECTION AT DISTANCES OF 0.258y°05048, 
0.75Q,AND & INCHES FROM THE CIRCUMFERENCE WHERE Q@:RADIVS /N INCHES 


Fie. 2 SxercH SHowine Distancn From Enp Surracn anp Distancr From 
CIRCUMFERENCE AT WHICH TEMPERATURES ARE COMPUTED 


function ‘ 

rd Gas + em 
in the general Equation [3] can be written in 
the form 

—heT (ae + ee) 


é 


and a similar substitution can be made for gq? 
in the terms of Equation [4] (Appendix). A 
change in the diffusivity factor for radial heat- 
ing will therefore represent a proportional 
change in the longitudinal diffusivity factor g*. 

Although the temperature curves given in 
the various charts were computed by assuming 
the particular diffusivities and temperature con- 
ditions cited, these curves may be used for any 
predetermined conditions of temperature and 
for any species having diffusivities different 
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from those assumed. The method of doing this is explained in 
the Appendix under the heading, “Method of Using the Fig- 
ures When the Diffusivity, Heating Conditions, or Both Are 
Different From Those Assumed in Computing the Time-Tem- 
perature Curves.’ Examples illustrating the use of the curves 
under different conditions will be given later. 


Data PLOTTED IN TEMPERATURE CURVES 


The temperature curves comprise three groups in which the 
initial wood temperature was taken as 70 F and the heating- 
medium temperature as 212 F. 

Group 1. Figs. 4 to 11, inclusive, show the computed tempera- 
ture distribution from the end to the mid-length of logs 4 and 8 ft 
long and having diameters of 1, 2, 3, 4, and 5 ft. Each set of 
curves is for a different heating period, in most cases varying in 
intervals of 5 to 10 hr, depending upon the diameter. 

Five temperature curves showing the temperature at the 
different distances from the circumference shown in Fig. 2 are 
given for each heating period. The top curve in Figs. 4 to 11, 
inclusive, shows the temperature in the zone or isothermal] region 


at a distance from the circumference equal to !/1) (a) where ais 
the radius. The second, third, fourth, and fifth curves from the 
top show the corresponding temperatures in the zones that are 
at distances of !/, (a), 1/2 (a), °/s (a), and a, respectively, from the 
circumference. In other words, if a is the radius of the log and r 
is the distance from the center or axis to the isothermal region at 
which the temperature is computed, then (a — r) is the distance 
from the surface to the point at which the temperature is given. 


The ratio (=) is therefore equal to 0.10, 0.25,'0.50, 0.60, 


and 1.0 for the curves shown. In the last case, r = 0; hence 
the curve for this distance from the circumference shows the 
temperature variation along the axis from the end to the mid- 
length of the timber. The distance from the circumference to the 
point on the radius at which the temperature is computed is 


evidently equal to the numerical value of (a2 multiplied by 


/ 


es) 


the radius. For example, if the radius is 12 in. and 
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0.10, the distance from the surface is 12(0.10) or 1.2 in. The 
curves are shown for only one half the length since the rate of 
temperature change is assumed to be the same from either end. 

Group 2. Figs. 12 to 14, inclusive, show the computed tem- 
perature distribution for various heating periods at all distances 
from the circumference to the center at the transverse mid-section 
of logs 4 and 8 ft long, for diameters of 1, 2,3,4,and5ft. Similar 
curves could be prepared for a transverse section at any distance 
from the end surface by plotting the temperature data from the 
curves shown in group 1 as was done in preparing the curves in 
group 2, The temperature at the transverse mid-section, how- 
ever, is generally of most interest, since the lowest temperatures 
are obtained in this region. The curves given in group 2 (Figs. 12 
to 14, inclusive) show the temperatures obtained at the end of 5- 
to 10-hr heating periods, except for particularly long heating 
periods or when the rate of temperature change was very slow. 
Under conditions in which this rate was unusually slow, somewhat 
longer time differences were used. 

The time required to obtain a given temperature (which is 
found between any two time-temperature curves) can be deter- 


mined with sufficient accuracy for practical purposes from the 
proportional distance on the vertical line intersecting the two 
curves between which the temperature is found. This assumes 
that the change in temperature between any two time-tempera- 
ture curves is directly proportional to the time heat is applied 
during this time interval. For example, assume the required 
temperature is 170 F at 8 in. from the surface of a log 36 in. diam 

— 8. or approximately 0.45. 


Ge 
18 


The vertical line from 0.45, Fig. 13B, intersects the horizontal 
line from 170 F at a point between the temperature curves for 
heating periods of 60 and 70 hr. The temperature after heating 
for 60 hr is 164 F, and that after heating for 70 hr is about 173 F, 
or an increase of 9 deg in 10 hr. When the proportional change 
in temperature is assumed to be the same as the proportional 
es ‘ . . 170— 164 

change in time of heating the temperature ratio is 173-164." 
6/, or */;. Since the time interval is 10 hr, the time is found from 
the relation, 2:3 = 2:10 and z (the time that is required to in- 


and 8ftlong. In this case 


JANUARY, 1946 


TRANSACTIONS OF THE A.S.M.E. 


WVIQ ‘NJ 09 NV DNOT LY § DOT JO ANG 


(4774) ONT WOYS IONVLSIG 
be al. Oe ae a CORT 2 yO. 


WVIq ‘NJ 09 aNV 
WOU] SHONVLSIC INGUAIIG, LV GANIVLAQ SHUALVUGANAT, TT “OI DNOT Li p DOT JO GN WOUT SAONVISIG INGUGIGIG LY GANIVLEO SANALVAGAWAT, OT “SIS 


(24774) ONF WO¥s FINYLSIC 


Z / 0 


G / 0 


s 
ass 09/ 
6 OL OLI 
OL/ i 
P 06 09/ 
@ ees 061 
$ ol 
/ 5 
card - Oe ; | = 002 ¢ 
SunoH ow GOWId NTH” = let Ke , ore 
€ 3 O/ 
s 
06/ ; i 0s! f 
* Be i ale im me SYNOH 00! 
& 2 JOlY7d ONILVIH 
} O/ j . 06/ 0S/ 
2 t ——— 1 2 
: 06/ SE =33 olz om 
i = 
Hi all SYNOH O11 SYNOH 06 CL 
T Ez Ole d0lAId INILVIH d0lek Id ONILVIH ali 
ISYNOH OF! OIF INILVIH [ ¢ =| ie 091 «2/2 
£ i = eae Se 
s A 06 06/ 
t= = SYNOH 08 
£ = 002 gold3d INILVIH 
e = 
a= S OW 
: iter = He (= 
+S OE 
— — 
+ 2 05) 
. Z | | = ONILYIH 
—— aur ie s aa 
v 
/ / / O6/ g £ } __}__ 
eal lille Ole ys Eau ea 
if ze [2 wale = ¢ 
SYNOH 08 SYNOH OL SYNOH 09 Leet 
COld Id INILVIH Old Id INILVIH Old QNILVIH | 
¢ aay 6 Sarat fase =e 
' 4 : a= 2 
v € 06 | 
£ + 
+ OW é 
[ [ Tei i] 
iS 2 | 
JEM a is oof 
Z emia li i Wi [ = 
OS / 
a \i - seesl ; i 
, I) 
/ / 
G6/ 
= 2/2 ele 
BEE rae Sara SYNCH OF SYNOH OF 
SYNOH OS SYNOH OF SYNOH OF 
GOldId INILY FH 01YFd INILVIH GOlYy7d INILVIH COlYId ONILVIH COlYId DNILVIH 


STHIM-09 YILIWVIC'*LIFS-@ HLINIT 


S YNOH OL SYNOH 09 
golydad QDNILVIH G0lYid INILVIH 


Dal 
Ss 
s 


06 


nr 
Ss 
~n 
(4) FYNLVYFIWIL 


ele oc/2 


SYNOH OF SYNOH 02 
O0ldId DNILVIH C0l¥sd DNILVAFH 


SFHINI-O9 HILIW VIC 'LIZI-& HLIONI7 


MacLEAN—RATE OF TEMPERATURE CHANGE IN SHORT-LENGTH ROUND TIMBERS 


WVIC ‘NI Sf INV 9€ GNY ONOT LY 8 INV $ SHOT ‘NOMOaS-a1 
40 UAINAQ OL GONGUAAWAONIO WOUT NOWMAamIsiq GUALvUadIWaAT, ET “OMT 


Sniavy Ag JIO/AIT 


Jovsens WO4S JINVLSIO~ Poe OILY 


or 60 90 L0 90 60 ¥0 £0 20 10 O 


snigvd AP GWIAIAIG 


JIWINS WOYS FINVISIO= %p OlYY F 


0) 60 90 LO 90 $0 #0 £0 20 10 O 


l 


@ NI-Gt YTLIWV/O 
1FF/9 HLINIT 


N/-9¢ YFLINVIC 
LIFTS & HLINIT 


Bi Ss 


Z 


| 


l 


> = 


L a 


cle 


+ + 


Bole 


bd 12 9! S/ ral 6 ] ¢ 0 
(SPHINI) JINIYFINN IID WOYS FINVLSIO 


g 


ah ah nee a 
oh Oe Rp Be, nee a 
(SIHINIIINTYTINNIAID WOYF FINYLSIO 


NS 


i 


i 
SPHINIGE SFTLINVIO 


SIHINISE YTLINVIC 
LFI4 & -HLINIT | 


Le LIFIF9 HLINIT 


II 


INN 


| 


+ elec ea Ge 


if i jmmezyzZ 


o S$! 2/ 6 9 £ 0 
(SHHIN) FINTATINNIGID WOYS FINVLSIC 


Factialcilimtieieimlon lla aL T 
8/ £4 2 6 
(SFHINI)FINTYISWNISII WOYF JINVLSIG 


OL 
08 
06 

00/ 


oll 


OL 


06 


9 £ oO 


(A) ZYNLVATIW IL 


(se) FYNLVATIWIL 


WVIC ‘NI $Z INV GI GNV ONO] LY 8 INV 7 SHOT ‘NOIMLOAS 
-C1JN 10 UTLNAD OL GONAUTAWADNID WONT NOMAAGIUISIG ATALVARAWAL, ST “OMT 


Snidvy AG _AIA/A/T 


sniagved AF JIAG/AIO 
FIVIONS WOAY TINVLSID=72% OILYY 2 


72vd9Ns W044 FINVLSIO =m Ollvd 


or 60 90 LO 90 50 ¥0 £0 20 10 O or 60 90 LO 90 S50 #0 £0 20 10 an 
AL F 
STHINKRZ YFILIWVIO SIHIN-G2 
tO IES EL! {| YILIWVITG. 06 
a LS%-HLINIT 
4 + T ool 
| oll 
S 
L O2l 
nN 
= OTR 
ic m 
Or! ® 
g 
= MEGS 
ons 
& 
| oLr~ 
=| Se { 
Wore Og! 
a 06/ 
}— 1 002 
| | alle 
Ze ea j 2/2 ee 
eae plete 
a ee ae a a aioe ls he ae ielical 
Zz 6 9 € 0 2/ 9 € 0 


6 
(S3H9W/) FINTATILNIGID WOY4F FINVLSIO 


(s9H9M/) FINTUFIWNIGII WOAS FINVLSIG 


SFHIN/-2! 
LIFI & -HLIN. 


iS Y0n og 
T T aa 


+ = 
a I 
9 & ¢ € Z / 0 9 Si e 5 z / 7) 
( $I HON) FJINTATIWNINMI W044 FINVLSID (s FHINI) FJINIYTSNNIY/I WOYS FINVLSID 


(J+) JUNLYYIdWIL 


10 


DISTANCE FROM CIRCUMFERENCE(INCHES) 
ONES SIGE ota 2.8/5; 


ad 
SEI Se Ae habe 
——t— — 2/ (Zo 


LENGTH 4-F EET 
DIAMETER 60-7NCHES 


ee 
Le 
Deer aesccssc2777 


NNSS 
WAAAY 


TEMPERATURE (¥) 


0 : 
0 O1 02 03 04 05 06 07 08 O09 


RATIO "= DISTANCE FROM SURFACE 
DIVIDED BY RADIUS 


10 


\ 


TRANSACTIONS OF THE A.S.M.E. 


VOMCL ELA CHL ESO, 


JANUARY, 1946 


DISTANCE FROM CIRCUMFERENCE (/NCHES) 
l2 15 


On eae O an, eh PE LR TAP Ly) 


LENGTH 8-FEET 
DIAMETER 60-7NC 


SS 
I RE 


02 03 04 O05 06 07 08 09 


RATIO S" =DISTANCE FROM SURFACE. 
DIVIDED BY RADIUS 


Fig. 14 Trmprrature DistrinuTIon From CIRCUMFERENCE TO CENTER OF Mip-Section; Loags 4 anv 8 Fr Lone anp 60 In. Diam 


crease the temperature from 164 F to 170 F) is (2/3) 10 or 6 2/s hr. 
The total time necessary to reach a temperature of 170 F is then 
60 + 6.7 = 66.7 hr. 

The same method would be used in finding the corresponding 
temperature when the heating period is known. If, for example, 
the heating period is 23.5 hr, and it is desired to find from Fig. 
12D the temperature obtained at a point in the log half way 
between the axis and circumference, the figure shows that after 
20 hr heating the temperature is 142 F, and after 30 hr heating 
it is 164 F, or an increase of 22 deg in 10 hr. The temperature 
increase above 142 F is then computed by the proportion (23.5 — 
20):10 2:22, from which it is found that the temperature 


“ 22 = 7.7 deg F. The temperature obtained 
after heating for 23.5 hr is then obtained by adding 7.7 deg to 142 
F and is approximately 150 F. 

If desired, the temperatures can be plotted against heating 
periods for any particular point in the timber, as for example at 
the center. Such curves would furnish a convenient method for 
reading directly the temperature obtained in a given time, or the 
corresponding heating period required to reach a definite tem- 
perature. 

Group 3, In order to obtain curves that would show the 
temperature at various distances from the circumference at the 
mid-length of logs of any diameter between 12 and 60 in., data 
used in plotting the curves in group 2 were also used in preparing 
the curves shown in Figs. 15 and 16. Each curve in Figs. 15 and 
16 was prepared by plotting the temperatures for each heating 
period for the logs of different diameters. Curves were then 
drawn through these plotted points for the same heating period. 
For heating periods over 50 hr, data were plotted only for timbers 
36 in. or more in diameter. 

Temperatures plotted in group 3, Figs. 15 and 16, are for 
points at distances from the circumferential surface of 1/4, 1/2, and 


change x 


oo 


3/, the radius, and for the center or where 1, When the 


required heating period or temperature is found between two 
time-temperature curves, the heating period or temperature can 
be found by using the proportional relation as explained for the 
curves in group 2. ¥ 

The data show that it is not feasible to heat the large-diameter 
timbers at a considerable distance from the surface because of the 
extremely long heating periods required. For example, under the 
conditions assumed in making the temperature computations, — 
Fig. 15B shows that if a temperature of 150 F is desired at a dis- 
tance of 12 in. from the circumference of a log 4 ft diam, the heat- 
ing period would be nearly 58 hr or 2.4 days. In such cases it 
would probably be better to heat to the minimum required tem- 
perature at a shorter distance from the circumference. After 
cutting to this depth the log could be reheated and with the re- 
duced diameter the rate of heating would be much faster. This 
procedure would also avoid overheating the wood near the surface 
and at the ends. 

The temperature curves shown in group 3 will be particularly 
useful for general purposes, since they show the relation of tem- 
perature and heating period for all diameters between 12 and 60 
in. If itis desired to find the temperature at distances from the 


—* = 0,25, 0.50, 


: 5 ‘ a 
circumference that are intermediate between 


a 
0.75, and 1.0, which are the distances shown in this group of 
curves, the data for any distance required may be found in the 
figures for group 2. These temperatures can then be plotted 
for the different diameter timbers as was done in preparing the 
curves for group 3. Temperatures for these four points can also 
be plotted against the proportional distance from the surface (as 
in Figs. 12 to 14, inclusive) for any diameter and any heating ~ 
period. Smooth curves drawn through these plotted points will 
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then show the approximate temperature for all distances from the 
surface to the center or axis of the timber. 

Method of Using Fig. 3. Fig. 3 is a chart that furnishes a con- 
venient method for finding the corresponding wood temperature 
when the heating conditions are different from those assumed in 
computing the temperature curves shown in the different groups. 

The initial wood temperature U, and the temperature U, 
obtained in the wood after any given heating period are read on 
the left vertical scale, while the heating-medium temperature U, 
is read on the right vertical scale. Values of U, the corresponding 
computed wood temperatures found Bee the figures, are given 
on the bottom scale. 

In using Fig. 3, it is merely necessary to place a straightedge 
connecting the initial wood temperature U, on the left scale with 
the heating-medium temperature U, on the right scale. This is 
illustrated by the line connecting the initial wood temperature 
U, = 50 F on the left scale with the heating temperature U, = 
180 F on the right scale. Assume, for example, that with the 
foregoing temperature conditions it is desired to find the corre- 
sponding computed temperature U (see Appendix for definition of 


symbols) at a point in the timber when the required temperature 
U,is 148 F. This is found as follows: 

Pass horizontally from 148 F on the left scale to the point of 
intersection with the line joining 50 and 180 F. Directly below 
this point of intersection. U is found on the horizontal scale to be 
177 F. The procedure is shown by the dotted line from 148 F. 
Again, let it be assumed that data from one of the charts shows 
that a temperature U of 133 F is obtained after heating for a 
period of 7’ hr, and it is desired to find the temperature U, that 
would be obtained at the same point, in the same time, and in the 
same timber when the initial wood temperature U, is 50 F and 
the heating-medium temperature U, is 180 F, instead of 70 and 
212 F, the temperatures assumed in computing the data used in 
plotting the curves. 

The temperature U, is found in this case by passing vertically 
from U = 133 F on the horizontal scale to the line connecting 50 
and 180 F. Directly across on the left scale U, is found to be 
108 F. These steps are shown in Fig. 3 by a dotted line connect- 
ing the temperature U = 133 F, and U;, = 108F. This, is sim- 
ply reversing the procedure for finding Uwhen U; is known. 
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Fig. 3 is convenient to use in finding how the wood tempera- 
tures vary when different heating conditions are employed. For 
example, in the illustrations given, heating under the lower 
temperatures assumed for U, and U, the wood temperature will 
be (177 — 148) or 29 deg F lower at the point under consideration, 
than if the initial wood temperature had been 70 F, and the heat 
ing-medium temperature had been 212 F, which were the tem- 
peratures assumed in computing the curves. Similarly when any 
point is at a temperature U of 133 F (with the temperature con- 
ditions assumed in computing the curves), when heating under the 
lower temperature with U, = 50 F and U, = 180 F, the same 
point would be at a temperature of 108 F or 25 deg F lower. 
This chart can be used for finding either U, or U instead of com- 
puting these temperatures from Equations [B,] and [C,] in the 
Appendix. : 

In general, problems relating to temperature changes will be 
either: 

1 To find the time 7, required to obtain a specified tempera- 
ture U, at some particular point in a timber; or 

2 To determine the temperature U, that will be obtained 


when the time 7, and temperature U are known. In any prob- 
lem the following variables are involved: Variables found by 
means of the temperature charts are U and 7, and variables 
relating to the timber in question are U,, U;, Uz, and Ty. The 
relation between T and 7), is shown by Equation [A] (Appendix), 
which is 74h,? = T(0.00027) - 

When U,, U,, and U, are assumed, U is found from Equation 
[C,] (Appendix) or from Fig. 3. The corresponding time 7 
required to obtain the temperature U can then be found from the 
proper chart, and the time 7’, can be determined from Equation 
[A]. 

When 7’,, Uz, and U, are given, 7 is found from Equation [A], 
and the corresponding temperature U can be found from the 
proper chart. The temperature U, can then be determined by 
means of Equation [B,] (Appendix) or from Fig. 3. 


ExampLes SHOWING How Cuarts Arr Usep? 


Example 1. Data are as follows: Diameter of log = 22 in.; 


3 Refer to Appendix for definitions of symbols used. 
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length of log = 8 ft; species, yellow birch; heating medium, 
water; water temperature during heating period = 180 F = U,; 
and initial wood temperature = 50 F = U,. 

Required to find the heating period (7) needed to obtain a 
temperature of 140 F (U;) halfway between the circumference 


There- 


— 1 
and center at the mid-length. In this case eae 2 
a 


fore, for a log of this diameter, the distance from the surface = 
Ve (1) = 5.5 in: 

From Table 1 the average specific gravity of yellow birch 
(based on weight when oven-dry and volume when green) is 
found to be 0.55, and the diffusivity h,? found in Fig. 1 for this 
specific gravity is about 0.000249 when the heating medium is 
water. In this example U,, U;, and U;, are given, while U, T, 
and 7, are to be determined. 

The first step is to find U, using Equation [C,] (Appendix), or 
Fig. 3. From Equation [C;] 


(180 — 140) (142) 
U = 212 | (180 — 50) | = 168 F 
This can be found readily from Fig. 3 by placing a straightedge 
connecting 50 F on the left scale with 180 F on the right scale, as 
shown by the straight line on this figure. Passing horizontally 
from 140 F on the left scale to this line the temperature U is 
found directly below on the lower scale to be 168 F, which is the 
temperature found by using Equation [C;]. 

The second step is to find from Fig. 16(B) the time 7 required 
to reach a temperature of 168 F in a timber 22 in. diam and 8 ft 
long. Passing from 168 F in Fig. 16(B) to a point directly over 
the diameter of 22 in. shown on the abscissa, it is found that 
this point is between the curves for heating periods of 20 and 30 
hr. The temperature after heating for 20 hr is found to be 153 F, 
and after heating for 30 hr is 174 F, or a temperature difference of 
21 deg F. The proportion of the vertical distance between the 
two curves for the temperature of 168 F is then computed as 


168—153° 15 
174 — 158 = a and since the time interval is 10 hr, the time T is 


15 
computed as | 20 + 31 co) | or about 27.1 hr. 


This is the time that would be required if the diffusivity were 
the same as that assumed in computing the temperature curves; 
but since the diffusivity for yellow birch was found to be about 
0.000249, the time 7, must be computed from the relation shown 
by Equation [A] (Appendix) that gives 


_ 7(0.00027) 
* 0.000249 


or 29.4 hr, which is the time to be determined. 

If steam at 212 F were used as the heating medium, the diffu- 
sivity as indicated in Fig. 1, would be about 0.000276. With 
U, = 50 F, U, = 212 F, and the required temperature U; = 140 
F, the corresponding temperature U is determined from Equation 
[C,] (Appendix) or from Fig. 3 tobe 149 F. From Fig. 16(B), the 
time 7’ required to reach a temperature of 149 F in a log 22 in. 
diam is found to be about 19 hr, and the required heating period 

. 19(0.00027) 

, is 0.000276 18.6 hr. 

Example 2. With the same data as given for example 1 find 
the temperature that wouldbe obtained at!/, the distance from the 
circumference to the center when the heating period T, is 25 hr. 

@—7r 


= (27.1) (1.084) 


In this example, = 1/,, and the distance from the circumfer- 


ence = 1/,(11) = 23/4 in. instead of 5!/,in.asinexamplel. In 
this example U,, U,, and 7, are given, and 7, U, and U, are to be 
found. 


First, find 7 from Equation [A] (Appendix). In this case 


pes T',(hz”) 


fal __ 25(0.000249) 


kh?’ © ‘0.00027 


= 23 hr, approximately 


Next, find U for this heating period. From Fig. 16(A) the 
temperature U is found to be about 185 F, after heating for 23 hr. 

Then find the temperature U, from Equation [B,] (Appendix), 
or from Fig. 3. 

Using Equation [B,] 


180 — 5 ee 
spree E 80 Ee =] HW eA 


which is the temperature to be determined. This temperature 
U,, can also be found from Fig. 3 by placing a straightedge on 50 
F on the left scale and on 180 F, right scale. Reading up from 
185 F on the bottom scale to the line connecting 50 F and 180 F, 
and directly across to the left scale, the temperature U, is found 
to be about 155 F as computed from Equation [Bi] (Appendix). 

If steam at 212 F were used as the heating medium the diffu- 
sivity would be 0.000276 (as found for example 1), and 7’ would be 
25 (0.000276) 

0.00027 
period is found from Fig. 16(A) to be about 188 F. When U, = 
50 F, and U, = 212 F, the corresponding temperature U, is found 
from Equation [B,], or from Fig. 3, to be about 185 F. 

Example 3. Data are as follows: Diameter of log = 48 in.; 
length of log = 8 ft; species, Coast Douglas fir; heating medium, 
water; water temperature during heating period = 190 F = U;; 
and initial wood temperature = 60 F = U,. 

Required to find the temperature U, obtained at a distance of 
15 in. (Z) from the end surface and at a distance from the circum- 
ference equal to 1/, the radius when the log section has been heated 


for 30 hr (T;). I= 


= 25.6 hr. The temperature U for this heating 


r . 
=1/,; hence the distance 


In this example 


from the circumference =!/, of 24 or6in. Inthisexample U,, U,, 
and 7’, are given, and 7, U, and U, are to be determined. 

From Table 1 the average specific gravity of green Coast 
Douglas fir is found to be 0.45, and from Fig. 1 the corresponding 
diffusivity for heating in water is found to be about 0.000293 = 
hse 

The first step is to find the heating period 7’ using Equation 
[A] (Appendix) 


_ Tr(hz?) _ 30(0.000293) 


= = = 32.6 hr 
0.00027 0.00027 


Because of the higher diffusivity of the Douglas fir, in comparison 
with the diffusivity of 0.00027 used in computing the temperature 
(0.000293 — 0.00027) 100 = 8.5 per 
0.00027 

cent faster, and the temperature attained in 30 hr at the point 
under consideration would be the same (for the same heating 
conditions) as that reached after 32.6 hr heating when the wood 
has a diffusivity of 0.00027. 

The second step is to find the temperature U in Fig. 9, From 


curves, the Douglas fir heats 


this figure the temperature U is found to be 155 F (to: Z = 15in, 
ee Sa i) after a heating period of 30 hr, and to be 169 F 


when the heating period is 40 hr. After 30 hr heating, the 
temperature at the point given would then be computed as 


and 
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32.6 — 30 
E + (226—20) (169 — 15) | = 159 F approximately 


This is the temperature U, or the temperature that would be 
obtained if U, = Up = 70 Fand Us = U; = 212 F; but since U, 
is assumed as 60 F, and U, as 190 F, U, is determined by using 
Equation [B:] (Appendix), or from Fig. 3, and is found to be 
about 141 F which is the temperature to be determined. 

This is similar to example 2, except that the point under con- 
sideration is between the end surface and the transverse mid- 
section, and the data are obtained from temperature curves 
shown in group 1. 

Example 4. With the data given in example 3 find the heating 
period 7, required to obtain a temperature of U, = 160 F for the 
same point in the log. In this example U,, U;, and U, are given, 
and U, T, and T, are to be determined. 

The first step is to find the corresponding temperature U when 
U, = 60F, U, = 190 F, and U, = 160 F. From Equation [C,] 
(Appendix), or from Fig. 3, U is found to be about 179 F. 

Fig. 9 shows that when the distance from the end = 15 in. 
and the distance from the circumference (a—r) = 6 in. or 1/4 the 
radius, the temperature of 179 F is between the temperature 
curves for heating periods of 50 and 60 hr. After heating for 50 
hr, the temperature is 177.5 F, and after heating for 60 hr, it is 
about 184 F. The heating period required to reach 179 F would 
then be computed as 


179 — 177.5 
Se | = 52. 
(2 = 1715 Sites 50 LEONE 
Since this is the time 7 required when the diffusivity is 0.00027 
the time 7; is as follows 


_ 52.3(0.00027) 


= = 48.2 
0.000293 oe 


This is similar to example 1 except that the point at which the 
temperature is to be determined is between the end surface and 
transverse mid-section. 
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Appendix 


MeErTHop oF CompuTING TEMPERATURES IN RouND TIMBERS OF 
Sxort Leneta WHERE Enp HEATING as WELL AS HEATING IN 
THE RapraL Direction Must Br Taxen Into ConsipERATION 


In order to derive an equation for computing the temperature 
U at any point r inches from the axis and Z inches from one of the 
end surfaces, it is necessary to solve the fallowing partial differ- 


ential equation 

6U &eU 1 6U 

= = Pb cde oe niger. 

6T E {2242 a, 
In this equation h? is the diffusivity factor for heating in the radial 
direction and g? is the diffusivity factor for heating in the longi- 
tudinal direction. 

Equation [1] must be solved, subject to the following boundary 
conditions: Let @ = (U — U,), where U; is the temperature of 
the heating medium and U is the temperature at the point under 
consideration at any time 7’, 

Let Uo = initial wood temperature; a = average radius of log, 


and LZ = length in inches. 
Ife = (U— U,), U = Up and @ = (Uy) — U;) when T = 0 


6 = Owhenr =a 
Yh, = (V) 
Z=L 


and T > 0 


6 = f(r,Z) when T = 0 
Substituting 6 in Equation [1] 


Ns, Be ate 
Sa 6r2 or Or T5773 || ame 


A solution of differential Equation [la] to meet the boundary 
conditions gives 


m=o n=O = 
» > =P hiner ee Z 
= Am,ne ( Lt ) sin os Jo(A,7). [2] 
m=1 n=1 


where 


4 s ot mary 
SSS BN i d 
Am,n i a f fru) Io ar) sin Zw 


When f(r,x) is a constant and equal to (U» — Ui), as in this case 


L 
A Te = Uadien fe i mau 

= d Jor, === 
Am,n ela.) Jo ip ‘ rJo(A,7) sin LT 


_ f=) L)\ 1 cos me 
* [oo 28 | (4) Givin ob eo 


oe Re es odd, and = 0 when mis even 
MarnAd (Apa) 


Substituting the value 


8(Uo — Ui) 
mar ,,aJ 1 (A,@) 
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for Am,n and the temperature difference (U — U,) for 6, Equation 
[2] may be written 


mM=o n=o@ 


5 Oe sD Daa a aa ) 


m=1 n=1 


= (hen? + a) 


. mrZe 
sin re Jonre 


In the foregoing equations 7’ is the time in seconds, Jo(Anr) and 
Ji(Ana@) are Bessel functions of the zeroth and first order, respec- 
tively, and (A,r) = x in Jo(x) and Ji(z).4 Term X, is a root of 
Jo(\,@) = 0 and eis the base of the Napierian logarithms. 

Writing the first few terms of Equation [3] 


1 Z 
U=U+ (W— oy §| ant L 
Jo(rrje 
972q? 
z —T (AA + 
I : sin — on Jo(Mrje ( 2 ) 


Madi (ia) 3 L 


1 1 
a Arad 1 (Ara) 5 sin 


242 2572 
roe (ane + Se") ai ee 


-T (iors? + t) 
+ sae jer = Jalna 


gees 1 1. 4 ork 
d2aJ1 (Aza) ae ot, 


as 
—T (hd? + 
Jo(rarje ( : *) + Bes 


1 -T (arse + 7) 
+ ———_ er x L 
a O5a) Hose 
Qn2q2 
1 1 Z —T ( h?3? + 
+ —— = sin — Be Jo(dar)e ( L? ) +... 


Aad, 1 (Asa) 3 L 


1 tt (irre 
deuerioua) 7, ole ( 7) 
1 1 37 -?T (ire a ne 


Tz Zolure 


may 


The number of terms that must be used in making temperature 
calculations will depend to a considerable extent upon the value 
of the negative exponent of the Napierian log base e. Since this 
exponent involves the time T and the dimensions of the timber (as 
well as the diffusivity factors h? and qg?), the series converges more 
rapidly as the time T increases, and as the diameter and length of 
the timber decrease, since the values of \ vary inversely with the 
radius, and the second term of the negative exponent varies in- 
versely as the square of the length L. 


oR had; 1(Asa) 3° 


4 Tables of Bessel’s functions for different values of z and roots of 
J,(z) =0 are given in ‘‘Smithsonian Physical Tables,’’ Smithsonian 
Institute, Washington, D. C., and in ‘‘Fourier Series and Spherical 
oe by W. E. Byerly, Ginn and Co., 1893, New York, 


Mertsop, or Usine Ficures Wuen Dirrusiviry, HEATING 
ConpITIons, oR Botu, ARE DirrERENT FRom THOSE ASSUMED IN 
CoMPUTING THE TIMB- TEMPERATURE CURVES 


Equation [A] which follows, shows the relation of time and 
diffusivity, and Equations [B] and [C]'show the method of finding 
the corresponding wood temperatures when the heating condi- 
tions are different from those assumed in computing the data for 
the temperature curves shown in the figures. 

Symbols applying to computed data in figures are as follows: 


Up = initial wood temperature, assumed as 70 F 

U, = heating-medium temperature, assumed as 212 F 

h? = diffusivity in radial direction = 0.00027 sq in. per sec 

q? = diffusivity in longitudinal direction = 2.52(h?) = 0.00068 

U = temperature obtained in given time T at some particular 
distance from surface. Values of U are temperatures 
within wood at distance from surface shown by figure 

T = time of heating (in hours) necessary to obtain tempera- 
ture U. This period is found from figure for diameter 
and length under consideration 


Symbols, applying to temperature computations when the 
diffusivity, heating conditions, or both, are different from those 
used in computing the data plotted in the figures, are as follows: 


U, = any initial wood temperature 
U, = any heating-medium temperature 
h,? = any value of diffusivity for transverse direction, sq in. 


per sec 


Since the diffusivity in the longitudinal direction (q?) is taken as 
approximately 2!/. times the diffusivity in the transverse direc- 
tion for the calculated temperatures, only values of h,? need to be 
considered. 


U, = temperature obtained at definite point in timber when 
initial wood temperature U,, heating temperature Us, 
or both, are different from Up = 70 F, and U; = 212 F, 
the initial wood and heating-medium temperatures 
assumed in computing curves shown in figure 

T, = time of heating (in hours) required to obtain temperature 
U,, when either diffusivity h,?, initial wood tempera- 
ture U,, the heating-medium temperature U,, or all 
three variables are different from corresponding values 
used in computing data for curves 


Relation of Heating Period and Diffusivity. The relation of the 
heating period 7’ and the heating period 7, is shown by the 
equation 

Dee) NE (O100022)) asthe aero ats « <fere ae [A] 


In this equation 7’, is the heating period assumed or to be 
determined when the wood has a diffusivity h,?, and T is the heat- 
ing period found from the figure for the timber in question. 
Values of h,” for green wood of any given specific gravity can be 
found in Fig. 1, when either water or steam is used as the heating 
medium. Specific-gravity values for different woods, based on 
weight when oven-dry and volume when green, are given in 
Table 1. 

Relation of Initial Wood Temperature, Heating-Medium Tem- 
perature, and Temperatures Obtained in Wood. A simple propor- 
tional relation exists between the temperatures obtained under 
any given heating conditions. For example, the relation of the 
temperatures obtained when the initial wood temperature U,, 
the heating-medium temperature U;, or both, are different from 
those assumed in making temperature computations are ex- 
pressed by the equation 


(Us) (Ui) 
(Uz =a U,) (UV, = Uo) 
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When U,, U;, and U are known it is necessary to deteymine U,, 
the wood temperature obtained when U,, U,, or both are different 
from Ujand U;. Solving for U, in terms of U;, U,, U1, Uo, and U 


gives 
_,, [=U — UV) 
U, =U, | (U,— Us | Wehweeices: [B] 


When U,, U,, and U, are known it is necessary to solve for U to 
find the heating period T from the proper figure. Solving for U 
in terms of U,, U;, U,, Uo, and U, gives 


o» © one [Op U,) Ci Ue) 
Cas | (aut) | Oot [C) 
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If the values of Up) = 70 F and U; = 212 F (the initial wood and 
heating-medium temperatures assumed in making the computa- 
tions) are substituted in Equations [B] and [C] 


U, =U,— [& — oe = ay ae [Bil 


U = 212 ee Ls aan 


(U, “im U,) 


Values of U, and U can also be determined by means of Fig. 3, 
as explained in the text. 


(Owing to travel emergency conditions existing when this paper was presented, discussion will be accepted until February 10, 1946) 


Flame-Temperature Measurements in 


Internal-Combustion Engines 


By O. A. UYEHARA,! P. S. MYERS,? K. M. WATSON,? ano L. A. WILSON! 


The theory, development, and calibration of an electro- 
optical pyrometer, capable of determining instantaneous 
true temperatures of luminous flames, are described. A 
description is also given of auxiliary apparatus used with 
the electro-optical pyrometer to obtain simultaneous 
values of pressure, temperature, and crank position in an 
internal-combustion engine. Typical data obtained with 
this instrumentation as applied to a Diesel engine are 
included. 


HE need for instantaneous and continuous temperature 

data throughout the combustion process in internal- 

combustion engines has long been recognized. Combus- 
tion rates, chemical equilibria, fuel vaporization, ignition lag, and 
heat transfer are all affected by temperature. Measurement of 
high temperatures in internal-combustion engines is complicated 
by the rapidity of temperature fluctuations and by variations in 
flame thickness and emissivity. 

About two years ago the chemical-engineering and mechanical- 
engineering departments of the University of Wisconsin initiated 
a co-operative Diesel-fuel research program, the objectives of 
which could best be achieved if precise combustion temperature 
and pressure data were obtainable. A study of methods used by 
previous investigators to measure combustion temperatures 
showed that thermocouples (1, 2, 3),5 the sodium-line reversal 
method (4, 5, 6), the total-radiation method (7 8), and the “two- 
wire’ method (9) are obviously not adaptable to instantaneous 
and continuous measurements. The “additive” method, as 
reported by Graff (10), is instantaneous, but it is useless when 
radiation from carbon particles becomes appreciable; and, in 
addition, the transformation of radiation to temperature is 
laborious and involves the use of some questionable assumptions. 

Hottel and Broughton (11) employed an optical pyrometer 
with filters, using two different wave lengths of radiation to 
measure the temperature of a steady luminous flame, and 
established the validity of the two-color principle. Although the 
Hottel-Broughton pyrometer was not directly applicable to the 
proposed combustion-engine research, the two-color principle 
appeared to be feasible. This led to the development of an 
electro-optical pyrometer which is instantaneous in response to 
flame-temperature variations and will determine the temperature 
at every point in the combustion process for an individual cycle 
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of the engine. Auxiliary instruments were also developed for the 
measurement of combustion-chamber pressures, injection-nozzle 
opening, and the indication of completion of combustion, all in 
relation to crank position. It is believed that this instrumenta- 
tion may prove useful in a variety of applications, particularly 
where rapidly fluctuating high temperatures and pressures are 
encountered. 


Tue Evecrro-OpricaL PYROMETER 


The relationship between the intensity of monochromatic 
radiation from a black body, the wave length of the radiation, 
and the temperature of the black body is expressed by Wien’s law 


—C2 
oe Sal ee ate ag ean se ott (1] 
where 
J = intensity of radiation at wave length 
C,, C2 = known constants 

\ = monochromatic wave length 

e = base for Naperian logarithms 

T = absolute black-body temperature 


A graphical representation of Equation [1] is shown in Fig. 1. 

For a nonblack body, it is necessary to introduce a new factor 
(Ey) called the monochromatic emissivity, which is defined as 
the ratio of the intensity of radiation from a nonblack body to 
that from a black body at the particular temperature and wave 
length involved. Thus, for a nonblack body 


es 
J = EC, °ed? 


The intensity of radiation from a nonblack body may also be 
represented mathematically by 


= 
: J = Cyn *PA 


where 7’, is the apparent temperature and is defined as the 
temperature at which a black body would have the same mono- 
chromatic intensity of radiation as the nonblack body. 

Equating [2] and [3] 


Ze Ge = 
iD es CN RE ie AR eee eS [4] 
or 
Calas 1 
ea (eee ra 4 
In Ey a (2 1) [4a] 


This is the form of equation that is generally used to correct 
optical pyrometers for departure from black-body conditions 
where the monochromatic emissivity is known. ‘These relations 
are shown graphically in Fig. 1, where the dotted line represents 
a nonblack body at 4000 deg R. The apparent temperature at 
4 would be 3000 deg R, and at »g 2500 deg R. If the mono- 
chromatic emissivity HE), were the same at all wave lengths, 
Equation [4a] might be written for each of two different wave 
lengths, \4 and Ag, and the monochromatic emissivities equated. 
Then, by determining the apparent temperatures at \4 and dz, 
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SOLID LINES REPRESENT SPECTRAL) 
DISTRIBUTION OF RADIATION 
FROM BLACK BODY AT GIVEN 
TEMPERATURE. DOTTED LINES 


INDICATE RADIATION FROM A 
NON-BLACK BODY AT 4000°R 


J-INTENSITY OF RADIATION 


> -WAVE LENGTH OF RADIATION 
Fie. 1 GrapxicaAL REPRESENTATION OF WIEN’s Law 


the true temperature of a nonblack body might be calculated 
without knowing its emissivity. 

In general; monochromatic emissivities are not equal at 
’ different wave lengths. It has been shown by Hottel and 
Broughton (11) that the monochromatic emissivity of a luminous 
flame is expressed by an equation of the form 


—KL 
TN ol ate Oe eee ee [5] 
where 
K = absorption constant per unit flame thickness and is inde- 


pendent of wave length 

L = flame thickness 
= constant for limited wave-length range 

It was also shown by Hottel and Broughton (11) that the true 
temperature of a flame could be determined by using this relation- 
ship between the monochromatic emissivities and measuring the 
apparent temperature at the two wave lengths involved. Fol- 
lowing the development of Hottel and Broughton, substitution of 
Equation [5] into Equation [4a] gives 


—KL 
Sabinay ;) eaeee 2 = ly 4) 
In (: e ) SNe T, 7 [6] 


Writing Equation [6] in exponential form and rearranging 


= ee (r,—7) ee [7] 


Taking the logarithm of each side and rearranging 


R 
I 


—KL 
eae 


—C2 1 seg 
=Kin= Wein tee ™ NTA 7) | pea Sieg [8] 


Two such equations may be written for two different wave 
lengths, denoted by the subscripts 1 and 2. Equating the two 
KL terms, which are equal over a limited wave-length range, and 
rearranging, gives the following 


(A1) 1 


== Ca (aul 1 == Cau a 
E aoa 7a 7) | = Fi asm Ce 7) | [9] 


Tf Ay, Az, on, and a2 are known, and the two apparent temperatures 
are measured, the true flame temperature T is the only unknown 
in Equation [9]. Once the true temperature is determined the 
term KL may be evaluated, thus giving an indication of the flame 
emissivity. Thus true flame temperatures may be obtained 
from two optical pyrometers which simultaneously record 
apparent temperatures corresponding to different wave lengths of 
radiation. 
Apparatus. 


In the use of this method for measuring true flame 
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temperatures, substantially monochromatic radiation of two 
different wave lengths must be obtained. Although filters have 
generally been used, a system of lenses and prisms was adopted in 
this work. This method does not require a division of the 
radiation as would the use of filters. Furthermore, the lenses and 
prisms are not affected by temperature and humidity changes. 
To obtain sufficient dispersion two prisms were used. 

A schematic diagram of the electro-optical pyrometer is shown 
in Fig.2. The radiation from the combustion chamber is focused 
by a lens system on a limiting diaphragm which has an opening 
somewhat smaller than the image of the source of radiation. 
Thus vibration or other small movements of the source of radia- 
tion do not cause a variation in the light permitted to pass 
through the limiting diaphragm. After the radiation passes‘the 
diaphragm it is paralleled by a lens, dispersed by the two prisms, 
and again focused by another lens. A barrier containing two 
slits is situated at the focal point of this lens. A phototube is 
located behind each slit in such a position that it is actuated by 
the radiation passing through the slit. A vacuum-type photo- 
tube with a caesium-oxide-silver surface was chosen because of its 
instantaneous response, its electrical stability, and its spectral 
response curve. The outputs from the phototubes are amplified 
and impressed on the vertical plates of two single-beam cathode- 
ray tubes. For this particular application the cathode rays are 
displaced horizontally in proportion to the engine crank position. 
In this manner a graph of apparent temperature as ordinates and 
crank angle as abscissas is repeatedly traced on each cathode-ray 
screen. 

The amplifiers are of the direct-coupled type designed for use 
with vacuum-type phototubes and have a practically flat fre- 
quency response curve from 0 to approximately 15,000 cycles per 
sec. Two stages of amplification are required for the desired 
sensitivity. Batteries are used in the first-stage amplifiers to 
insure stability of operation. Voltage for the second-stage push- 
pull amplifiers is supplied from a voltage-regulated power pack. 


PHOTO 
TUBES 


PRISMS PHRAG! 
DIA ME oe COMBUSTION 
LENS 
’ 


Fie. 2 ScHematic Diaaram or Evectrro-OpricaL PyYROMETER 


Black-Body Calibration. Care was necessary to insure an 
accurate calibration of the electro-optical pyrometer since the 
method depends on the measurement of two apparent tem- 
peratures of nearly equal values and the use of this small tem- 
perature difference in calculating the true temperature. Two 
methods of calibration were considered; (a) calibration by the 
use of a broad-ribbon tungsten lamp with known emissivity 
values and, (b) calibration by sighting directly upon a black body. 

A search of the literature on the emissivity of tungsten (12, 13, 
14) revealed wide disagreement both as to the values to be used 
for emissivities and as to their variation with temperature and 
with wave length. Accordingly, the black body was used for 
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calibration purposes and the broad-ribbon tungsten lamp was 
retained as a secondary standard. Its use in this connection is 
described later. 

Black-body conditions were established by drilling a hole with 
a conical bottom in a graphite block. The graphite block was 
heated in the crucible of an induction furnace, and an optical 
pyrometer with a lamp calibrated by the Bureau of Standards 
was used to determine the black-body temperature. 

Above 2800 F, the formation of smoke from impurities in the 
graphite began to cause a departure from black-body conditions. 
In order to eliminate this difficulty the graphite block was 
blanketed with an atmosphere of inert nitrogen. This prevented 
the formation of smoke until a temperature of about 3100 F was 
reached, 

The procedure during the calibration process was as follows: 
The current to the induction furnace was adjusted and when the 
black body had reached a constant temperature a reading was 
then taken with the standard optical pyrometer. The electro- 
optical pyrometer was then placed in position viewing the black 
body. The corresponding cathode-ray deflections at the two 
wave lengths were measured by means of dividers and a steel 
scale graduated in hundredths of an inch. The electro-optical 
pyrometer was then removed from its position viewing the black 
body, and another reading of the black-body temperature was 
taken with the standard optical pyrometer. The average of the 
temperatures as measured by the standard optical pyrometer was 
taken as the correct temperature corresponding to the deflections 
noted. The foregoing procedure was repeated at various tem- 
peratures until the entire temperature range, representing zero to 
maximum deflection on the cathode-ray tubes had been covered. 
These data are shown in Fig. 3. : 

After the final calibration of the electro-optical pyrometer had 
been completed, it was used to calibrate the broad-ribbon tung- 
sten lamp which was to be employed as a secondary standard. 
A rigid holder for the tungsten lamp was necessary to obtain re- 
producible results since the temperature of the filament was not 
uniform over its various portions. The resulting graph of de- 
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flection versus current is shown in Fig. 4. Since these data are 
from several different runs, the reproducibility of results is 
clearly established. 

The original tungsten lamp proved defective and the filament 
broke after a short period of use. Consequently, it was necessary 
to repeat the black-body calibration at a few temperatures in 
order to standardize a second tungsten lamp. No effort was 
made to reproduce the original calibration and since the calibra- 
tion data plot as straight lines in Fig. 3, the entire temperature 
range was not covered. 

Determination of Equation Constants. In order to calculate 
true temperatures from Equation [9] it was necessary that \; and 
»: and the constants a1 and az be evaluated. 

Two factors had to be considered in choosing the wave lengths 
to be used. Approximately equal oscillograph deflections were 
desired for accuracy in scaling. Assuming equal electronic 
amplification, the deflection is a function of the spectral sensi- 
tivity of the phototube and the intensity of radiation at the 
temperature and wave length involved. Rassweiler and With- 
row have shown (15, 16) that in the range of wave lengths below 
0.5 micron a small amount of spectral energy from intermediate 
products of combustion is present and might cause a departure 
from the relationships established between the monochromatic 
emissivities of a luminous flame. There is also a number of 
strong water-vapor absorption and emission bands (17, 18, 19) 
from 0.190 micron to 0.36 micron, and from 0.95 micron to 314 
microns, and strong CO: bands (18, 19) in the region of 14 to 16 
microns. However, Rassweiler and Withrow (15, 16), Graff 
(10), and Landen and Blanc (20) found that black-body spectral 
distribution existed in the flames studied at wave lengths above 
0.6 micron. Accordingly, the phototubes were placed in the 
region of 0.6 to 1.0 micron and experimentally shifted in this 
region to obtain approximately equal responses under engine 
operating conditions. 

It was first attempted to determine the values of the wave 
lengths used by means of a sodium-vapor lamp and a mercury- 
vapor lamp. The locations of the sodium and mercury lines 
were easily identified in the visible region of the spectrum by 
inspection and by the response of the phototube. However, one 
of the phototubes was located in the infrared region of the spec- 
trum where it was necessary to identify the spectral lines by the 
response of the phototube alone. Difficulty was experienced in 
identifying positively the reference lines because of the fact that 
the response of the phototube varied considerably in the infrared 
region. Hence, the wave lengths determined by the phototube 
alone were considered unreliable. 

A second method of wave-length determination utilized data 
obtained during the black-body calibration of the electro-optical 
pyrometer. Assuming constant amplification from the elec- 
tronic equipment, the calibration data, when plotted as in Fig. 3, 


; ‘ ¢ 
should plot as a straight line of slope oe where C; is a known 


constant. From the slope of the calibration curve a value of 
0.696 micron was obtained for the shorter wave length. This is 
in close agreement with the value of 0.699 micron determined by 
thg use of the sodium and mercury lines. However, the value for 
the longer wave length of 0.902 micron was not in agreement 
with 0.860 micron estimated by the use of the mercury and 
sodium lines. Since the latter value was in doubt, and since the 
values obtained during the black-body calibration are “‘effective’”* 
wave lengths, the black-body values were taken as correct. 

The range and relative intensities of the wave lengths per- 
mitted to pass through each of the two barrier slits, Fig. 2, to the 
phototubes are of interest since the narrowness of the peaks of 
the relative-intensity curves are an indication of the mono- 
chromacity of the radiation reaching the phototubes. The 
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are shown graphically by the relative-intensity 
curves in Fig. 5. The accuracy of the equal-inten- 
sity assumption is discussed in the next section. 


CHECK RUNS ON DIFFERENT DAYS, 


The values to be used for a; and a» in Equation 
[9] depend on the wave lengths \; and». It was 


+ 1 


h 


suggested by Hottel and Broughton (11) that the 
range of wave lengths ordinarily used might be 
covered by two constant values of a. The values 


WwW 


recommended are 1.39 for the range from 0.3 to 0.8 
micron, and 0.95 for the range from 0.8 to 10 mi- 


cron. These values were obtained by plotting ex- 
perimental transmissivity data in such a manner 
that the slope of the line represented a. Inspec- 
tion of the data indicated that each of the recom- 


mended values corresponds to the average slope of a 
portion of the transmissivity curve. Inasmuch as 
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following method was used to construct the relative-intensity 
curves shown in Fig. 5: 

Data obtained during the wave-length determination were 
plotted showing the distances from a reference point to the points 
where known wave lengths were intercepted by the barrier. By 
measuring the distance from the two edges of the slit to the refer- 
ence point, and using the curve of Fig. 5, the range of wave 
lengths permitted to pass through each barrier slit was obtained. 

To obtain the relative intensity of the wave lengths passing 
through the barrier slits, the size of the image of the limiting dia- 
phragm was measured at the barrier with monochromatic radia- 
tion entering the electro-optical pyrometer. The wave length at 
the center of each slit was determined from Fig. 5, and that por- 
tion of the radiation at this wave length permitted to pats 
through the barrier slit was calculated from the width of the slit 
and the dimensions of the limiting-diaphragm image. This pro- 
cedure was then repeated for thirteen equidistant points across 
the width of each barrier slit. If the intensity of radiation is 
assumed to be the same for the wave lengths represented at the 
various points taken, the relative intensities for the wave lengths 
passing through each barrier slit are then given by the ratio of the 
respective portions of the monochromatic radiation passing 
through each barrier slit. The results of the analysis outlined 


the wave lengths employed in the electro optical 
pyrometer had been evaluated, it was considered 
more accurate to determine a values by measuring 
the slope of the curve at the points corresponding to 
these particular wave lengths. The values so de- 
termined are a; = 1.82 and a, = 1.05. The effect of any er- 
ror in these constants is discussed in the following section. 

Accuracy. Since this method requires very accurate measure- 
ments of two apparent temperatures, the effect of all possible 
experimental errors must be minimized. One possible source of 
error is the electronic equipment. Direct-coupled amplifiers 
with two stages of amplification, batteries, and a voltage-regu- 
lated power pack were used. At first, slight inconsistencies 
were observed, but after waxing resistors and connections, 
“burning-in” resistors, etc., reproducible results were obtained 
as shown by the close agreement of the data in Fig. 4. 

A small error’is possible in measuring deflections on the 
cathode-ray tubes since the spot is of finite size. Measurements 
taken by various individuals differed by a maximum of only 
0.02 in. on a 5-in. cathode-ray tube. This error was minimized 
by measuring each time from a given point on the spot, and by 
using the average of several readings. 

In measuring the black-body temperature with the standard 
optical pyrometer observations by various individuals differed 
by amaximum of 15 F. Maximum deviation of the final calibra- 
tion data from a straight line when plotted, Fig. 3,was only 10 F. ~ 

The effect of errors in evaluating a is discussed in reference 
(11), where calculations are presented showing that the error in 
true temperature due to a change in @ from 1.7 to 1.39 is only 13.7 
deg F for moderately thick flames. Thus a small error in deter- 
mining a@,will not seriously affect the calculated temperature. 

Forsythe (21) discusses the magnitude of the errors involved in 
assuming that the effective wave length of radiation used remains 
constant with changes in temperature, and he states that for most 
work the correction to be applied is negligible, amounting approxi- 
mately to 4 deg at 5400 deg R.. : 

The fundamental accuracy of this method of determining true 
flame temperatures has been demonstrated convincingly by 
Hottel and Broughton (11). However, in order to verify both 
the calibration and the calculations, a steady high-temperature 
flame was set up so that its temperature as determined by the 
electro-optical pyrometer might be compared to that measured 
by a standard method. 

In order that the required high temperatures might be reached 
a burner was constructed to use acetylene, with air and oxygen 
supplied in any desired proportions. The experimental flame 
was approximately 1/,in. X 1 in. in cross section. The burner 
was set up inside a water-cooled jacket in order to simulate the 
conditions in the engine. Two peepholes were provided through - 
which the flame might be observed. 
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The method used to determine the true temperature of the 
experimental flame was as follows: A broad-ribbon tungsten 
lamp of controllable brightness was placed behind the flame in the 
line of sight of the standard optical pyrometer. Due to the 
partial transparency of the flame, the standard pyrometer could 
see’ the tungsten lamp. The temperature of the tungsten 
lamp was adjusted until there was only a small difference 
between the reading of the standard optical pyrometer when 
sighted through the flame onto the tungsten lamp and 
when sighted directly on the tungsten lamp. Graphical inter- 
polation of the readings taken in this manner permitted 
estimation of the temperature at which the apparent tem- 
perature of the tungsten lamp and the true temperature of the 
flame were equal. 

The results of the determination of the true flame temperatures 
are given in Table 1. It was found that the flame temperature 
varied as much as 100 F over various portions of the flame 
visible through the peepholes; therefore, care was exercised to 
sight on the same portion of the flame with the standard optical 
pyrometer as with the electro-optical pyrometer. 

To determine the effect of flame thickness on the electro-optical 
pyrometer temperature readings, the burner was rotated 90 deg 
so that the flame thickness was approximately 1 in. The lens 
aperture opening in the electro-optical pyrometer was stopped 
down so that the oscillograph deflections were not more than full 
seale. The data in Table 2 were then obtained. 

The foregoing experimental evidence was accepted as satis- 
factory proof of the accuracy of the method and apparatus for the 
proposed internal-combustion-engine studies. 

Working Charts. Since the solution of Equation [9] is a tedious 
process, a chart was constructed to express true temperature as a 


function of cathode-ray deflections. This chart was derived 
from the following rearrangement of Equation [9] 
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In the solution of this equation values were selected for 7’ and for 
T ao and the equation was solved for T4;. By means of Fig. 3, 


TABLE 1 DETERMINATION OF TRUE FLAME TEMPERATURES 


I First Comparison 
Standard Optical Pyrometer (thin axis of flame): 
Temperature of lamp through flame (average of two readings), 


ODER AD Se eG Suite Sig, eM hte ercet CE hand adessoaiie, oeeaciane 4025 
Temporncers of lamp (average of two readings), degR........ 4025 
Electro-optical pyrometer readings (thin axis of flame): 
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Dehsetow AG Naiay sone cosines nate Gamer ere tenes iene eareM 2.84 


True LEMperature, GOSakeiy. cc 2 aiciercscsaetsle rcs eaters ra oO 4020 
II Sxrconp ComPaRIsoN 
Standard Optical Pyrometer (thin axis of flame): 
First trial: 
Temperature of lamp through flame (average of three read- 
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TemperacUure OL LAMP AEQiwW os cisecclewirc wee tins cule semiaiens 4145 
Second trial (lower lamp current): 

Temperature of lamp through flame (average of three read- 
POZE) Op Rome eter ical ctsnctsi creetarcte hctetmeenrc(e kisccleneee re 4025 
Tom peraturevomla mp; eg Riese sralelardan cups cue Hlewtereiere mua eferane 4005 
Using a graphical interpolation, the temperature of the 
lamp when sighted through the flame equals the temperature 
of the lamp at 4035 deg R 
Electro-Optical Pyrometer Readings (thin axis of flame): 
Deflection at i, in. 1.86 
Deflection at dz, in. oe aks 
True temperature, ROBE erie obo tie. aileveliecd erstotueleeee exe whens ledigiece 4020 


TABLE 2 ELECTRO-OPTICAL PYROMETER READINGS (THICK 
AXIS OF FLAME) 
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values of 74, and T'a2 were translated into inches of deflection on 
the oscillographs. The completed chart is shown in Fig. 6. 

Since a knowledge of the variations in the KL factor is of inter- 
est, a similar chart, Fig. 7, was constructed, showing lines of con- 
stant KL factor... 

Quartz Windows. Quartz was chosen as a suitable material to 
transmit combustion radiation because of the uniformity of its 
transmissivity over the range of wave lengths used in the electro- 
optical pyrometer (22). A polished quartz rod is fitted into a 
cylindrical steel holder which screws into one of the openings to 
the combustion chamber. To hold the rod against explosion 
pressures, a shoulder, forming a hole slightly smaller than 
the quartz rod, was machined at the outer end of the holder. The 
quartz was cemented in with litharge and glycerine to prevent 
leakage. 

Periodically it is necessary to remove carbon deposits from the 
inner end of the quartz rod because the vertical deflections of the 
cathode rays vary with the intensity of the radiation transmitted 
through the quartz window. However, since the presence of a 
coating of carbon has the same effect as a reduction in flame 
emissivity, the true temperature value obtained is affected only 
to the extent that the scaling accuracy is impaired by the reduc- 
tion in deflection. 


AuxttiaRY EQuipMENT Usrp WitH ELxcrro-OpticaL Py- 
ROMETER 


Sweep and Timing Arrangement. The two apparent tempera- 
tures to be used in calculating the instantaneous true temperature 
of a rapidly fluctuating flame must of necessity be measured at 
precisely the same instant. In the conventional internal-com- 
bustion engine this is best accomplished by an arrangement 
whereby the rotation of the engine controls the horizontal motion 
of the cathode rays and simultaneously provides timing marks at 
frequent intervals corresponding to known crank positions. In 
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addition, it should be readily adjustable so that any small portion 
of the combustion process can be expanded to reveal minute de- 
tails which are not clearly discernible when the entire process is 
shown on the cathode-ray screen, and it should function equally 
well at all engine speeds without alteration. These features were 
incorporated in the sweep and timing arrangement designed for 
use with the electro-optical pyrometer. } 

The sweep and timing arrangement was built as a separate unit 
and coupled to the end of the dynamometer shaft by an Oldham- 
type coupling. The characteristics of this type of coupling are 
such that the angular velocity of the driven shaft is the same at 
every instant as that of the driving shaft, even if the shafts are 
not in perfect alignment. 

The arrangement of the working parts of the unit is shown in 
Fig. 8. Secured to the shaft are two cam-shaped disks, H and J. 
One of these disks has a constant increase in radius per degree for 
330 deg, while the other has the same total increase in 90 deg. 
The remainder of the periphery in each case is at a constant 
radius until finally a radial decrease in diameter reduces the 
radius to the minimum value. Using either disk, light from a 
source A is projected through a radial slit C, behind which the 
edge of the rotating disk gradually cuts off the beam of light. A 
mirror # reflects the remaining light into a phototube D. Thus 
the amount of light reaching the phototube at any instant is 
directly related to the position of the disk at that instant. The 
phototube is connected to a direct-coupled two-stage amplifier, 
and the output voltage is impressed on the horizontal plates of all 
the cathode-ray tubes. Inasmuch as the voltage required to 
sweep the cathode rays a given distance is not the same for all 


tubes, potentiometers were installed as a means of making com- 


pensatory adjustments. Since the horizontal movement of the 
cathode rays is controlled by the variation in radius of a disk, the 
sweep is directly related to engine rotation during the period 
corresponding to K in Fig. 8. After a pause occasioned by the 
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constant-radius part of the disk, the abrupt decrease in radius 
instantaneously returns each ray to the opposite side of its screen 
to begin another sweep cycle. 

The portion of the combustion process viewed on the cathode- 
ray screen can be varied, either by changing from one disk to the 
other by switching lights A, or by varying the intensity of 
the light. Using the first method, the horizontal sweep of the 
cathode rays can be made to correspond either to 330 deg of 
crank rotation, or to 90 deg, depending on which of the lights A 
is illuminated. With equal ease, the total sweep of the cathode 
rays can be adjusted at will by varying the intensity of the light 
source A by means of a rheostat. These adjustments permit as 
much or as little of the combustion process as desired to be 
spread across the fluorescent screens of the cathode-ray tubes. 
The range, in degrees of crankshaft rotation, may be any amount 
from 9 to 330 deg. The 9-deg range is particularly advantageous 
where pressure and temperature changes occur most rapidly. 

The timing device is integral with the sweep arrangement. 
Using a milling-machine head, small holes were drilled in each 
disk on an arc 5 in. from the center of rotation and extending over 
that portion of the disk where the constant increase in radius 
occurs. The holes are 41/2 deg apart in the 90-deg disk, and 12 
deg apart in the 330-deg disk. The center hole of the series was 
drilled slightly larger than the others. 

A light source A’, a lens B’, the mirror #, and a phototube D’, 
Fig. 8, are so placed that light reaches the phototube inter- 
mittently through the holes in the disk as the disk rotates. The 
phototube is connected through amplifiers to the vertical plates 
of one of the cathode-ray tubes. This arrangement causes the 
cathode ray to be deflected vertically at intervals corresponding 
to the spacing of the holes in the disk. The position of each disk 
on the shaft was fixed in assembling so that the large hole passes 
light to the phototube when the engine is exactly at top center, 
causing a higher deflection which clearly marks the top-center 
position. 

In order to bring into view portions of the cycle outside the 
normal range, the holder for the lights, lenses, mirror, and photo- 
tubes was made movable about the axis of the shaft, and a ring 
graduated in degrees was mounted adjacent to the holder. 
When the holder is shifted from its normal position the long tim- 
ing mark on the screen no longer corresponds to the top-center 
position. However, since all the timing marks are shifted the 
same number of degrees as the holder, the position on the screen 
corresponding to top center can readily be found by means of the 
graduated ring. 


It was found necessary to balance the disk-and-shaft assembly 
dynamically and to mount the entire assembly on a heavy cast- 
iron block insulated from floor vibration. It was also found 
necessary to use bronze bearings rather than ball bearings since 
the latter, when tried, set up vibrations of a magnitude sufficient 
to cause noticeable distortion of the cathode-ray diagrams. The 
entire unit is painted black and inclosed in a lightproof housing 
to prevent stray light from reaching the two phototubes. The 
unit is shown in Fig. 9, ~ 

Pressure Indicator. An electro-optic pressure indicator de- 
veloped by Robertson (23) at the University of Wisconsin is em- 
ployed to obtain pressure data. A schematic diagram of the 
pressure pickup unit is shown in Fig. 10, in which light from source 
A is made parallel by lens B, then is reflected by a mirror C 
through a limiting diaphragm D onto a highly polished pressure 
diaphragm Z, which is screwed into an opening in the wall of the 
combustion chamber. Pressure variations in the combustion 
chamber cause the diaphragm # to change shape, dispersing 
the light reflected back by the polished surface. A portion of this 
light passes through perforations in the reflecting surface of the 
mirror C, and through a limiting diaphragm F to a phototube G. 
The output of the phototube is amplified and impressed on the 
vertical plates of a cathode-ray tube. The cathode ray is actu- 
ated horizontally by the sweep circuit previously described. 
Since the intensity of the light reaching the phototube is varied 
by changes in pressure, the cathode ray traces a pressure versus 
crank-angle diagram on the tube screen. 
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The electro-optic pressure indicator was calibrated in place on 
the engine under actual operating conditions by using a supple- 
mentary balanced-diaphragm indicator. This method eliminates 
the possibility of errors due to changes in the modulus of elas- 
ticity of the polished diaphragm, or in the reflectivity of its sur- 
face, when exposed to high-temperature conditions. The bal- 
anced-diaphragm indicator was similar to one developed by the 
Bureau of Standards, but was modified slightly to avoid any 
change in the engine-clearance volume when it was screwed into 
an opening normally oceupied by a quartz window. A tank of 
compressed nitrogen was connected to the outer side of the dia- 
phragm to balance the engine pressure, and a calibrated pressure 
gage was used to measure the applied pressure. When the 
engine pressure exceeded the nitrogen pressure, the diaphragm 
moved outward, making contact with the electrically insulated 
center portion of the indicator which was coupled electrically, by 
a condenser, to the grid of the first-stage amplifier used to amplify 
the output of the phototube of the electro-optic pressure indi- 
cator. Whenever this contact was made or broken, the cathode 
ray was temporarily displaced, producing a jog in the pressure 
diagram. A reproduction of a pressure diagram taken with the 
balanced-diaphragm indicator in operation is shown in Fig. 11. 
The effect of engine vibration on the condenser coupling is visible 
in this picture. Nevertheless, the point at which the diaphragm 
made contact is clearly indicated. By varying the pressure sup- 
plied to the outside of the balanced diaphragm and measuring the 
corresponding deflections, a calibration curve of deflection versus 
pressure was obtained as shown in Fig. 12. The intake pressure 
was taken as atmospheric since the error in this assumption was 
negligible. 

Since the magnitude of the deflections obtained depended on 
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the light intensity as well as the pressure variations, it was found 
necessary to supply closely regulated direct current to the light 
source. This was accomplished by the use of storage batteries 
which were constantly charged at a rate slightly higher than the 
current required for the light source. By the use of this method 
to obtain constant light intensity, results reproducible within the 
limits of measuring accuracy were obtained. 

Injection Indicator. The engine is equipped with a Bosch solid- 
injection system with the conventional commercial method of 
varying the time of injection. A shutter is attached to the nozzle 
plunger so that when the nozzle valve is lifted from its seat and 
injection started, the amount of shutter opening is increased. A 
light on one side of the shutter actuates a phototube on the other 
side whenever the shutter opening changes. The output of the 
phototube is amplified and impressed on the vertical plates of 
one of the cathode-ray tubes. This cathode ray is actuated hori- 
zontally by the sweep circuit previously described. Thus the 
beginning of injection is indicated by a vertical deflection of the 
cathode ray, and the end of injection is indicated by the return 
of the cathode ray to its original position. In addition, any 
pressure waves of sufficient magnitude to move the nozzle plunger 
and shutter are shown. 

End-of-Combustion Indicator. The end of combustion occurs 
either in the precombustion chamber or in the main chamber, 
depending on the load on the engine and the point at which 
injection occurs. One of the quartz windows is located in the 
cylinder head so that combustion in the main chamber can be 
viewed through it. By observation, it was found that at light 
engine loads combustion apparently does not carry over from the 
precombustion chamber into the main chamber. Therefore; at 
light loads, the end of combustion is indicated by the electro- 
optical pyrometer. Under heavy-load conditions combustion 
carries over into the main chamber. Accordingly, a phototube 
was placed in such a position that it is actuated by the radiation 
from the quartz window viewing the main chamber. The output 
of this phototube is amplified and impressed on the vertical plates 
of one of the cathode-ray tubes. This cathode ray is actuated 
horizontally by the sweep circuit. Thus the end of combustion 
can be determined regardless of where combustion ceases. 

Camera-Shutter Synchronizer. A permanent record of the data 
is obtained by photographing simultaneously all of the oscillo- 
graphs while all the devices which activate the cathode rays are 
in operation. In order that single cycles with their cyclic pres- — 
sure and temperature fluctuations can be studied, as well as 
general trends of numerous cycles, photographs of single cycles 
are desired. 

In order to photograph these single cycles, a shutter synchro- 
nizer is required which will open and close the shutter at the 
proper crank position regardless of engine speed, and which will do 
so one time only without being reset. 

The device, shown in Fig, 13, was developed to meet these re- 
quirements. A wooden disk J is coupled to the end of the sweep- 
and-timing mechanism shaft through reduction gears having a 
ratio of 40 to 1. Two small blocks A are attached to one side of 
the disk to form an adjustable cam. This cam actuates a fol- 
lower C which closes the breaker points D, permitting current to 
flow through the solenoid F# which opens the camera shutter. 
When the cam leaves the follower the breaker points open the 
solenoid circuit and the camera shutter closes. The combined 
length of the two blocks forming the adjustable cam corresponds 
to one engine cycle. Therefore the camera shutter remains open 
for one cycle, the point at which it opens and closes being deter- 
mined by the location of the adjustable blocks on the disk. A 
manually operated switch H in the solenoid circuit is closed when 
a picture is desired and opened immediately after the shutter 
closes. Thus one combustion cycle is photographed at a time, — 
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and as many additional photographs as desired can be taken by 
repeating the operation with fresh film. To prevent the inad- 
vertent closing of this switch when the breaker points are already 
closed, a light # is connected in parallel with the solenoid. 

In order to obtain compact and complete photographic records 
of original data, the cathode-ray tubes were placed in a vertical 
row beside a black title board’on which pertinent identifying data 
are shown by white letters. A view of this arrangement is shown 
in Fig. 14. The board is illuminated when taking photographs. 
Also, two small illuminated holes, 10 in. apart, to be used for align- 
ment and enlargement purposes, were provided adjacent to the 
bank of cathode-ray tubes and are photographed simultaneously 
with the oscillographs. Before photographs are taken, the hori- 
zontal travel of the cathode rays is adjusted by means of the 
potentiometers so that, at any time, a vertical line through all 
the cathode-ray ‘‘spots’” is parallel to the two illuminated holes. 
When the image of the oscillographs is projected onto a ground- 
glass screen with the proper enlargement, the alignment ‘‘dots’’ 
are 10 in. apart. The deflections corresponding to particular 
crank positions are determined by projecting lines from the 
corresponding timing marks parallel to the two dots. A 35-mm 
camera fitted with an f 1.9 lens is used to photograph the oscillo- 
graphs. 

Engine and Auziliary Apparatus. The engine to which the 
previously described apparatus is attached is a commercial-model 
Diesel engine, the specifications of which appear in Table 3. 
It has a high-turbulence precombustion-chamber type of head 


which has openings as shown in Fig. 15. The engine is coupled 
to an electric dynamometer. A view of the engine with the elec- 
tro-optical pyrometer, the pressure indicator, the injection indi- 
cator, and the end-of-combustion indicator in position is shown in 
Fig. 16. 

The length of time necessary for the engine to consume a, pre- 
determined amount of fuel, and the total revolutions during that 
time, are automatically measured by means of an electromagnetic 
tripping device. 

As part of the energy balance, the heat loss to the jacket cooling 
water is measured by means of calibrated thermometers, a direct- 
connected positive-displacement pump, and a calibrated #/s-in. 
sharp-edged stainless-steel orifice. 

In order to determine the average air-fuel ratio, the air con- 
sumption is measured by means of a thin-plate orifice constructed 
according to the A.S.M.E. Code. Air pulsations are smoothed 
out by volumes and pressure drops as recommended by the Code. 

Two smoke meters are attached to the exhaust line, one a con- 
tinuous-flow type and the other a sampling type. 

A Bureau of Mines gas-analysis apparatus is used to analyze 
the exhaust gases for oxygen, carbon dioxide, and carbon mon- 
oxide. The remainder is assumed to be nitrogen. 


PRELIMINARY RESULTS 


The engine head, as shown in Fig. 15, has a radial opening and 
a tangential opening in the precombustion chamber through 
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which temperature observations could be made. In order to 
determine which one would be best suited for temperature meas- 
urements, quartz windows were installed and observations were 
made through each of these openings. The procedure for ob- 
taining comparative data for the two openings was as follows: 
The engine was run at rated speed and load until equilibrium 
conditions were reached. Then the engine was momentarily 
stopped, the quartz windows cleaned and reinserted, and the 
engine started and allowed to run for about 5 min in order to again 
reach equilibrium conditions. Temperature observations were 
then made by sighting alternately the electro-optical pyrometer 
through the two windows, and photographing the oscillographs of 
apparent temperature variations. ‘ 

The temperatures determined through the tangential window 
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Fie. 16 Dissp, ENcGine anp AssociaTEp EQUIPMENT 
(1 Electro-optical pyrometer; 2 Pressure indicator; 3 Injection indicator; 4 End-of-combustion indicator; 5 Quartz window.) 


were slightly lower than those determined through the radial 
window, particularly during the first part of combustion. This 
was undoubtedly a result of the cooling effect of the combustion- 
chamber walls. During the last half of combustion the two 
temperatures were found to be nearly equal. 

At light loads, rapid variations in flame emissivity, with a 
consequent rapid variation in the apparent temperatures, were 
observed when the electro-optical pyrometer was sighted on the 
opening which views the precombustion chamber tangentially. 
At full load these variations in flame emissivity tended to de- 
crease both in rapidity and in magnitude. The flame-emissivity 
variations observed through the radial opening were always less 


rapid and of smaller magnitude than those viewed through the 


tangential opening. 

The analysis of the photographic records from the two openings 
failed to show any appreciable difference in ignition lag. 

After consideration of these factors, the quartz window for 
temperature measurements was left in the radial opening and the 
pressure indicator was inserted in the tangential opening, as 
indicated in Fig. 15. 

Having settled on a suitable location for the quartz window, 
and having previously established the reproducibility of tempera- 
tures measured with the electro-optical pyrometer, Fig. 4, the 
cyclic reproducibility of engine temperatures under constant 
operating conditions was next investigated. Some of the results 
of this investigation are shown in Fig. 17. These data show that 
the average trend is reproduced, but the temperature at any 
particular point in the cycle is not necessarily reproduced from 
cycle to cycle. Data taken on different days under the same 
operating conditions showed similar temperature fluctuations, 
but the average trend was the same as shown in Fig. 17. The 
pressure measurements showed closer agreement. 

As a preliminary to investigating the effects of fuel composition 
on the combustion process, tests were in progress at the time this 
paper was prepared to establish for a single fuel the effect of 
different operating conditions on the temperatures during the 


combustion process. Operating variables such as engine speed, 
engine load, point of injection, jacket-water temperature, etc., 
are individually varied while the other operating variables are 
maintained constant. 

For the preliminary tests reported in this paper, the operating 
variable was the point of injection, which was adjusted to occur at 
top center, at 18 deg before top center, and at 24 deg before top 
center, rated speed and load being maintained in each case. 
Photographs were taken with the horizontal scale expanded 
different amounts, so that particular portions of the combustion 
process might be studied in detail. The photographs were 
analyzed in the manner previously described. 
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Fig, 21 


OscILLoGRAM OF ENTIRE ComMBUSTION Process WITH 
InsEcTION AT Top CENTER 


Views taken during the run with injection occurring at 24 deg 
before top center are shown in Figs. 18, 19, and 20. In Fig. 18 
the horizontal scale is compressed so that the whole combustion 
process is shown. The timing marks are 12 deg apart, the high 
mark indicating top center. The effect of expanding the hori- 
zontal scale by increasing the light intensity passing through the 
slit is shown in Fig. 19, The timing marks are still 12 deg apart. 
The 90-deg disk was in use when the photograph reproduced in 
Fig. 20 was taken, consequently the timing marks in this picture 
are 4!/, deg apart. Fig. 21 shows the effect of late injection. 
The alignment dots are visible on the right-hand side of all the 
photographs. 

The true temperature, pressure, and KL factor fluctuations, as 
recorded on the illustrations Figs. 18 and 21, are shown in Fig. 22, 
together with data for a single cycle with injection occurring at 
18 deg before top center. Each of the dotted lines in Fig. 22 
was obtained by analyzing 10 cycles, similar to Figs. 18, 19, and 
20, for each point of injection and drawing average curves for 
these data. 

The data in Fig. 22 are presented at this time primarily to show 
that the electro-optical pyrometer is capable of following and re- 
cording the temperature variations inside the cylinder of an in- 
ternal-combustion engine. Not enough data have been obtained 
to explain the cyclic temperature fluctuations or the variations in 
temperature with point of injection. However, additional tests 
are now in progress to establish more completely the effect upon 
the temperature of varying the point of injection, as well as the 
effect of the other operating variables and the fuel composition. 
The results of these tests will be reported at some future time. 
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Discussion 


H. B. Norraaer.® As a tribute to experimental ingenuity and 
skillin tackling a difficult problem of practical interest, this paper 
is of outstanding merit. Even with slight imperfections, just so 
they would be consistent and recognized, the results of a compre- 
hensive engine study with the author’s equipment should funda- 
mentally advance our recognition of the perverse idiosyncracies 
of combustion in an engine. 

It is admittedly difficult to make a paper of this nature com- 
plete in its coverage of all details of technique. A further word 
from the authors would be appreciated concerning the extent to 
which this equipment is independent of the diagnostic and manip- 
ulative sensitivities of the individuals who have made it work. 
What pitfalls would await the uninitiated adventurer following in 
their path? 

Is the present form of the equipment the first one assembled? 
What further improvements are needed or expected? What is 
the useful life of the mechanical parts? How quickly do the 
photocells deteriorate? How frequently must the photocells be 
recalibrated? Are separate means advisable to check on the in- 
terchangeability of the photocells? How high an engine speed 
would this equipment stand with maintained accuracy? What 
maximum and minimum pressure and temperature limits can be 
handled? Will this equipment show up detonating frequencies 
independent of ‘‘hash?” 

Is it worth while to consider building in a differentiating circuit 
to determine the rates of pressure and temperature variations? 

How much of the lack of reproducibility between cycles might 
be due to the engine-operating controls themselves? How long 
is the engine run in before data are taken? Is ordinary commer- 
cial injection equipment believed to be satisfactory here? An- 


6 Project Engineer, Pratt & Whitney Aircraft, East Hartford, ~ 
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other oscillogram of fuel injections rate versus time might be 
imagined as an additional complication. 

Then, from an analytical viewpoint, are there believed to be any 
limitations on Equation [5], as applied to an engine combustion 
chamber? What is the magnitude of the monochromatic trans- 
missivity of the flame? What would be the corresponding re- 
flectivity and absorptivity of the combustion-chamber surfaces in 
operation? What kind of temperature does one measure along a 
radiant beam passing through a gas mass in violent physical and 
chemical agitation? Has this temperature been employed to cal- 
culate the radiant component of the heat transfer to the combus- 
tion-chamber walls? 

Further data are most certainly awaited with keen interest. 


J. P. Mzeroor.? In studies such as the author has made, the 
apparatus involves a small window inset in a cylinder. Beyond 
doubt the coolest part of the cylinder is the wall. The question 
then arises, can a really good idea of the temperature inside the 
cylinder be obtained from measurement through such a window? 

Why was the measurement of pressure taken with the pre- 
sented method instead of with the piezo-electric crystal method 
in which the indicating parts are inertia-free and therefore more 
correct? ° 


Hans Bouustav. How many pictures could be taken before 
the window had to be removed and cleaned? Were the measure- 
ments influenced by the radiation of the combustion-chamber 
walls? Was it possible to measure the combustion-chamber-wall 
temperature and determine the effect of an increase in such tem- 
peratures on the cleanliness of combustion and flame temperature? 
It is commonly known that after an engine has run for a while 
and is properly warmed up, the combustion becomes cleaner; 
therefore, the temperature of the combustion-chamber wall 
would have an influence on the speed of combustion and tem- 
perature of the flame. 


L. R. Turner.? What is the possibility of using this method 
to measure the temperature of gases which have no radiant par- 
ticles in them? In their experiments, have the authors used addi- 
tives in sufficient quantities such that the radiation is along one 
wave length? It is the writer’s understanding that a method has 
been worked out but he is unfamiliar with the quantities required 
to get black-body conditions at one wave length. 


Leropoip Srravuss.” The authors’ illustrations show the pre- 
combustion chamber as spherical. What would be the effect of a 
differently shaped chamber on the temperature measurements? 
Also, what effect does turbulence within the chamber have on 
temperature measurement? 


AvurHoRs’ CLOSURE 


The discussion submitted by H. B. Nottage includes several 
pertinent questions which one thoroughly conversant with the 
technical theories and principles involved might naturally raise. 
The authors welcome the opportunity to discuss these questions 
at this time but not all of them can be answered conclusively until 
after further experience and experimentation. 

The authors believe that the equipment is quite independent 
of the personal element, once the operators become familiar with 
the basic principles and technique. The equipment is still in its 
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10 Sterling Engine Company, Buffalo, N. Y. Jun. A.S.M.E. 


original form except for a few minor changes. Simplification is 
possible and certainly would be necessary before it could be used 
successfully by nontechnical operators, 

Probably the most desirable improvement in the electro- 
optical temperature indicator would be realized if the equations 
involved were solved electronically so as to produce true tempera- 
tures directly on a single oscillograph screen. This would elimi- 
nate the laborious, time-consuming routine inherent in the present 
arrangement and should improve the accuracy. Needless to say, 
some means of attaining this objective is vigorously being sought. 

The mechanical parts are few and not subject to rapid wear or 
deterioration. The phototubes are used at low light levels, and at 
moderate voltages, consequently they should be usable for a long 
period of time. The original phototubes are still in use and show 
no appreciable deterioration after many hours of service. Any 
change that might occur in the characteristics of the phototubes, 
occasioned by deterioration or replacement, is taken care of by 
strict adherence to the practice of checking the calibration of the 
apparatus each time it is used. 

The authors are of the opinion that the instrumentation will 
operate successfully at any speeds employed in present-day re- 
ciprocating engine practice if proper precautions are observed in 
designing the electrical circuits and the mechanical auxiliaries. 
Since there are no mechanical connections between the engine and 
the electronic indicators, these instruments should follow ac- 
curately, without interference from engine vibrations, any fre- 
quency up to the limit of frequency response of the amplifiers. 

As a matter of possible interest, the electro-optical temperature 
indicator was once sighted toward the sun and the temperature 
thus obtained checked the generally accepted value very closely. 

If, as is hoped, the true flame temperature can be indicated 
directly on a single oscillograph screen as a function of crank 
angle, temperature rate curves will definitely be of interest. 
With the present equipment, however, apparent temperature rate 
curves do not seem to be in order as both flame emissivity and 
temperature variation affect the rate. 

Lack of cycle reproducibility is believed to be inherent in the 
combustion process, since the injector rack is usually locked in a 
fixed position while data are being obtained. This conclusion 
seems reasonable because of the manner in which ignition origi- 
nates and combustion is propagated in a Diesel engine, and more 
particularly so in the high turbulence, precombustion-chamber 
type used in the present investigation. 

The engine is run at the desired operating conditions for a 
period of approximately one hour before data are taken. Special 
injection equipment may prove useful, but the final answer to 
this question must await further experimental work. The use of 
injection rate diagrams has been considered. While such addi- 
tional data would undoubtedly contribute useful information, the 
necessary equipment has not been developed as far as the authors 
are aware and consequently this improvement must await future 
development work. 

The authors consider that Equation [5] is just as applicable 
to luminous flames in engine combustion chambers as it is to open 
flames. The equation must fail, however, when a luminous flame 
becomes nonluminous and the radiation becomes discontinuous 
with respect to wave length. In the present experiments the 
luminosity of the Diesel flame is of such magnitude that the equa- 
tion is considered applicable. Further experimental work will be 
necessary to establish its validity in the case of the less luminous 
Otto-cycle flames. 

The monochromati¢ transmissivity varies between wide limits. 
As yet, the authors have made no attempt to separate the trans- 
missivity of the flame and that of the sooty quartz window. 
KL values for the flame-window combination range from 0.001 to 
1.0. The corresponding absorptivity of the combustion-chamber 
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walls is unknown but must be quite high due to the surface coat- 
ing of soot. It should be noted that the temperature obtained 
with the electro-optical pyrometer is exponentially weighted 
toward the highest temperature existing in the radiant particles 
visible to the instrument. This fact, together with the high 
turbulence and the suspected lack of thermal equilibrium, may 
be largely responsible for the lack of reproducibility. 

Separation of the flame absorptivity and the window absorp- 
tivity would be necessary before the radiant heat transfer could 
be calculated. A new type of quartz window is now being em- 
ployed which may make such separation feasible. 

With respect to Mr. Megroot’s question concerning the effect 
of window temperature, it is unlikely that radiation from the 
window itself can appreciably alter the electro-optical pyrometer 
readings because the window temperature is far less than the 
flame temperatures, and the exponential averaging previously 
mentioned tends to nullify minor radiation effects. 

Experience with the piezo-electric pressure indicator in previous 
combustion studies at the University of Wisconsin indicated that 
it does not hold to its calibration well. The electro-optical 
pressure indicator has proved entirely satisfactory in that re- 
spect. Furthermore, the natural frequency of the diaphragm 
used in the electro-optical pressure indicator is approximately 
200,000 cycles. Since the ability of the diaphragm to respond to 
rapid changes in pressure is determined by its natural frequency, 
the response is satisfactory up to the limit of the frequency 
response of the amplifiers used. 

Replying to Hans Bohuslav’s questions, the frequency with 
which the quartz window must be cleaned is a function of the 
engine-operating conditions, with a minimum of about 15 minutes. 
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Since the paper was presented, a self-cleaning window has been 
developed which apparently remains clean indefinitely. No at- 
tempt has been made in the present investigation to measure 
combustion-wall temperatures. In view of the exponential 
averaging, previously mentioned, radiation from the relatively 
cool wall is negligible as far as flame-temperature measurement 
is concerned. 

Concerning the discussion submitted by Mr. Turner, since the 
method in question measures the temperature of the soot particles 
in the flame it is not applicable to nonluminous flames. Inasmuch 
as the Diesel flame is luminous, there has been no occasion to 
employ additives in the present investigation. The possibility of 
additives affecting the combustion characteristics would seem 
to make their use undesirable if it can be avoided. The authors 
have not experimented with the amount of additives required to 
approach black-body conditions, nor do they know of any refer- 
ence to experiments of this nature. Reference (10) describes a 
method whereby the spectral radiation of an additive was used to 
measure flame temperatures, but black-body conditions were not 
approximated. 

The question by Leopold Strauss concerning the effect of 
combustion-chamber shape must await further experimentation 
because, as yet, the apparatus has been employed on only one 
type of combustion chamber. Turbulence should have no effect 
upon the temperature measurements if thermal equilibrium exists. 
Lacking thermal equilibrium, however, temperature variations 
would be indicated if “hot” and ‘cold’? masses were swept 
alternately in front of the combustion-chamber window. 

The authors wish to express their appreciation to those who 
contributed to the discussion of the paper. 


Polar Diagram for Tuning of Exhaust Pipes 


By TROELS WARMING,! MILWAUKEE, WIS. 


In this paper is shown how a simple polar diagram, not 
requiring any higher mathematics, can be used for calcu- 
lation of the natural frequency of the exhaust system on a 
Diesel engine. In general this diagram can be used for 
determination of all the natural frequencies of any sys- 
tem consisting of gases or fluids in containers of complex 
form. 


URRENTLY there seems to be an increasing interest in 

( , tuning of exhaust pipes on Diesel engines, see, for in- 

stance, a paper by P. H. Schweitzer;? and it may there- 

fore be of advantage to have an easy method to determine the 

natural frequencies of such systems. A polar diagram, similar to 

the diagram developed some years ago by F. M. Lewis? for cal- 

culation of crankshaft systems, may be used for this purpose. 

It will be found, however, that gas systems are much simpler to 
deal with than shaft systems. 
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AMPLITUDE 


Fie. 1 OscinLaTIoN IN Prez OPEN at ONE END 

Consider first the system shown in Fig. 1, consisting of a pipe 
with constant area and open at oneend. The natural frequency 
of this system is 


v l5v 
N = — - 60 = — vibrations per min 
~4h lo ? 
where ly is the length of the pipe in inches, and v the sound velocity 
in inches per second. The exhaust from a Diesel engine may 
with good approximation be considered as atmospheric air for 
which the sound velocity is 


. v = 585 V1 + 460 in. per sec 


where t is the temperature in degrees F. 


1 Mechanical Engineer, Nordberg Manufacturing Company. Jun. 
A.S.M.E 

2‘Tmproving Engine Performance by Exhaust-Pipe Tuning,” 
by P. H. Schweitzer, Journal of the American Society of Naval Engi- 
neers, vol. 56, May, 1944, p. 185. 

2“‘Torsional Vibration in the Diesel Engine,’’ by F. M. Lewis, 
Trans. Society of Naval Architects and Marine Engineers, vol. 33, 
1925 p 109. 

Contributed by the Oil and Gas Power Division and presented 
at a meeting of the Cleveland Section, Cleveland, Ohio, May, 1945, 
of Tur AMERICAN SOCIETY OF MECHANICAL ENGINEERS. 

Note: Statements and opinions advanced in papers are to be 
understood as individual expressions of their authors and not those of 
the Society. 
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Fig. 2 Poutar Diagram ror System IN Fria. 1 


When the air in this pipe oscillates, amplitude and pressure 
at the various points are given by a sine and a cosine curve, Fig. 1. 
Amplitude and pressure at a point B in the pipe are therefore 
the co-ordinates for the corresponding point B on the circle 
shown in Fig. 2, where 


arc AB = ; - 90 deg 


The scales used in this polar diagram are of no importance for 
the determination of natural frequencies since only the relative 
values are needed. Amplitude and pressure are therefore 
both measured in inches on the diagram. It is of importance, 
however, that these scales are not influenced by the area of the 
pipe. 

Consider now length PQ = I, of the pipe. 
circle arc PQ = arc, is 


The corresponding 


l 
arcy = A =O) GOP a tecuerstatratty hus eee [3] 
Substituting 4 from Equation [1] 
6N1 
arc, = Ei dan oUt, Tete co te ae {4] 


The air particles between P and Q oscillate independently of 
what takes place in the rest of the pipe, with the exception that at 
P and Q certain pulsating pressures are supplied and certain 
volumes of air are pumped in and out. If these pressures and 
volumes were provided by oscillation at the same frequency but 
in pipes of different diameter, the oscillations between P and Q 
would not be affected. 

. This makes it possible to treat more complicated systems, 
Figs. 3and 4. Each length of pipe with constant area will have a 
corresponding arc of a circle determined by Equation [4]. Where 
two pipes of different area are connected, the pressure does not 
change, but the two amplitudes, a, and an4,, at the point of 
connection are inversely proportional to the areas, S, and Sp+4; 


Se 


It may be added that the condition is different from that of 
constant flow, in which case the pressure changes where the area 
changes. This pressure change is necessary to accelerate the air 
passing to a different velocity. In a vibratory system, as con- 


s 
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ha pacoonacileea oak 


Fie. 3 Morr Compricatep System 
(Letters S indicate pipe areas.) 


AMPLITUDE 


‘ PRESSURE 


Fic. 4 Pouar DiaGRAM For SYSTEM IN Fic. 3 


sidered here, the air in the whole pipe is accelerated simultane- 
ously, thus not causing any sudden pressure changes at any point 
in the system. 


CONSTRUCTION OF POLAR DIAGRAM 


In order to explain the construction of the polar diagram in 
Fig. 4, it is first assumed that the natural frequency N of the 
system is known. The first point A; of the diagram is chosen 
arbitrarily on the abscissa axis because it gives the condition at 
the closed end A of the system where the amplitude is zero. 
Point B,, which gives the condition at the end B of the first pipe, 
is located on a circle arc through A; determined by Equation 
[4] and with center in point O. Point B, gives the condition 
at the beginning of the second pipe, where the pressure is the 
same as at B,, while the amplitude is determined by Equation [5]. 
Point B, is therefore found by simply multiplying the ordinate for 
B,, as measured on the diagram, by this proportion. The circle 
arc B,C, has center in point O and is determined by Equation [4]. 

The construction is continued through all the steps of the 
system, but, with atmospheric pressure at the open end, the 
last point must fall on the ordinate axis. However, if the dia- 
gram had been drawn for a frequency different from the natural 
frequency of the system, the last point would not have fallen on 
this axis. We have therefore in the polar diagram a method of 
checking whether a given frequency is the natural frequency. 

It is added that the diagram may be drawn just as well by first 
choosing point H, arbitrarily on the ordinate axis and then pro- 
ceeding the other way through the system. If point A, then 
falls on the abscissa axis, the frequency used was the natural 
frequency. 

So far only the fundamental vibration of the system has been 
considered, but for higher frequencies the polar diagram, in- 
stead of 90 deg, will go through 270 deg, 450 deg, etc. Fora sys- 
tem closed at both ends, the diagram will go through 180 deg 
(or 360 deg, 540 deg, etc.) thus beginning and ending on the 
abscissa axis. If the system is open at both ends, the diagram 
will begin and end on the ordinate axis. 

Before calculating the natural frequency of an actual exhaust 
system, it must be simplified to a system of pipes with constant 
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area, see, for instance, Figs. 3, 5, or 7. Parts of complicated shape 
are changed so that they keep the same volume and the same 
length, the length being measured as the distance the gases flow. 
Exhaust pots and mufflers are thus considered as short pipes with 
large area. At the open end of a pipe is added a length equal to 
0.4 times the pipe diameter in order to provide for the fact that 
the air particles,just outside the pipe take some part in the vibra- 
tion. : 

A frequently used exhaust system has a pipe, with constant 
area S and length lJ, and some additional volume V as, for in- 
stance, valves or exhaust pots, near the closed end. If V/S <1/2, 
the system may be calculated as one straight pipe, Fig. 1, of 
length h = V/S +1. The error involved is less than 3 per cent. 


CALCULATING NATURAL FREQUENCY 


As an example, the natural frequency will be calculated of the 
system shown in Fig. 5, assumed filled with exhaust gas at an 
average temperature of 500 F. The sound velocity, Equation 
[2], then is 


» = 585 V/500 + 460 = 18,100 in. per sec 


A rough estimate of the natural frequency is obtamed from 
Equation [1] by considering the system as one straight pipe 
159 15-18,100 


= 1700 vibrations per min 


ty 2 100" 


t= 160" 


Fie. 5 Exuaust System CaLcuLaTEep IN EXAMPLE 


~ The first polar diagram in Fig. 6 is drawn for this frequency. 
The angle of the circle arcs, all of which have center in point O, 
is from Equation [4] 


6Nly,  6°1700- 
Shh Sea 0.564:1, deg 


aha ir 18,100 


where l, is the corresponding pipe length. 
Point A; is chosen somewhere on the abscissa axis. 
found by 


Point B, is 


arc A,B, = 0.564:l = 0.564 + 20 = 11.28 deg 
Point B, isfound by multiplying the ordinate for B,; by 


Si _ 200, 
Sh 7 


Point C2 is found by 
arc B,C, = 0.564 +1, = 0.564 + 40 = 22.56 deg 
Point C3 is found by multiplying the ordinate for Cz by 


S: _ 100 
= —— = 2.5 
S: 40 


Point D; is found by 
arc C,D; = 0.564-1, = 0.564- 100 = 56.4 deg 
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The diagram shows that the estimate of 1700 vibrations per 
min was too high since the total angle described is 124.1 deg 
instead of 90 deg. A second estimate is therefore made 


90 
N = Pree 1700 = 1280 vibrations per min 
For this frequency a new polar diagram is constructed with the 
angle of the circle arcs 


6+ 1230 - I, 


= 0.407 -1 
18,100 0.407 - l, deg 


arc, = 
The procedure is the same as before and the result this time is a 
total described angle of 98.2 deg. 
By extrapolation from the first two estimates is found 


1700 — 1230 


SCAT alAsinl 210810 


(98.2 — 90) 


= 1080 vibrations per min 


A polar diagram for this frequency has a total described angle 
of 89.4 deg which is quite close to 90 deg. However, by interpo- 
lation between the last two tries will be found 


1230 — 1080 


=i pee sari SS 
ny, 080 + 98.2— 80.4 


(90,— 89.4) = 1090 vibrations per min 


which is the natural frequency of the system with good accuracy. 
During the construction of the polar diagram, difficulties may 
arise with points falling either outside the paper or too close to 


the center to give an accurate construction. The following rule 
will then be helpful: Any point in the diagram may be moved 
on a radius, closer to or away from the center, and the construc- 
tion may then be continued from the new location. 


BRANCHED SYSTEMS 


Branched systems need a separate polar diagram for each 
branch and for the main pipe, all diagrams using the same fre- 
quency, Figs. 7 and 8. The diagrams for the branches are begun 
at the end of the branch and worked toward the junction. At 
this point the pressure p must be the same in both branches so the 
diagram for one of them, in this case branch II, is reduced ac- 
cordingly. The diagram for the main pipe now can be drawn 
beginning at the junction point. The pressure p must be the 
same as in the branches, and the amplitude a7 is determined by 


G7 * Sz = a4 * Sy + ag? Se 


where a4, a, and a7 are measured in inches on the diagram, and 
S4, Se, and S; are pipe areas. The diagram is then completed 
and, if the last point falls on the ordinate axis, we have the natural 
frequency of the branched system. 
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Discussion 


T. M. Rosiz.‘ What is the effect of discharging into a concrete 
pit and taking the exhaust pipe out somewhere before the end of 
the pit? In other words, supposing the exhaust pipe enters the 
pit horizontally at one end and leaves the pit vertically at some 
point before the other end of the pit. Would the vibrations in 
the pit, due to the exhaust gases hitting the end wall and bouncing 
back, affect the calculations and if so, what correction should be 
made? 


R. O. Bernauer.§ The author has been a co-worker with the 
writer at Nordberg for the past 5 or 6 years and in that time has 
worked on the critical speeds of crankshafts, parallel operation of 
synchronous generators, and various other vibration problems. 
He has always had a way of presenting anything very simply and 
has used polar diagrams a great deal in solving problems. This 
paper is his latest application of the polar diagram. 

The writer is not going into a technical discussion regarding the 
use of it. In 1925, Professor Lewis used the polar diagram in the 
determination of critical speeds of crankshafts and the author is 
using it in an analogous manner in this paper. In parallel-opera- 
tions problems he uses vector diagrams which are closely related. 
The writer congratulates the author on his unique and simple 
paper and thinks it should be of great value to engineers because 
of its simplicity. The paper deals with relative lengths from 
which pressures and amplitudes in any system may be secured. 

In his presentation the author did not go into branch or net- 
work systems at any length but did give examples in his paper. 
These are obtained by the confluence of pipes and are very diffi- 
cult to handle by other methods. His method is important in 
that it_ presents a simple way of solving branch or networks 
problems. 


L. R. Turner.® In connection with the polar diagrams, what 
happens where there is a restriction at the end of the pipe? 


L. M. Ticuvinsxy.’ The determination of natural frequencies 
even of simple members of a Diesel engine is required in order to 
obviate a possible condition of resonance. Their calculations are 
rather lengthy and require accurate procedures. 

The simple and elegant method presented by the author for the 
same purpose invites a question as to its accuracy. It is believed 
that the method is practical for designers; however, it would be 
interesting to know the error when comparing the calculated and 
the graphically obtained values of natural frequencies of a system 
such as shown in Fig. 7 of the paper. 


V. L. Maurry.? The writer checked the method described in 
this paper for an exhaust system connecting several Diesel en- 
gines to a common muffler with a long tail pipe. The system con- 
tains chambers of various volumes and pipes of various diameters 


4 General Diesel Sales Division, Fairbanks Morse & Company, 
Chicago, Ill. Mem. A.S.M.E. 

’ Nordberg Manufacturing Company, Milwaukee, Wis. 

‘ National Advisory Committee for Aeronautics, Cleveland Air- 
port, Cleveland, Ohio. 

7 Vice-President in Charge of Engineering, American Bearing Cor- 
poration, St. Louis, Mo., and Indianapolis, Ind. Mem. A.S8.M.E. 
Rae Division, Harvey Machine Company, Inc., Los Angeles, 
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and lengths and it was a tedious job to find by usual calculations 
the natural frequency of the gases in this system. The polar- 
diagram method gave a great saving of time without any sacrifice 
of accuracy. In addition, this method shows very clearly what 
changes should be made in a system if a certain natural fre- 
quency is to be obtained. 

The writer believes that, where due to the difficulty of calcu- 
lating the natural frequency, the latter was often only estimated, 
with the polar-diagram method available the natural frequency 
actually will be computed and this will lead to better design. 
Furthermore, the writer believes that this method will be very 
helpful not only for tuning exhaust pipes of internal-combustion 
engines, but also for designing exhaust mufflers which up to date 
are not too satisfactory, and possibly in designing injection sys- 
tems of compression-ignition engines. 

The only criticism which pertains, however, not to the content 
of the paper, is in respect to the use of the designation ‘‘are” in- 
stead of the word ‘‘angle;”’ 90 deg multiplied by a ratio gives a 
certain number of degrees and not an are. In laying out the 
diagrams, again one uses angles on the protractor and not arcs. 
The polar diagram itself deals with angles, not ares. 

The paper is written very clearly and will be appreciated by 
all who have to deal with natural frequencies of gases or fluids in 
various systems. It constitutes a valuable contribution to the 
engineering literature. ; 


AutHor’s CLOSURE 


Mr. Robie asks about a concrete pit. Usually it will be sufficient 
to substitute a straight pipe, of the same length as the average 
path of the gases through the pit, and with an area giving same 
volume as the pit. If the part of the pit beyond the vertical out- 
let is very long, it may be necéssary to consider it as a branch to 
the system. This also should be done if higher orders of vibra- 
tion are considered. 

Mr. Bernauer’s kind remarks are very much appreciated. The 
association with him, through the past years, has always been 
very pleasant. 

In reply to Mr. Turner, it can be said that if the open end of a 
pipe is restricted by a plate with a hole, this hole is considered 
as the open end. We then must add a pipe, of same diameter d 
as the hole, and of length 0.4 d, on both sides of the plate. It will 
be found this way, that the smaller the hole, the more effective 
is the mass in the added pipe, and consequently the lower is the 
natural frequency. 

Mr. Tichvinsky questions whether the simplicity of the polar 
diagram is at the expense of accuracy. This is not the case. 
With a careful graphical construction the error should be less 
than 1 per cent and if that is not satisfactory, it is quite simple to 
calculate the steps of the diagram, if necessary with logarithms. 
Usually that does not serve any purpose, however, because the 
engine does not behave that accurately. 

Compared to other published methods, the polar diagram is 
more accurate because it takes account of the inertia and the 
compressibility of the air throughout the system. In other 
methods, the inertia in exhaust pots or other large volumes is dis- 
regarded. The error involved is small, but the polar diagram 
takes care of the inertia automatically. 

It is interesting that Professor Maleev has found the polar dia- 
gram useful and time-saving. That is what the author hoped to 
accomplish with this paper. 


Substitution of Lower-Quality Industrial 
Diamonds in Diamond Dresser Tools 


By H. WHITTAKER,! CRANFORD, N. J. 


A shortage of grinding-wheel dresser diamonds, which 
threatened the war effort early in 1942, led toaninvestiga- 
tion of lower-quality diamonds hitherto not used for 
dresser tools. According to estimates made available by 
the British Ministry of Supply to the War Production 
Board, the supply of those qualities and sizes of industrial 
diamonds preferred and used in large quantities for truing 
and dressing abrasive grinding wheels was expected to be 
critical. The War Metallurgy Committee on the recom- 
mendation of the Office of Production Research and De- 
velopment of the War Production Board undertook a 
study, as “‘restricted’’ Project NRC-507, to determine if 
lower-quality, more plentiful stones could be substituted 
for the preferred qualities and sizes for some or all wheel- 
truing and dressing operations. The investigation was 
carried out in the Research and Development Laboratories 
of the Crane Company, Chicago. This paper, based on 
that work, has been released for publication by the OPRD. 
In the first part of this investigation, a survey of users of 
diamond dresser tools established definitely that the 
higher qualities were preferred and used almost exclusively 
fer dressing grinding wheels. The results of this survey 
and a review of the literature on the practices and the 
factors which influence the use of dresser tools were re- 
ported in a previous paper.? In order to provide positive 
information on the possibilities of using lower-quality 
stones, it was necessary to classify the seemingly infinite 
variety of industrial diamonds according to distinguishing 
characteristics and to test each type under controlled 
conditions which simulated grinding practices. 


EXAMINATION AND CLASSIFICATION OF INDUSTRIAL DIAMONDS 


RACTICALLY the entire production of industrial diamonds 
P: distributed in bulk by the Diamond Trading Company, 

London, in the form of six general classified assortments. 
These assortments, from highest to lowest in value, are known 
commercially as “industrial serie,” ‘finest industrials” ‘‘in- 
dustrial smalls and mine rejections,” “industrial cleavages and 
chips,” “better bort”? (not distributed in quantity at present), 
and “crushing bort.”” 

The ‘industrial serie” contains the better, whole, mostly trans- 
parent stones from all producing areas. The stones are classified 
according to size, color, shape, and quality, and range in size 
from about 0.25 carat to 7 carat. 

The assortment known as “finest industrials” is made up from 
the production of Consolidated African Selection Trust and Gold 
Coast fields. It contains whole stones, twinned crystals, called 


1 Consulting Engineer, Mem. A.S.M.E. 

“Substitution of Lower Quality Diamonds,’’ War Metallurgy 
Committee Progress Report W-57, Nov. 6, 19438. 

Contributed by Special Research Committees on Cutting Fluids 
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macles, flats, and a fairly large number of chips, or fragmented 
stones all in small sizes averaging about 0.05 carat, with some 
stones as large as 0.25 to 0.5 carat. 

“Industrial smalls and mine rejections” is roughly 70 per cent 
C.A.S.T. production with the balance from Forminiere and 
Angola. The stones are all under 0.1 carat in size down to fine 
“sand” 60 per carat. Over one half the assortment is made up 
of chips and fragments dark in color, and the remainder is whole 
stones transparent to dark colored. 

“Industrial cleavages” range from transparent to colored frag- 
ments, as opposed to single, whole stones, in all sizes. Darkest 
colors and frosted and coated fragments are usually not included. 

All of the remaining materials from all producing areas not in- 
cluded in the foregoing assortments are lumped together in 
“crushing bort,’’ which is intended for crushing to powder. The 
dark-colored coated fragments and stones characteristic of the 
Belgian Congo predominate in crushing bort; and for this reason 
the term crushing bort is often used synonymously with Belgian 
Congo. “Better bort’’ consists of the single, whole, and sound 
stones selected from crushing bort. 

The main source of dresser-tool stones is the industria! serie. 
The most plentiful sources, and practically the only sources, of 
substitutes for dresser-tool stones over 0.25 carat in size to sup- 
plement the supply are crushing-bort and better-bort assort- 
ments. Samples representing the preferred types of dresser 
stones and the possible substitutes were obtained for testing. 
Stones were selected from 15 of the 32 categories in the industrial 
serie. These categories were second yellow, first and second 
colored, first gray, first brown, dark and darkest brown, first and 
second gray brown, second river, and common gray colored 
rounds, first yellow, first brown, first gray, and second gray 
shapes. A few samples of ‘‘cleavages’’ were obtained, representa- 
tive of the classifications known as colored, first and common 
gray cleavages. A 5000-carat parcel of bort completed the 
sampling of materials. It was made up of selections from 
crushing bort and better bort and it covered most of the varieties 
found in these assortments, such as common coated stones, 
common coated cleavages, black-bort stones, and colored-bort 
stones. 

Representative samples of each of the different categories 
were examined with a binocular microscope, a petrographic 
microscope, and other methods for soundness and distinguishing 
characteristics. The higher-quality stones in the industrial serie 
were transparent whole stones of definite crystal habit. Most 
of the stones contained small inclusions and several stones were 
fractured internally. ‘The lower-quality stones in the serie were 
darker in color, of low transparency, and they contained inclusions 
in fairly large amounts. The crystal habit of the lower qualities 
was not so well defined, and these stones had more striated and 
irregular surfaces than had the other qualities. Three quarters 
of the transparent stones showed fluorescence, which varied in 
degree from vivid to pale shades of blue and yellow, when the 
stones were exposed to the light of a mercury-vapor lamp. About 
one half of the stones exhibited phosphorescence. A majority 
of the diamonds examined under polarized light showed the 
presence of internal strain. Density measurements made by 
the immersion method using carbon tetrachloride on about 50 
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of the larger stones gave values of density ranging from 3.510 to 
3.526 on the basis of water at 4 C. 

In contrast with the samples from the serie, a majority of the 
stones examined from the bort assortment were translucent to 
opaque; the degree of opacity appearing to vary with structure 
and crystal habit. Most of the stones were well defined in 
erystal habit. A considerable number of the stones consisted of 
cleavages, fragments, and aggregated crystal masses. Those 
stones that were sufficiently translucent to enable internal in- 
spection were observed to contain large areas of inclusions. 
Only about one fifth of the opaque coated stones examined 
showed fluorescence of dull intensity, and a few showed phos- 
phorescence. A few of the more opaque stones were radiographed 
to determine whether internal defects could be detected by this 
method. The stones were immersed in Russian mineral oil in 
cellophane containers and exposed for 5 sec at 5 in. focal distance 
to unfiltered X rays produced by a molybdenum target operated 
at 20 kv and 18 ma. The resulting shadowgraphs, several of 
which are given in Fig. 1, indicated that it may be possible to 
detect flaws and inclusions in opaque stones by this method. 
However, the detection of small inclusions of foreign elements 
having low absorption coefficients comparable to carbon is not 
possible; and thus limitations may be placed on the usefulness 
of radiography. 

It was desirable to supplement the present commercial classifi- 
cation, which is rather loose, with a classification which would 
make possible closer distinctions among the various types of 
diamonds for testing purposes. While the physical properties 
of diamond establish its suitability for wheel dressers, varia- 
tions in structural characteristics and soundness are the criteria 
for selection of individual stones. 

The structural characteristics and factors to be considered are 
as follows: 


Wholeness; whole stones or fragments. 
Crystal habit; largely determines shape. 
Opacity; degree of impurities. 
Homogeneity; continuity of structure. 
Soundness; number of fractures and cracks. 
Size. 
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Compared with the present loose classification, crystal habit 
is analogous to shape, quality is a combination of wholeness, 
opacity, homogeneity, and soundness, and color is partially analo- 
gous to opacity. Practically all of the industrial serie stones 
were whole single crystals of definite crystal habit. Of the 5000 
carats of better and crushing bort available for this investigation, 
about 3000 carats were whole stones. Well over 90 per cent of 
the bort exhibited definite crystal habit. One half of the whole 
stones were single crystals, 17 per cent were twinned crystals, 
30 per cent were aggregated crystal masses, and 3 per cent were 
of no discernible crystal habit. 

It was observed that the stones occurred in one of three crystal 
forms: Octahedron, dodecahedron, and cube, and combinations 
of these forms. No other crystal forms were evident except 
modification to tetrahexahedron. Thirteen classes, grouped 
under five general headings, illustrated in Figs. 2 and 3, were 
defined which served to include all of the observed forms and 
combinations, as follows: 

Dodecahedron. Class No. 1, ‘distorted dodecahedron:” Ellipti- 
cal, dodecahedron crystals lengthened in the direction of trigonal 
axis of symmetry with faces more or less modified to tetrahexa- 
hedron. 

Class No. 2, ‘‘normal dodecahedron:’’ Symmetrical dodeca- 
hedron crystals with or without curved faces resulting from 
tetrahexahedral influence. 

Octahedron. Class No. 8, “octahedron - dodecahedron:’’ 
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Twenty-faced combination octa-dodec crystals with octahedron 
faces usually curved and shieldlike. 

Classes No. 4 and No. 5, ‘‘octahedron:” Octahedron crystals 
with or without rounded faces caused by transparent or glassy 
growth plates, and giving rise to striated edges. No. 5 is same 
except growth plates are translucent. 

Combination Octahedron, Dodecahedron, and Cube. Classes 
Nos. 6 to 9 constitute a series between an octahedron and a cube 
with the sizes of the octahedron faces decreasing as the cube faces 
increase, the dodecahedron faces varying accordingly: No. 6 
octahedron faces large, others small; No. 7 octahedron faces 
large, other faces larger, but not as large as the octahedron faces; 
No. 8 equal-sized octa-dodec-cube faces, which configuration 
produces a spheroidal stone; No. 9 minor octahedron and do- 
decahedron faces relative to cube faces. 

Cube. Class No. 10 and No. 11, “‘cube:” Six-sided crystals 
with rough terraced growths on faces. No. 11 has convex, con- 
cave, pyramidal, or flat faces. 

No Definite Crystal Habit. Classes No. 12 and No. 18, no habit: 
Crystal habit not evident with glassy or transparent rind. No. 13 
has translucent to opaque coating or rind. 

Opacity, degree of diaphaneity or ability to transmit light, is a 
measure of the amount of impurities in the form of dispersed 
included matter. A scale was set up based on the amount of 
light transmitted through a stone using a 10-w Mazda lamp con- 
tained in a white-walled box covered with a 1/s-in-thick sheet of 
opal glass. An opaque stone that did not transmit any light 
under these conditions was assigned a value of 7 and a water-— 
white stone was assigned a value of 1. Five other stones having 
equal intermediate steps of increasing opacity were selected to 
serve as standards of comparison to give a complete graded 
range of opacities. 

Homogeneity, or continuity of structure was observed to be 
related somewhat to crystal habit and opacity. A large number 
of stones which fell in classes 6 to 9, ‘combination octahedron- 
dodecahedron-cube,” and which had high opacities 4 to 7, con- 
sisted of a transparent octahedron crystal nucleus surrounded 
by a secondary growth containing large amounts of inclusions, 
thus giving rise to a discontinuity in structure at the interface 
between core and rind. The thickness of the secondary growth 
appeared to be rather definitely related to the width of the dodeca- 
hedron face in the ratio of 0.4-0.7 to 1, depending on the extent 
of rounding of the face. Transparent crystals, classes 1 to 5 in 
particular, appeared to be built up of layers on a basic octa- 
hedron form. The interfacial boundaries in this case do not 
result in as marked a discontinuity in structure because. the 
primary and secondary growths are about equal in composition. 

During the detailed examination of over 5000 carats of bort 
and higher-quality material, a great deal of evidence was found 
in support of the theory of diamond-crystal formation stated by 
Polinard, “‘one might explain that all forms of diamonds are due 
to the development and to the piling up of octahedron plates 
which superimpose themselves on the faces of an octahedral 
nucleus.” 

Soundness is used to indicate the number of weaknesses which 
may be contained in a stone as a result of cracks, fractures, and 
fissures. The soundness of transparent stones can be judged by 
inspection with a microscope using suitable lighting. Radiog- 
raphy may be of some assistance in detecting flaws in opaque 
stones which cannot be inspected by ordinary means. The 
Norton Company has used an impact test to grade bort for suita- 
bility for dresser tools with fair success. That firm has tested 
samples of most of the classes herein described. According to 
the results of these tests, the impact resistance of single crystals 


\ 


averaged higher than either twinned or aggregated stones in like - 


classes. Discernible differences in impact resistance were ob- 
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Class No. 11 
Opacity—7 
Opaque Cube 
Diamond inclusion 


shown 
Class No. 6 
Opacity—4 


Opaque Combination 


Class No. 4 
Opacity—5 
Opaque Octahedron 
Inclusions shown 


Class No. 11 
Opacity—3 
Translucent Cube 
No inclusions 


Fic. 1 RapiograpHs or Opaque DIAMONDS 
(‘The black outline which surrounds the radiographs is caused by aluminum foil which was used to support the stones.) 
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Class No. 1 Class No. 2 Class No. 3 
Opacity—1 Opacity—1 Opacity—1 
Distorted dodecahedron Normal dodecahedron Octahedron-dodecahedron 


Class No. 4 Class No. 5 Class No. 6 
Opacity—2 Opacity—3 Opacity—4 
Octahedron Octahedron Combination habit 


Transparent growth 


Class No. 7 Class No. 8 Class No. 9 
Opacity—4 Opacity—3 Opacity—4 
Combination habit Combination habit Combination habit 


Fig. 2 ILLusrrations or CLassés OF DIAMONDS 
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Class No. 10 Class No. 11 Class No. 11 
Opacity—3 Opacity—2 Opacity—1 
Cube Cube Cube 


Class No. 11 Class No. 6 Class No. 6 
Opacity—7 Opacity—4 Opacity—5 
Cube Combination habit Combination habit 


Coating broken away 


Fie. 3 


tained for six classes tested. The more rounded stones ex- 
hibited greater resistance to impact than did stones having sharp 
corners. The classes fell in the following order: Classes Nos. 9, 
8, 7, 6, 10, and 11, from 92 per cent unbroken grain to a low of 
50 per cent unbroken grain. It may be possible to adjust the 
impacting force to just destroy flawed and weak stones without 
damaging sound stones and thus produce a sharper and more 
efficient separation. 

The classification herein proposed and used experimentally on 
over 5000 carats of industrial diamonds accounts for the most im- 
portant factors which are thought to influence the utility of in- 
dustrial diamonds for dressing purposes. Whereas, further 
refinements may be needed, it is flexible and allows for all possible 
types and combinations of crystal form, structural differences, 
and degrees of impurities. The author has examined, in addi- 
tion to the materials herein described, original parcels of all 
assortments among which he has not found any stones that could 
not be fitted into this classification, in particular with respect to 
crystal habit. : 

A stone may be completely defined for the purpose by designat- 
ing (1) the class number, from 1 to 13, which defines its crystal 
habit and shape; (2) the opacity number, from 1 to 7, which 
indicates degree of impurities; (3) the extent of flaws; and in 
addition, (4) by indicating whether a whole stone, a fragment, a 
twin or an aggregated crystal mass. Industrial-serie materials 


Coating broken away 


ILLUSTRATIONS OF CLASSES OF DIAMONDS AND CoaTED STONES 


fall in classes 1 to 4 in the lower opacities, numbers from 1 to 8, 
with some exceptions; for example, dark-brown rounds having an 
opacity of 5. The major portion of material in crushing bort 
falls in classes 6 to 9 and higher opacities, numbers from 3 to 7. 
Minor amounts fall in classes 3 to 5 and 10 to 138. 


GRINDING- WHEEL DressEeR-Toot Test MretHop 


To compare the different varieties of diamonds as to their rela- 
tive suitability for dressing grinding wheels, a test method was 
developed for the rapid quantitative determination of the rate of 
wear under simulated shop-dressing conditions. The test con- 
sisted of the continuous dressing of a l-in. face-grinding wheel, 
the periodic measurement of the amounts worn off the diamond 
and the abrasive wheel, and careful observations of changes ii 
shape and chipping of the diamond point. The essential parts 
of the test apparatus, shown in Fig. 4, were a heavy tool holder, 
specially designed to allow for any desired setting of drag, rota- 
tion, and traverse angles, mounted on the table of a Brown and 
Sharpe No. 2 cylindrical grinder with automatic feed, counter, 
speed control, and fluid system for wet dressing. The diamonds, 
which ranged from 1/2 to 1 carat, were set by the powder hot- 
press method in a large-head steel shank. The shank was 
7/, in. long, 1/:in. diam at the head, and '/,-in-diam shaft. 

To save time and also to enable an inspection of the stone as 
wear took place, the loss on the diamond was measured by means 
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Fie. 4 Werar-Tresr APPARATUS 


of a graphical method. The diamond point in its shank was 
projected on a screen at X100 magnification and the profile 
recorded. After testing, a second profile was made in the same 
position over the first projection. The difference between the 
two profiles was a measure of the loss in height because of wear. 
The area of the flat generated on the diamond point was traced 
at X100 magnification using a reflection microscope and a 
reflex camera and was measured with a planimeter. Height loss 
and area measurements are all that are needed to compute the 
volume loss, provided the shape of the removed portion is known. 
The volume loss was then converted to weight loss using 3.52 
for the specific gravity of the diamond. The computed weight 
loss of the diamond in carats was divided by the volume in cubic 
inches removed from the grinding wheel to give the “unit of wear 
rate.” Since the order of magnitude of the unit of wear rate 
was ten hundred thousandths of a carat per cubic inch, the nota- 
tion ct. X 10~5/cu in. was used, that is, the unit was multiplied 
by 100,000 in order to express the result in terms of whole num- 
bers. 

The test procedure had the advantage of savings in time and 
materials, but precision was sacrificed to some extent. By 
means of this testing procedure a diamond could be tested under 
almost any combination of conditions that might be encountered 
in practice in about a day and a half. To do the same job by 
weighing would consume 5 times as many test wheels and would 
take a week, because of the larger amount of diamond that must 
be worn off to enable an accurate weight determination. The 
greatest source of error was in the assumption of the shape of the 
solid removed in making the volume-loss computation. 

Since the actual shape of the solid figure removed from the 
point of the diamond would be some shape between that of a 
sharp-sided pyramid and a segment of a sphere of equal heights 
and base areas, depending on how close the assumed shape is 
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to the actual shape, the error can be as great as 20 per cent. This 
was checked experimentally by testing a number of stones suffi- 
ciently to enable accurate weight-loss determinations, and in 
most cases the wear rates determined by weight and by com- 
puted volume were well within the possible error of 20 per cent. 
Therefore, for purposes of comparison in this investigation, dif- 
ferences less than 25 per cent, roughly 5 X 107® ct. /cuin. (0.00005) 
were not considered significant; two stones of 15 and 20 X 
10-5 ct./cu in., respectively, would be considered alike in so far 
as could be determined. 

Of interest in the establishment of a set of standard test condi- 
tions for comparison of diamonds and also of practical value were 
the tests that were made to learn the effects of the many factors 
which may influence the rate of wear of diamond dressers in the 
dressing operation. These factors had to be controlled in order 
to place the test results on a comparable basis: 

1 One set of factors is determined by the conditions of opera- 
tion, such as loading of the wheel with metal, wheel speed, rate 
of traverse, feed, and use of a coolant. In practical grinding 
operations the grinding wheel is dressed and trued after the 
wheel has been loaded with metal and the edges are broken down 
from the grinding of work. No appreciable difference in wear 
rates was found when a wheel was dressed after it was loaded 
with metal by grinding a number of pieces of hard tool steel and 
when the wheel was dressed without loading. Changing wheel 
speed from 4000 to 6000 fpm had very little effect on wear rate. 
Nor did varying the rate of traverse between 0.001 and 0.008 in. 
per revolution affect rate of wear to any great extent. Contrary 
to expectations, a 0.001-in. feed caused 1!/; times the wear rate 
of a 0.005-in. feed. However, it will be noted that for each pass 
a 5 times greater volume of grinding wheel is removed with a 
0.005-in. feed. Thus the average loss per pass was 3 times as 
great for the 0.005-in. cut as for the 0.001-in. cut. Observations 
on the durability of diamond dressers in the shop are usually 
based on the number of dressings. The feed used would 
have an important bearing on the comparisons. A considerable 
difference in wear rate was observed between wet and dry dress- 
ing. The average rate of wear on nine stones when the wheel 
was dressed dry was 20 X 1075 ct./cu in. and was only one fifth as 
great, 4 X 10~*ct./cuin., when the tools were flooded with water 
at the rate of 7200 cc/min. The difference tended to become 
slightly greater as grit size was decreased. Average values ob- 
tained for 36, 60, and 120-grit vitreous wheels were 27, 18, and 
8 for dry dressing and 6, 4, and less than 1 for wet dressing, 
respectively (all X 1075 ct./cu in.). 

2 Another set of factors is introduced by the variations in 
types of grinding wheels. The grinding-wheel factors are grit 
size, grade, type of bond, type of abrasive, and size. The rate 
of wear may be affected as much as tenfold over the range of 
types of grinding wheel. The rate of wear increased as grit size 
increased. A soft grade (Norton 60I5BE) wheel abraded the 
test diamond at the rate of 7;, a medium grade (Norton 60M- 
5BE) at 11; and a hard grade (Norton 60P5BE) at 13, X 1075 
ct./cuin. Of all factors in the wheel, the type of bond had the 
greatest effect on wear rate; for example, vitreous-bonded wheels 
are up to 10 times as abrasive to the diamond as are resin- 
bonded wheels. The effect is less marked for fine grits, soft 
grades, and wet dressing, but is still of the order of fivefold. The 
two types of abrasive aluminum oxide and silicon carbide have 
similar effects on wear rate. That the size of the grinding wheel 
influences wear on the diamond was apparent from the results of 
tests on a 20-in. X 6-in. wheel. The average wear rate of 10 
stones was 19 X 1075 ct./cuin. as opposed to 4 X 1075 ct./cuin. 
obtained ona 12-in. X 1-in. wheel, both dressed wet. The increase 
in wear rate with size of wheel is due, possibly, to the increased 
time of continuous contact of the wearing surfaces, and the 
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TABLE 1 
Tests made by 

Tool nO. eatin oe ee ee Te op 6689 
Initial weight, carats Sal 
Final weight, carats. 2.59 
Loss in weight, carats. bee 0.51 — 
Makerof grindérics dsc scsiews valde ss Cincinnati 
IDV DC jOL STURN) eiclaieitc/alaney oss acdhe lone Centerless 
WWiheelisizoyinies seca -rsvsre ceetetetev et aete 20 X 6 


Wheeligrit siz0s 00. ccne vse s ccd cle 36 


Number Of Gressings). 0... a ee ei 1532 
Passes per dressing.... 2 
Depth; of cut, any. .<s0 +40. 0.0025 
Volume of wheel removed. 14 


Wear rate, ct. X 1075/cu in. 


Mee 20.1 
Wiork' eroumdine: ie Sete ces. Seche es 


Valve stems 


tendency of the area of the wearing surface of the diamond to 
increase, both resulting in higher temperatures on the stone. 
The effect of wheel size was further demonstrated in results of 
shop tests given in Table 1. These results also served to indi- 
cate that the order of magnitude of wear rate under the test 
conditions was the same as under shop conditions. 

Thus it appears that the rate of wear of diamond dressers is 
considerably influenced by the conditions of use. The tests, 
while preliminary, served to establish in a general way the wide 
range of wear rates that may be encountered. The wear rate 
may vary between 1 and 40 X 10~*ct./cu in. depending upon the 
combination of conditions and may be as high as 100 X 1075 
ct./cuin. under particularly severe conditions. For most dressing 
operations the average combination of conditions produces a wear 
rate of 15-25 X 10~* ct./cuin., and this would be a good figure to 
keep in mind, or about 20 X 10-5 carat per dressing, since be- 
tween !/; to 1 cu in. is removed per dressing from medium-sized 
wheels. Consumption of one half of a 1-carat stone should pro- 
duce about 2500 dressings. To produce a reasonable number of 
dressings per tool, the size of the stone must bear some relation 
to the rate of wear for the conditions of use. This does not mean 
that for double the rate of wear, the stone should necessarily be 
doubled in size, but a certain minimum size is required. How- 
ever, factors other than the rate of wear must be considered in 
determining the proper size to be used. 

Another factor which has received very little attention is the 
effect of the shape and size of the stone on the finish imparted to 
the wheel for different classes of grinding wheels and work. Also, 
the size and shape of the stone are important in connection with 
form-grinding. When these effects are known more precisely, 
size and shape specifications can be developed which should lead 
to more efficient and economical use of diamond dressers. Of 
interest in this connection was the observation that in the dressing 
operation the wear on the diamond progresses across the flat in 
the form of a ridge as the diamond is traversed across the wheel, 
starting at the point of first contact with the wheel (see Fig. 5). 
The ridge may be compared with the broad cutting edge of a 
chipping tool. Thus the main working area of diamond is 
roughly the length of the developed ridge, which extends across 
the flat, < the width of the ridge, which usually measures 
less than 1 mm. The remaining areas do, however, maintain 
some rubbing contact. Recognition of this mechanism of cut- 
* ting and wear involved may lead to improved designs of dresser 
tools. 


WEAR TEsTs ON VARIETIES OF INDUSTRIAL DIAMONDS 


Representative samples of the various classes of industrial 
diamonds were tested under controlled standard conditions 
designed in conformity with the foregoing findings to reproduce 
average shop conditions. A l-in-wide by 12 in. to 7-in-diam, 
60-mesh-grit medium-grade vitreous-bonded aluminum-oxide 
grinding wheel was dressed dry; the speed of the wheel was 
maintained at 5000 fpm; and the tools were traversed, with a 
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3.9 0.725 0.688 
3.83 0.720 0.685 
0.07 0.005 0.003 
Mattson Heald eald 
Surface Internal Internal 
20 X 3 31/2 X 2 31/2 X 2 
36 60- 60- 
243 810 516 
2 3 and 1 3 and 1 
0.002 0.003 and0.001 0.003 and 0.001 
184 141 90 
38 3.6 3.3 
Toggle-valve 6-In. gears 6-In. gears 
links 
RIDGE 
DEVELOPED 
ON FLAT 


TOP VIEW OF TOOL 


WEAR TAKES 
PLACE AT 
RIDGE 


~F 


DIAMOND POINT 
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drag angle of 15 deg across the 1-in. face of the wheel at a rate of 
0.0015 in. per revolution for six to eight runs of 50 passes each. 
The amount of wear was measured after each 50 passes, and the 
tool was rotated to a new position for each run. 

The results of the tests on 60 diamonds are listed in Table 2 
on an “‘adjusted” basis and on an “over-all’’ basis. The ‘‘ad- 
justed”’ value of wear rate is derived by excluding values of runs 
in which chipping occurred and in which “‘soft’’ directions hap- 
pened to be struck to give abrupt increases in wear. This latter 
effect was an indication that the variation in hardness of diamond 
influences the performance of dresser tools; for, when the 
orientation of the stone during a run happened to coincide 
roughly with the so-called soft direction, almost invariably 
the wear rate for the run was high and the flat usually had a 
polished or shiny appearance. The soft directions are crystal- 
lographically the same directions that the gem cutter uses in 
lapping on the facets of a cut stone. Analysis of a large number 
of runs on a large number of stones indicated that the rate of wear 
on the soft directions was on the order of 3 times as high as on 
the other or “‘hard”’ directions. 

Incidentally, when a slight drag angle is used in the dressing 
operation, setting of the diamond with a crystallographic axis 
parallel with the axis of the tool would tend to avoid the possi- 
bility of striking soft directions as the tool is rotated to new 
positions; and thus the best results would be obtained. 

Reference to Table 2 will indicate that the results on the 
“adjusted” basis show definitely that many of the lower-quality 
diamonds are as resistant to abrasion by grinding wheels as are 
the higher-quality stones. Taking each class as a whole, there 
were no distinct differences among the classes. ‘The wear rates 
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GROUP 


A. Dodecahedron 


B. Octahedron 


C. Combination Octa-Dodec-Cube 


D. Cube 


1 Adjusted wear rate in carats loss per cubic inch of wheel removed X 1075, 
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TABLE2 WEAR RATE OF DIFFERENT TYPES OF INDUSTRIAL DIAMONDS 


CRYSTAL 
CLASS _NO. 


NONNNNNANINNA FEF 


lst Brown Round 
lst Yellow Shape 
Common Grey 


OPACITY i. fpstlonean meee | 
ve D. T.C. CATEGORY SOG CHIPPING 


lst Grey Round 
lst Colored Round 
lst Colored Round 
Common Coated 

2nd Grey Brown 
2nd River Round 


2nd Yellow Round 
2nd Grey Shape 
lst Grey Round 
2nd Yellow Round 
lst Cleavage 
Dark Brown Round 
Black Rejection 
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Black Bort Stone 
Black Bort Stone 
Black Bort Stone 
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Common Coated 
Common Coated 
Common Coated 
Common Coated 
Common Coated 
Colored Bort 
Black Rejection 
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Common Coated 
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1 
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ae aL lst Brown Shape 23 hip 
1 Colored Cleavage 108 Many chips 
i Cut Gem Stone 18 Chip 
2 lst Grey Shape 3h Crumbled 
2 lst Grey Brown 20 Chips 
2 2nd Grey Brown al None 
5 Bleck Rejection 8 None 
6 Black Rejection : 30 Chip 
5 Common Grey Colored 28 1S 31 Large Chip 
5 Black Bort Stone 36 52 Large Chip 
2 Common Coated -- -- Large Chip 
3 Common Coated 12 18 None 
4 Colored Bort Stone 18 h2 None 
4 Common Coatec 18 10 None 
L, Colored Bort 20 -- None 
6 Colored Bort Stone 13 2k None 
7 Black Bort Stone 9 | 9 | None 
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4 Common Coated 10 37 Crumble 
4 Common Coated 113 US) Chip 
4 Common Coated 30 53 Chipping 
5 Common Coated 13} 38 Crumble 
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5 Black Rejection 26 None 
5 Common Coated 356 Chipping 
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7 Black Rejection 160 Large Chip 
4 78 Chip 
"/ 290 Large Chip 
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2 Over-all wear rates in carats loss per cubic inch of wheel removed X 107-8, 
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adjusted basis, in each group varied from about 10 to 25 and 
averaged 18 X 1075 ct./cuin., with a few scattered values above 
30 X 10-5 ct./cu in. in each group. 

The ‘over-all’ wear rate differs from the adjusted wear rate 
in that it was based on the total amount of diamond lost and 
included loss due to wear in all positions, chipping, and crumbling. 
The difference in rate between over-all and adjusted wear rates 
might be thought of as a measure of the structural failure which 
takes place. Excessive over-all losses, values over 40 X 1075 
ct./cu in. in Table 2, shown by a number of stones in all classes, 
were caused in most cases by chipping and crumbling of localized 
areas. This excessive over-all rate was more prevalent among 
the lower qualities, classes 4 to 11, than among other qualities. 
On the other hand, many of the lower qualities did not exhibit 
chipping or localized increases in wear, and the over-all wear 
rates were essentially equal to the adjusted wear rates. 

In the progress of this research it became more and more appar- 
ent that the chief factor involved in the use of diamonds for wheel- 
dressing purposes is the ability of the stone to withstand fracture. 
Structural weaknesses of diamonds which may be potential 
sources of failure are the same as those for other crystalline ma- 
terials, such as cleavage planes, twinning planes, internal frac- 
tures, incipient cracks, inclusions which act as loci of stress 
concentrations. To this list should be added the interface be- 
tween the interior crystal and the exterior, or later growth, which 
is often of a dissimilar color, characteristic of coated stones. 
Relatively easy parting in the octahedron plane known as cleav- 
age follows from the diamond structure. All that is needed is a 
sufficient force applied in the proper direction to cause cleavage. 
Incipient cracks in the cleavage directions increase the tendency 
to cleave and are potential sources of failure. Large inclusions 
of dissimilar material may be causes of failure, in particular, 
under thermal stresses. Flaws and internal fractures are poten- 
tial starting points of failure. 

What has been said of cleavage planes is also true of twinning 
planes. Diamond twins are common with the twinning plane 
parallel to the octahedron face. These are known as contact or 
spinel twins. In the trade the stones are called macles. Judging 
from the structure of the so-called coated stones, the interface 
between the interior crystal and the coating, or secondary 
growth, may be a plane of weakness. Pieces of the coating are 
likely to be pulled out when the stone is worn to expose the in- 
terior crystal. It would be expected that the position of the 
plane of weakness in the stone would be important. The further 
from the working point that the flaw, or plane of weakness, is 
removed, the less should be its influence. The setting metal and 
method is important in that a stone well supported at all points 
is less apt to fail through a plane of weakness. 

To test for strength, in particular, the effect of weaknesses in 
structure, a large number of runs must be made to wear the stone 
down and expose the planes of weakness; and also, long use may 
develop cracks from continuous shocking. While the laboratory 
test was quite severe in that no coolant was used, the amount of 
abrasive wheel removed was small, and very little of the diamond 
was consumed in a test. In the life of a tool in the shop the 
number of passes is far greater by 50 to 100 times. Hence the 
chances of striking a plane of weakness in a stone are increased 
proportionately. In order to approach the requirement of a 
large number of dressings, tests were made with ten tools repre- 
senting the most important classes of stones on a Cincinnati 
centerless grinder. The conditions-of test were essentially the 
same as those of the standard laboratory test except that a 
coolant was used and the wheel was much larger. Each tool was 
given at least 850 passes over the wheel, removing between 250 
to 270 cu in. of the wheel in a test. The amount of diamond 
consumed in each case was between 0.05 and 0.1 carat which 


amounted to 5 to 20 per cent of the stone, depending on the 
initial size. The results are listed in Table 3. 


TABLE 3 WEAR RATE ON VARIOUS CLASSES OF STONES WHEN 
TESTED ON CINCINNATI CENTERLESS GRINDER USING 20-IN. 
X 6-IN. WHEEL, WET DRESSING 


Wear rate—carat loss X 10-5 per cu in. 
wheel removed 
Cincinnati centerless 
Standard labora- grinder test, 20-in. X 6-in. 
tory test whee 
12-in. X l-in., Computed from Weight- 


Class 


no. wheel volume loss loss check Remarks 

2.2 12.5 11.1 8.4 

4.2 10.3 17.5 16.1 

4.3A 9.8 25.5 19.7 Some chipping 

4.5 10.5 23.4 17.6 Some chipping 

6.7 9.4 23.8 22.3 pew chip at inter- 
ace 

7.4 12.9 15.8 94.27 Large chip — 0.14 
carat 

7.5A 25.8 14.0 30.45 See chip at inter- 
ace 

10.5 21.8 15.5 15.7 

10.6T 24.5 14.7 12.2 Fissure opened 

11.6 24.1 32.4 26.9 Some chipping 

Average 16.2 19.4 18.8 


® Large chips caused by high infeed on last run accounts for high wear 
rate on weight check. If chip of 0.14 carat is allowed for, wear rate is about 


’ Difference between computed and weight check can be accounted for 
with 0.03-carat chip. 


Here again, there was very little difference in wear rate, chip- 
ping excluded, between the types of stones tested. Two out of 
the eight lower-quality stones chipped quite badly. The chip- 
ping occurred when the stones were worn sufficiently to expose 
the core under the coating. Fig. 6 shows the chipping pro- 
gressively. Other coated stones were worn down to the point of 
exposing the plane between the primary and secondary growth 
but did not chip. Aside from failure by chipping it was observed 
that the wear rate did not change as the stones were penetrated 
from their exterior natural skins inward to depths equal to 50 
per cent of the radius. This observation is contrary to the 
opinion of the trade, which has maintained that the natural skin 
of the diamond is harder than its interior. 

Thus it appears that the performance of coated bort stones 
may be influenced by the thickness of the coating. When the 
coating is very thick in relation to the diameter of the stone, a 
large number of dressings at normal wear may be made before 
the interface is penetrated at which point chipping may or may 
not occur. Such chipping as may occur would then reduce the 
life of a stone by only a small amount. When the coating is 

relatively thin, flaking may occur after only a few dressings, but 
the life is not materially shortened since the chipping consumes 
only a small portion of the stone. 

The practicability of bort for dressing purposes was further 
demonstrated by the results of tests which were made by Norton 
Company. These results were especially helpful in providing 
practical confirmation of the laboratory tests, which confirmation 
would have entailed considerable time and effort to obtain. 
Forty-eight tools made with selected bort, similar to classes 6 to 9, 
were used under regular shop conditions on most types of grinders 
in general use for over a year. The tools were used to the limit 
of the setting, no attempt being made to obtain maximum use 
by resetting. Thus the tools were worn down until the stone 
dropped out of the setting. In one half of the tools, the residual 
stone was lost. The results of these tests are compiled in Table 
4 and are analyzed in Table 5. 

The average size of stone used in the tools was 1.28 carats. An 
average of 0.6 carat was consumed which produced an average 
of 2000 dressings. In all, upward of 100,000 dressings were made 
with the forty-eight tools and six are still in use, with about 20 to 
be reset. Fifty per cent of the tools performed as well as first- 
quality stones to give over-all wear rates of less than 40 X 1075 
ct./cuin., and an average of 26 X 10~*ct./cuin. Chipping proba- 
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TABLE 4 RESULTS OF SHOP TESTS ON DIAMOND DRESSERS USING SELECTED BORT BY NORTON CO. 


Weights in cre Weight Loss 
Tool pea a ae No. Carats x 10-5 | Wear2 
No. | Initial | Final| Loss Dressings| per Dressing Rate Rema 
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* Final weight estimated from cavity in tool nib since tool was used to the limit of the setting and the residual stone 
was pulled out and lost. All other values are weights of residual stone. 

2 Wear rate in carats lost X 1075 per cubic inch of wheel removed. Wear per dressing divided by 0.75 cu. in., the 
average amount of wheel removed per dressing as estimated by Norton Company. 
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high average loss of 152 X 10-5 ct./cu in, because the stones broke No. of | Per varie’ Average number 
up into fragments. In all cases a fair number of dressings was Re’ OSIRIS LS ot ieee ee 
obtained before failure occurred. All things considered, the apr obey fragments... ; iG 152 750 


. a Wear rate above 40, indicat- 

coated Congo stones which were tested in the Norton shops per- ing chipplude endtirectass 
formed exceptionally well. Costwise, the comparison is in favor Dee ices rene eee nen aes 14 30 53 1450 

atk : Wear rate below 40, indicat- 
of bort many times over. ing normal wear......... 25 52 26 3000 


Oar ive 


After 637 passes 


After 756 passes 


After 900 pastes 


Fie. 6 Toot Wits CoaTep Strong, SHOWING PROGRESSIVE FRACTURING OF CoATING IN UsgE; AsBoutT X6 


(Test conditions: 


Cincinnati centerless grinder, wet dressing, 20-in. X 6-in., 60M30 Carborundum Company wheel, 0.001 in. feed, 


15-deg drag angle.) 


CoNCLUSIONS 


The order of magnitude of the rate at which diamond wears 
when used to dress grinding wheels was established over the range 
of conditions encountered in dressing operations. A dressing of 
the average size and type of grinding wheel in the shop consumes 
about 0.0002 carat of diamond. A tremendous quantity of 
grinding wheel can be removed with a very small consumption of 
diamond; the ratio is on the order of 1 part diamond to 1,000,000 
parts grinding wheel by weight. The consumption of diamond 
may vary between less than 0.00001 to 0.0004 carat per cu in. of 
wheel removed, depending on the combination of conditions. 
Fine-grained wheels and resin-bonded wheels abrade diamond at 

“a much lower rate than coarse-grained wheels and vitreous- 
bonded wheels. To obtain the best and most economical per- 
formance with diamond dressers, particular attention should be 
given to the selection of the most efficient feed consistent with 
requirements and the adjustment of the coolant in wet dressings, 
since these two factors have the greatest effect on wear. 

Grinding wheels can be dressed satisfactorily with lower- 
quality industrial diamonds. Results of laboratory and practi- 
cal shop tests indicated that many of the so-called lower-quality 
types of stones, of which large quantities are available in better- 
bort and crushing-bort assortments, may perform as well as the 
higher-quality stones of the industrial-serie variety in a majority 
of standard dressing operations. The different varieties and 
types of diamond have essentially the same resistance to wear 
by abrasion of grinding wheels, when all factors’are considered. 

Factors other than the resistance to wear are of greater im- 
portance in determining the suitability of individual stones for 
dressing purposes. Soundness and structural strength and 
the ability of individual stones to withstand chipping and crum- 
bling appear to be the limiting factors in their utilization. 

The selection of sound stones from bort is difficult, because 
of the wide variation and the opacity of bort. The uncertainty 


in performance usually associated with bort arises from this 
difficulty in selection, and this difficulty probably accounts for 
the preference in the trade for the so-called higher qualities 
which can be more readily inspected because they are trans- 
parent. Even so, it should be possible to select sufficient quanti- 
ties of single-coated stones from bort categories to make up 
shortages which may occur in the thirty-two categories of the 
industrial serie. 

Care is necessary in the selection of suitable stones from better 
bort and crushing bort. Solid appearing more or less sym- 
metrical whole stones, free from surface flaws, fissures, and 
granulated areas, octahedron to cubic in shape, may be expected 
to produce consistent results. Aggregated crystal masses, frac- 
tured stones, fragments and cleavages, can be used to dress 
grinding wheels, but their performance may not be entirely 
satisfactory because of a tendency to chip. 
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Measurements of ‘Temperatures in 


Metal Cutting 


By A. O. SCHMIDT,! O. W. BOSTON,? anv W. W. GILBERT? 


Measurements of the heat developed in a metal-cutting 
operation, as well as investigation of cutting and chip 
temperatures by the authors have brought out that the 
chip temperature is uniform at the higher cutting speeds 
when all other conditions are kept constant. The tem- 
perature of the tool increases with the cutting speed, and 
the amount of metal removed before tool failure will 
generally be inversely proportional to the cutting speed. 


REVIEW AND DIscussION OF PREVIOUS INVESTIGATIONS 


N a previous paper (1)4 by the present authors, it was shown 
I that the work done by a metal-cutting tool is transformed 
into heat which is carried away by the tool, workpiece, and 
chips. When the time of cutting is very short, the losses due to 
radiation, convection, and conduction are so small that they can 
be disregarded in most metal-cutting tests in which the calori- 
metric principle is applied. A comparison of the horsepower 
values in drilling, as computed from torque and thrust with those 
computed from the heat generated during cutting, showed that 
the calorimetric method gives results comparable to those of a 
well-calibrated dynamometer, 
When the quantity of heat in the chips only is measured with 
a calorimeter (2), the average temperature increase of the chip re- 
sulting from a cutting operation can be computed by the thermal- 
balance method: 


(AT,) (chip weight) (specific heat) = (A7’) (water equivalent 
of calorimeter, chips, and water) 
where. AT, = average temperature rise of chips, deg F; and 
AT = temperature rise in calorimeter, deg F. 


The foregoing formula does not take into account the change 
of specific heat with increasing temperatures. With the specific- 
heat values of gamma iron, a corrected maximum average chip- 
temperature value was computed for cutting tests in which the 
test-bar material was 8.A.E. 1055. 

Boston and Kraus (3) evaluated cutting temperatures when 
machining an annealed low-carbon steel in a lathe at cutting 
speeds between 20 and 320 fpm. A constant depth of cut of 
0.150 in. and a constant feed of 0.030 in. were used. The cutting 
temperatures were measured with a thermocouple arrangement 
in which one element was the tool and the other the workpiece. 
These were calibrated in a furnace and the temperatures evalu- 
ated in millivolts. The temperatures obtained during cutting 
operations at different speeds are plotted in their paper and 


1 Research Engineer in Charge of Metal Cutting Research, Kearney 
& Trecker Corporation, Milwaukee, Wis. Mem. A.S.M.B. 

2 Professor of Metal Processing, Chairman of Department, Uni- 
versity of Michigan, Ann Arbor, Mich. Mem. A.S.M.E. 

§ Associate Professor, Metal Processing, University of Michigan, 
Ann Arbor, Mich. Mem. A.S.M.E. 

4 Numbers in parentheses refer to the Bibliography at the end of 
the paper. 
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indicate an increase in temperature from 650 to 900 F between 56 
and 115 fpm cutting speed.® For the cutting speeds between 
115 and 320 fpm, the cutting temperatures are reported as vary- 
ing between 900 and 1000 F,, the maximum temperature of 1000 
F being obtained at a cutting speed of 210 fpm. 

Plotted on the same co-ordinates are chip temperatures de- 
termined from temper colors. Very little variation could be 
detected for the cutting speeds between 100 and 320 fpm, and the 
authors state in their conclusions: ‘‘The temperature of the chip 
as determined by temper colors was changed less than 5 per cent 
in the range of cutting speed from 90 to 320 fpm.” 

This finding is supported by the measurements of chip tem- 
peratures with the calorimeter. The temperatures were de- 
termined as almost constant for the higher cutting speeds. This 
is not the case at the lower cutting speeds. Even though the work 
to remove the chip is the same, at a low cutting speed the work- 
piece and tool will conduct a higher percentage of the heat gen- 
erated in the chip because of the slow separation of the chip from 
the workpiece. 

Even though work, which is transformed into heat, must be 
done in metal cutting, it is possible to remove metal at extremely 
low speeds with a relatively low temperature increase in the chip. 
This is because the chip is separated slowly, and the heat in the 
chip is given a chance to flow into the tool and workpiece. It is 
recognized that at the low cutting speeds less heat per unit time 
is generated, although a larger proportion of this smaller amount 
of heat generated at low speeds goes by conduction into the tool 
and workpiece. At higher speeds the work done is nearly pro- 
portional to the increased cutting speed, and most of the heat 
developed is carried away with the chip because of the short time 
the chip is in contact with the workpiece and tool. Assuming 
that the energy or work required to remove a given amount of 
metal is the same regardless of the time involved in removing it, 
it follows that for a longer period of time, that is, at low cutting 
speeds, the heat per unit time is less. At higher cutting speeds 
the time to remove the chip is reduced, and the amount of heat 
developed at any one instant is higher. It is therefore held that 
the temperature of the tool is much higher at higher speeds than 
at low speeds because the work done per unit time is much greater. 
For this reason a tool will generally remove a greater amount of 
metal at low cutting speeds before failure than at high cutting 
speeds. 

The average chip temperature measured at higher cutting 
speeds is more uniform because there is less chance for the heat to 
travel into the tool and workpiece. On the other hand, as soon 
as the temperature at the cutting edge is equal to that in the 
chip, there will be no heat flow from the chip into the tool. 
The higher amount of work which is done per unit time at higher 
cutting speeds and the resulting greater amount of heat is also 
carried by a correspondingly larger weight of chips. 

When the cutting speed is increased, the tool will have to doa 
proportionately greater amount of work. This results in a tem- 
perature increase in the tool material. The heat will not be uni- 
form in the tool, but most intense near the point at which the 
cutting takes place. The zone which is near the cutting edge 


' Bibliography (38), Fig. 12. 
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800 


CHIP TEMPERATURE °F 


100 
CUTTING SPEED FPM. 


150 


Fic. 1 Curves or Cuip TEMPERATURE Vs. CuTTiNne SPEED 


(From diagram by Boston and Kraus, reference 3, Fig. 12. The upper tem- 

perature line was determined by the tool-work thermocouple while the lower 

ones were estimated from the temper colors of the chips. Annealed 0.21C 

steel, depth of cut 0.150 in., feed 0.030 in. per revolution. The dotted line 
indicates average temperature values at the low speeds.) 


will reach the highest temperature and will break down first. 

When using the tool-work thermocouple, Boston and Gilbert 
(4) came to the conclusion that the temperatures measured may 
be affected by changed electrical properties of the metal in the 
test log which is compressed and deformed by the tool point. 
However, this change of electrical properties is brought about 
not only by the changes in the test log but also by those in the 
tool. Because the thermocouple was calibrated with the tool and 
workpiece at the same temperature, the calibration was not made 
under the condition which generally exists in a metal-cutting 
operation. During cutting, tool and workpiece temperatures 
will vary, and at the higher cutting speeds the temperature of the 
tool will be higher than that of the chip. This is a condition 
which can hardly be obtained during the calibration of the tool- 
work thermocouple, and its readings in an actual test will not 
reveal true cutting temperatures. 


MEASUREMENTS OF CHIP TEMPERATURES 


From the investigation carried out on the calorimetric prin- 
ciple, it was found that the cutting force for a particular radial 
rake angle and chip cross section has substantially the same 
value at any cutting speed (5). In the foregoing tests the cut- 
ting speed was varied between 100 and 800 fpm. 

At a cutting speed of 300 fpm, the tool forces were the same‘as 
those at 600 fpm. Since the maximum volume of metal re- 
moved with any test cutter was 0.6 cu in., the cutters could be 
considered sharp. Therefore, any error introduced by a change 
in the cutting edge could be disregarded entirely. However, 
when cutting at 600 fpm, twice as much metal is removed per 
unit time as at 300 fpm and, consequently, the power required 
for cutting is doubled. In computing the chip temperature 
with the formula given, the same value was obtained for both 
cutting speeds. Although twice as much work was done per 
unit time when cutting at a speed of 600 fpm as compared with 
300 fpm, the temperature of the chips remained the same because 
the work was done in removing twice the amount of metal. This 
means there is twice as much chip material to carry away twice 
as much heat. (The amount of heat per unit weight of chips 
remains constant.) 

Most of the work expended in a metal-cutting operation is 
carried away as heat in the chips. The remainder of the heat is 
-absorbed by the tool and workpiece. 

The heat in the tool is of greatest importance because it is the 
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800 


AVERAGE CHIP TEMPERATURE °F. 


400 600 800 
CUTTING SPEED FPM. 


1000 =: 1200 


Fic. 2. AvpraGE CHIP TEMPERATURE IN RELATION TO CUTTING 
SPEED 


(These chip temperatures were measured with a calorimeter when machining 

l1-in-diam test bars of S.A.E, 1055, normalized, with a face mill which had two 

tungsten-carbide tips of 6-deg negative radial rake. Depth of cut 0.125 in., 

feed per tooth 0.008 in. The dotted line indicates estimated temperature 
values at lower cutting speeds.) 


Fie. 3 CrMENTED CarBIDE Curtine Too, Arrex MacuInING 

10.8 Cu In. or Steet at 580 Fpm Curttine Speep, Depts or Cut 
0.150 In., Freep 0.010 In. pER Toor 

(Two passes were taken over the workpiece and the chips are arranged in 

columns in the order in which they were removed. Chip formation is uni- 


form throughout the cut. Abrasive action of chips on blade face has formed 
a shallow crater 0.006 in. deep.) 


accumulation and concentration of intense heat at the cutting 
edge of the tool which hastens failure. Comparing tool life 
with cutting speed, as brought out in numerous investigations 
(6), it is evident that the amount of metal removed before tool 
failure will generally be inversely proportional to the cutting speed 
when the chip cross section is kept constant. This happens de- 
spite the fact that the chip temperature will remain substantially 
the same. Assuming an ideal case in which the tool does not 
wear and thus change the cutting forces, the same volume of 
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Fic. 4 Crmpntep Carpipe Curtinc Toout Arter Macuinine 
10.8 Cu In. or Stent UnprER Conprtions Ipenticat Witn THosE 
In Fic. 3 Excerpt THat THE CuTtrine SPEED Was 2650 Fem 


(Note that at the start the chip formation is similar to Fig. 3 but, due to 
the rapid increase in heat at the cutting edge, wear of the tip has made a 
crater 0.050 in. deep through which the chip formation is affected.) 


too] metal has to dissipate twice the amount of heat per unit 
time when the cutting speed is doubled. For this reason, a tool 
will break down faster at higher cutting speeds, and remove a 
smaller volume of metal than at a lower cutting speed. There 
is only a small zone including the cutting edge of the tool which 
will have a high heat increase with an increase in cutting speed. 
It is in this zone also that most of the friction between the chip 
and the tool occurs (7). 

The temperature of the workpiece is affected mainly by its 
size, its heat conductivity, and specific heat. Other factors 
which also influence the workpiece temperature include cutting 
speed, cross section of chip, tool design, and wear at the cutting 
edge. All these must be taken into consideration. 
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Fig.5 Surrace TEMPERATURE RISE OF WORKPIECE 


(This graph is based on the evaluation of 500 readings taken with an 

nor low-range thermocouple when machining workpieces of S,A.E. 

1055 of identical properties and dimensions at speeds and feeds as indicated 

and a constant depth of cut of 0.125in, The cutter was the same as the one 
used in tests plotted in Fig, 2.) 


However, if all conditions except cutting speed are kept con- 
stant, the following can be stated with regard to the chips, 
tool, and workpiece. 


1 The temperature of chips of uniform cross section remains 
approximately the same in the high range of cutting speeds. 

2 The temperature at the cutting edge of a tool increases 
with the cutting speed. 

3 The temperature of the workpiece is lower at the higher 
cutting speeds. ‘ 
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Charts for Fuselage Torsion Versus 
Control-Surface Flutter 


By W. T. THOMSON,! SAN DIEGO, ‘CALIF. 


Charts for evaluating the allowable unbalance in the 
control-surface product of inertia for tail flutter involving 
fuselage torsion are presented. The basic theory, as- 
sumptions, and limitations are fully discussed. 


NOMENCLATURE 


The following nomenclature is used in the paper: 


6.= torsional displacement of tail, radians 
6, = torsional displacement of wings, radians 
8 = control-surface rotation with respect to fixed surface, 


radians 

semi-chord of tail surface, in. 

location of control-surface hinge axis measured from 
tail-surface mid-chord, as a fraction of tail-surface 
semi-chord 

location of control-surface leading edge, measured from 
tail-surface mid-chord, as a fraction of tail surface 
semi-chord 


xg = distance from hinge line to center of gravity of the con- 
trol surface (positive if CG is behind hinge), in. 
y = span distance measured from torsional axis, in. 
L = control-surface span, in 
J. = mass moment of inertia of wings with, respect to tor- 
sional axis, lb-in. sec? 
J; = mass moment of inertia of tail with respect to torsional 
axis, lb-in. sec? 
- 2 
JeHJrt (*) J,, = effective torsional moment of inertia, 
lb-in. sec? 
Jg = mass moment of inertia of control surface with respect 
to hinge line (rudder or two elevators), lb-in. sec* 
mg = mass of rudder or two elevators, lb-in.~! sec? 
p = mass density of air, lb-in.~* sec? 
v = flight velocity, in. sec—? 
w = circular frequency of flutter, radians sec™+ 
we = natural torsional frequency of fuselage, radian sec~? 
wg = natural circular frequency of control-surface system, 
radians sec” 
g = structural damping coefficient @ times logarithmtic 
. decrement 
INTRODUCTION 


Flutter *of control surfaces may be avoided by balancing. 
The complete balancing of control surfaces, however, is costly in 
weight and often unnecessary in avoiding flutter. 

The study of varying the amount of balance by flutter calcula- 
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tions is long and tedious. To avoid such calculations, charts 
enabling the determination of the allowable static unbalance of 
the control surfaces have been presented by the Army Air Force 
(1)? for flutter modes involving fuselage bending. The purpose 
of this paper is to present a similar type of labor-saving device 
for the flutter modes involving fuselage torsion. 


DEVELOPMENT OF EQuATIONS 


In the following analysis with two degrees of freedom, La- 
grange’s equation 


will be used, where 7 and U are kinetic and potential energies, 
respectively, Q, the generalized force associated with the work 
done in a virtual displacement 6q,, and gq, the generalized co- 
ordinate for the particular mode. 

The kinetic energy of the airplane in the torsional oscillation 
then becomes 


fH 


L 
Le 1 S F 
Bi Gunn hh [(y® + xB)? — (y8)?]dmg 
0 


L 
1 
+= oh 62(dJg — xpg?dmg) 
2/0 


Differentiating with respect to the generalized co-ordinates and 
time, and assuming harmonic motion, the Lagrangian inertia 


forces become 
=n 
06 
op 


The potential energy represents the strain energy stored in the 
structure. Assuming Hooke’s law, it may be written as 


d 


dt 
d 


dt 


—w? (Je + S-xeydmg - 8) 


—w?( J xgydmg - 6 + Jg- 8) 


U = 5 Kee6? 


z 
2 


The Lagrangian elastic forces then become 


oU 

a = Kop0 = we2J 90 : 
Senet pe a [5] 
6 = KggB = wg S apydmg + B 


where the K coefficients have been replaced by the inertia co- 
efficients by considering uncoupled motion in each co-ordinate. 
For terms which remove energy from the system, structural 
damping which is nearly independent of frequency and propor- 
tional to the amplitude must be considered. This can be accom- 
plished by multiplying the elastic coefficients K by (1 + 29). 


Numbers in parentheses refer to the Bibliography at the end of 
the paper. 
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The aerodynamic forces and moments supplying the energy for 
flutter are defined by Theodorsen (2) and Kussner and Schwarz 
(3), for the two-dimensional flow. These equations have been 
organized in a form convenient for flutter analysis by B. Smilg (4). 

The expressions applicable to this analysis are ones for the lift 
per unit span and the torque about the control-surface hinge per 
unit span, and may be written as follows 


bY 


mpb3u” (4m 6 + Agg- a) 


Ly py 


lI 


mpbta? (45 u + Age: a) 


where the aerodynamic coefficients A are complex functions of 
v 

o and have been evaluated in terms of the following coefficients 
a) 

appearing in reference (4) 


Aw = Ly | 
Asg = Lg — L,(¢ — e) 
Age = T, — P(e — e) 
Ape = Tg — Pe(c — e) 


. [7] 
€) 4- B2(c—e)2 


LAGS 


The generalized forces Q are determined from the work associ- 
ated with the virtual displacements as follows 


L L 
iW = di L’5(y@)dy + yh T’5(6)dy 
0 0 
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THE STaBILITy EQUATION 


Substituting the expressions of the previous section into Equa- 
tion [1], the following two equations are obtained 


where 
3 


L 
E = —J6[w? — we2(1 + t9@)] — mpb’w?A og a 
L 
, bL? 
F = —w? i zeydmg — mpb?wA og vy 
0 


L 
bL? 
G = —w? i xeydmg — mpb*w?A gg ry 
0 


H = —Jg[w? — w*g(1 + igg)] — mpb’w?A geb?L 


Dividing each term by —zpb’w*L? and equating the determinant 
to zero, the stability equation becomes 
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Equation [10] is a quadratic in w, resulting in two roots corre- 
sponding to the two flutter modes. The equation is nondimen- 
sional, the three parameters being as follows 


(7) 


SAxpbayLe 


which involves the moment of inertia of the airplane about the 


torsional axis. 
i 
L 


a (pb?) b? A 


which involves the moment of inertia of the control surface about 


its hinge line. 
f E azeydme 
0 L 


(rpb?)bL 


which involves the product of inertia of the control surface about 
the hinge and torsional axis. 

Considerable simplification results from assuming the control- 
surface natural frequency to be zero, under which condition the 
solution of the foregoing equation results in a real and an imagi- 
nary part as follows 


Ww = 
of the structural damping g, necessary to sustain steady-state 
oscillation, may be determined from the imaginary part. 


5 og\. é r anh 
Since (“*) is near unity in the critical-flutter ranges, the value 


CHARTS FOR ALLOWABLE PRopUCT OF INERTIA 


The numerical evaluation of Equation [10] was carried out 
under the following assumptions: 


1 The natural frequency of the control surface at zero air 
speed equals zero (wg = 0). ; 

2 The leading edge of the control surface is at mid-chord 
(e = 0). The effect of varying the leading-edge position is 
indicated by one set of calculations for 40 per cent control surface 
(e = 0.2). 

3 The effective torsional axis of the airplane is along the cen- 
ter line of the fuselage. 

4 The control-surface aerodynamic balance equals 20 per cent 
or 33 per cent (chosen arbitrarily). 

5 Aerodynamic vectors for infinite aspect ratio are applicable. 


The real and imaginary parts of Equation [10], as indicated 
by Equations [11], are tabulated as shown herewith. The 
charts, Figs. 1 to 9, represent the imaginary part of the solu- 
tion as function of aes 

bw 
structural damping, while the abscissa is proportional to the true 


Since the ordinate is proportional to the 
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FUSELAGE TORSION VERSUS CONTROL-SURFACE FLUTTER (wf8 = 0) COMPUTATION FORM 


Hie (Aee)a B+ (Ags)e a 
3 (Ace), (Av): 
ie 
ie 
aly | P ; eae a | is | 
e (Ao), fi (Ase); 
iE | (Aue (Awe | -| : 
¢ (A) (Ave )r 
2 7 
E + (Ag) | + (Avs) | = | e | 
B+ (Ay)e c-tsc-a]=] f | 
-(Ay): 
We es ¢ R olen 
al’ ) J- a he eae e . 
e = 
A(3)9 = Le se Bs 
Pas 


5.12 
715.10 
20.15 
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= +5493 
= 6 Thb 1 


+8538 
+8833 1 


+3972 
~2.3916 1 


wee 0459 
+9812 4 


--1950 
4433 1 


712.4643 
-+2001 14 


--5520 
568242 1 


~20.3737 
+8505 4 


1.3450 
9.5350 4 


49.6014 
5.0100 4 


2.4460 
-16.6400 1 


137-4490 
16.4656 1 


= 95988 
5805 4 


045968 
~-007378 1 


3 wl 0467 04217 
-2.3916 14 --7268 1 --01998 4 
3+33 ~+1950 12.4022 +0372 
40433 4 ©2715 2 =-03703 4 
4.17 -+5520 -20.2736 20342 
“548242 4 1.4701 4 --0487 2 
6.25 1.3450 49.4170 20276 
945350 1 6.0244 4 2079 4 


10.0 | 


=2 4460 
-16.6400 1 


7137-1474 
(18.2358 1 


20184 
01390 1 


-27.689 
3.0608 1 


71.3450 
9.5350 1 


45.1814 
58428 1 


air speed, the charts represent a series of flutter-stability curves. 
The real part of the solution presented in the tabulation may be 
used to determine the flutter frequency if this information is 
desired. All computations were carried out in the tabular form 
shown. 


Errect or Finire-ContrRoLt-SuRFACE FREQUENCY 


Examination of Equation [10] indicates that the effect of finite- 
control-surface frequency may be readily determined from the 
charts presented by assuming g = 0, and using 


ee Fi 
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» Due to neglecting the friction in the control-surface system, the 


results will be on the conservative side. 

Illustrative Example. Determine the allowable rudder product 
of inertia from the following data: 

Average rudder chord = 22 in. 

Average fin chord = 22 in. (i.e., e = 0) 


c—€ 
= 0.20 


Aerodynamic balance 
—e 


Moment of inertia of rudder Jg = 3.0 lb-in. sec? 

Rudder span = fin span = LZ = 80 in. 

Moment of inertia of tail about torsional axis = 1500 lb-in. sec? 

Moment of inertia of wings about torsional axis = 50,000 Ib-in. 
sec? 


6 
Ratio of torsional deflections ‘@) = (0:20 


Torsional natural frequency of airplane = 1000 cycles per min 

Maximum flight speed of airplane = 500 mph 

The calculation of the nondimensional parameters are as 
follows: ; 


mpb? = 1(8.45 X 10-8)22% = 0.0001285 (10,000 ft altitude 
assumed) 2 
Jo = 1500 + (0.20)? 50,000 = 3500 lb-in. sec? 
0 0. A 
Ag = sos = 53.2 X .03 = 1.60 (3 per cent damp- 


~ 0.0001285 x 803 
ing assumed) 
3.0 


B= = 0.603 
0.0001285 x 80 X 22? 


500 x 5 x 12 


From the stability curve for B = 0.6 find P = 0.18 corresponding 
to Ag = 1.60. The allowable product of inertia is then 


JS xeydmg = 0.18 X 0.0001285 X 22 X 80? = 32.6 lb-in. sec? 


From the real part of the solution R = 0.2 corresponds to the 
foregoing values. The frequency ratio then becomes 


we R 0.2 
=) exrql) Line cma ea 
(=) " A J Pe aig ee 


and the error in neglecting this correction is less than 1 per cent. 


SUMMARY AND CONCLUSIONS 


Control-surface flutter involving fuselage torsion has been 
investigated in terms of nondimensional parameters and presented 
in chart form. The results are applicable to conventional tail 
surfaces with 20 or 33 per cent aerodynamic balance. The effect 
of fin or stabilizer bending was not included in the investigation 
since this mode is generally covered in the three-dimensional 
analysis involving fuselage bending. The real part of the solu- 
tion was included in the tabulation as an indication of the fre- 
) which was assumed to be unity. 

w 

Several conclusions may be drawn by examination of the 
charts, as follows: 

1 The stability of the control surface for this form of flutter is 
mainly a function of its product of inertia and aerodynamic 
balance. 

2 For a given product of inertia, the flutter speed will be 
higher for a control surface with larger aerodynamic balance. 

3 Large aerodynamic balance increases the explosiveness of 
flutter, as indicated by higher peaks of the stability curves. 


quency ratio ( 
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4 Reducing the control-surface chord for the same aero- 
dynamic balance increases the flutter speed and its explosiveness, 
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as indicated by the shift of the stability peak to higher values of 
v 

— and Ag. 

5 and Ag 


5 The effect of the finite-control-surface frequency is to shift 
the stability curve to a lower value of B, the main effect of which 
is to decrease the stability peaks to lower values. The flutter 
speed is also lowered slightly. 

The charts presented are intended to aid in the preliminary 
determination of the required product of inertia. They also 
enable a study to be made of the effect of varying the design 
parameters. 
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Thermal Accommodation Coefficients 


By M. L. WIEDMANN? AND P. R. TRUMPLER? 


This investigation deals with the conduction of heat 
_through a gas at low pressure. The theoretical work of 
M. Knudsen is extended to cover the case of conduction 
between two infinite concentric cylindrical surfaces of 
different properties. Experiments to determine the 
accommodation coefficients of air on flat black lacquer and 
nine metal surfaces indicate that the values lie between 


0.87 and 0.97. 


NOMENCLATURE 


The following nomenclature is used in the paper: 


A = surface area, sq ft 
a = radiant emissivity 
E, = energy of a gas stream in equilibrium with surface 1 
E,’ = energy ofa gas stream after collision at surface 1 
eet 
2 
L = length of cylindrical test section, ft 
M = molecular weight of the gas (Oy is 32) 
p = gas pressure, mm mercury 
Q, = net conduction heat transfer, Btu per hr 
Qrz = net radiation heat transfer, Btu per hr 
r = radius, ft 
T = gas temperature, deg F abs 
T, = temperature of surface 1, deg F abs 
a = accommodation coefficient 
Y= : = ratio of specific heats of the gas. 
p = length of a path between inner and outer surfaces, ft 
o = Boltzmann constant = 0.172 X 10-8 Btu per hr per 
sq ft per (deg F)# 
B, 6,0 = angles 


Subscript 1 refers to inner surface, subscript 2 to outer surface. 


INTRODUCTION 


In a sense there are only two methods by which heat may be 
transferred between solids separated by a gas, i.e., radiation and 
molecular collision. Collisions may take place between the gas 
molecules themselves or between gas molecules and the surface of 
the containing vessel. The latter type of collision is the subject 
of this investigation. 

At ordinary pressures the effect of surface collisions is lost in 
the large number of collisions involving only gas molecules. But 
if the pressure is reduced far enough there will be only a compara- 
tively small number of gas molecules in the containing vessel, and 
most of the collisions will be made at the bounding surfaces. 


1 This paper summarizes a dissertation submitted by one of the 
authors’ to Yale University in partial fulfillment of requirements for 
the degree of Doctor of Philosophy. 

2 Associate Professor of Mechanical Engineering, Yale University, 
New Haven, Conn. : 

3 Development Engineer, M. W. Kellogg Company, New York, 
N.Y. Jun. A.S.M.E. 

Contributed by the Heat Transfer Division and presented at the 
Annual Meeting, New York, N. Y., Nov. 29-Dec. 3, 1943, of Tux 
AMERICAN SOCIETY OF MECHANICAL ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be un- 
derstood as individual expressions of their authors and not those of 
the Society. 
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These collisions will then exert a controlling influence on the heat 
transfer. 

Kinetic theory predicts that the conductivity of a perfect gas 
far from a solid surface is independent of the pressure. If certain 
dimensions of the containing vessel are of the order of magnitude 
of the mean free path,‘ the heat conduction is a function of the 
pressure and the surface characteristics. For the special case of 
gas pressure so low that collisions between gas molecules are in- 
frequent compared with surface collisions, the heat conduction is 
proportional to the pressure (1, 2, 3).5 The theoretical treatment 
to be presented applies to a perfect gas at this low pressure, here- 
after called a “rarefied gas,” and the experimental work was con- 
ducted with air in or near this pressure region. 

Air at 32 F and 0.000001 atm pressure has a mean free path of 
2.3 in. (4). It is apparent that the results here presented are ap- 
plicable primarily to problems of heat conduction through a gas 
at very low pressure or in a very small space. 

The characteristic properties of a surface which determine the 
radiation are the emissivity and the absorptivity. In heat trans- 
fer by surface collisions the characteristic property is the accom- 
modation coefficient, which expresses the extent to which a par- 
ticle stream impinging upon a solid surface comes to thermal equi- 
librium with the solid. 

If a particle stream approaches a surface at energy EH,’ and 
leaves at energy H,’, and if the energy of the stream at thermal 
equilibrium with the solid would be Z;, then the accommodation 


coefficient is defined as 
RB he 
Limit eae [1] 
E,/—E, _E2’ — Ey 


If only kinetic-energy changes are considered and the particle 
stream has a maxwellian velocity distribution before and after 
collision, the accommodation coefficient may be expressed in 
terms of temperatures 


: IN he 
Si UAV |) een dent one as 2 
i Limit | I =i | 


Values of accommodation coefficient reported in the literature 
range from less than 0.1 to unity. 


fs 


REsuLtTs or Tests 


The equation for heat conduction between two infinite con- 
centric cylinders separated by a rarefied gas is 


; 1 poe Pal 
eI App ee) Bere Ee) Bl 
6 y—1 S/T 1 ( ) 1 
= {Sy} \) oh 
T2 \ ae ay 
where 
@ = ‘heat transfer by conduction, per unit cylinder length, 
L Btu hr-!ft-} 
y = ~ = ratio of specific heats of the gas 
Cy 
p = gas pressure, mm mercury 


4 The mean free path may be defined as the average distance a gas 
molecule travels between collisions with other gas molecules. 

5 Numbers in parentheses refer to the Bibliography at the end of 
the paper. 
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T = gas temperature at place where the pressure is p, deg 
FA 
M = molecular weight of the gas (Oz is 32) 
T;,T, = temperatures of inner and outer cylinders, respec- 
tively, deg FA 
r= radius, ft 
accommodation coefficient 


R 
ll 


Subscripts 1 and 2 refer to inner and outer cylindrical surfaces 
respectively. 

Experiments were made to determine the accommodation co- 
efficients and emissivities of aluminum, cast iron, and bronze. 
Polished, machined, and etched surfaces of each metal were 


tested. The values of Q- are plotted against pressure for 


Tells 
nine metallic surfaces and a reference surface of flat black lacquer 
in Figs. 1-3, inclusive. Maximum and minimum values of ac- 

QV T 
p(T: = Jf) 
as a function of pressure p, Fig. 4, and are tabulated with the em- 
issivities in Table 1. 


commodation coefficients were taken from a plot of 


TABLE 1 SUMMARY OF RESULTS 


Accommodation 
coefficient ——. 


Surface Description Emissivity Minimum Maximum 
x Flat black lacquer on bronze 0.932 0.881 0.894 
A Bronze, polished 0.103 0.91 0.94 
F Bronze, machined 0.110 0.89 0.93 
D Bronze, etched 0.192 0.93 0.95 
B Cast iron, polished 0.184 0.87 0.93 
E Cast iron, machined 0.391 0.87 0.88 
G Cast iron, etched 0.697 0.89 0.96 
(& Aluminum, polished 0.221 0.87 0.95 
J Aluminum, machined 0.0909 0.95 0.97 
H Aluminum, etched 0.775 (0.89) (0.97) 


Two general conclusions may be drawn as follows: 


1 The values of all measured accommodation coefficients lie 
between 0.87 and 0.97. 

2 The effect of different technical surfaces on the accommo- 
dation coefficient, if such an effect exists, is small and is not shown 
by the present tests. 


In evaluating the experimental results the abnormal behavior 
of two of the surfaces should be noted. The emissivities of the 
etched aluminum dropped from 0.833 to 0.753 during the tests, 
and the value of the accommodation coefficient is therefore 
unreliable. In testing the etched cast-iron surface G it was ob- 
served that the air pressure in the sealed system decreased ap- 
preciably in the course of each conduction test. The emissivity 
remained unchanged. 

It is undoubtedly true that all the metallic surfaces were oxi- 
dized to some extent. The appearance of each surface was the 
same before and after the tests, and, except for the etched alu- 
minum, the emissivity was constant. But even if the oxide coat- 
ings were stable their influence on the properties of the metallic 
surfaces must be considered in interpreting the data. 

The values of accommodation coefficient found by the authors 
are of the same magnitude as those given by Knudsen (1) for 
oxygen in contact with polished platinum, platinum coated 
lightly and then heavily with platinum black. His values are 
0.808, 0.910, and 0.934, respectively. 


THEORY 


The apparatus employed by the authors for determining ac- 
commodation coefficient requires the heat-transfer equations for 
radiation and for rarefied-gas conduction between two infinite 
concentric cylinders having different surface properties. The ra- 
diation equation, assuming diffuse reflection, is well known (5) as 


JANUARY, 1946 


REFERENCE SURFACE 
BRONZE, POLISHEO 
BRONZE, MACHINED 
BRONZE, ETCHED 


CONDUCTION DATA 
SS FOR 
BRONZE SURFACES 


ce) 10 3 20 26 
PRESSURE, MM H6 X 10 


Fig. 1 Conpuction rigs FoR Bronze SURFACES 


REFERENCE SURFACE 
ALUMINUM, POLISHED 
ALUMINUM, MACHINED 
ALUMINUM, ETCHED 


10 20 268 
PRESSURE, MM Ho X 10° = 


Fic. 2 Conpuction Data ror ALUMINUM SURFACES 


REFERENCE 


CAST IRON, 
CAST IRON, 


CAST IRON, 


SURFACE 
POLISHED 
MACHINED 
ETCHED 


CONDUCTION DATA 
FOR 
CAST IRON SURFACES 


20 26 


PRESSURE, MM He X 10° 


Fie. 3 Conpuction Data ror Cast Iron SurFaces 
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FLAT BLACK 

BRONZE , POLISHED 
BRONZE, MACHINED 
BRONZE, ETCHED 
CAST IRON, POLISHED 


CAST IRON, MACHINED 
CAST IRON, ETCHED 
ALUMINUM, POLISHED 
ALUMINUM, MACHINED 
ALUMINUM, ETCHED 


Fie. 5 


SLOPE 
OF 


particles comes from the inner surface. The re- 
mainder (1 — F) come from the outer surface at 


CONDUCTION CURVES 


energy H,’; therefore 


4 6 8 te} 12 14 16. 8, 
PRESSURE P, MM H, X10 


Fie. 4 Storr or ConpucTIon CuRVES 


oAi(T14 =. T 4) 
Qr = aan WA AAT a : % pao (4] 


= heat transfer by radiation, per unit cylinder length, Btu 
hr-!ft-! 

area of cylindrical surface, sq ft 

= radiant emissivity 

= Boltzmann constant = 0.172 x 107° 


ll 


aa mem | 


Although Equation [5] applies only to gray bodies, its use with 
metallic surfaces is justified as a means of separating radiation 
from conduction in the authors’ experiments. Hmissivities cal- 
culated from Equation [5] are not accurate if the surfaces differ 
appreciably from those of gray bodies. 

The equation for the conduction through a rarefied gas (almost 
all particle collisions at the solid surfaces) may be derived by ex- 
tending the work of M. Knudsen (1), the validity of which has 
been fairly well established (6, 7, 8, 9). His equation for maxi- 
mum conduction between two concentric infinite cylinders (ac- 
commodation coefficients unity) is 

Q Mae oh (Cit 2 1) 


~* = 3600 r. ee 
L ‘y—1 /MT 


With reference to Fig. 5, consider a large number of gas mole- 
cules striking a unit area of surface 1 per unit time having a total 
energy E;’ before, and H,’ after collision. Then the energy trans- 
ferred from surface 1 is H;’— HE,’ per unit area per unit time 
and the accommodation coefficient is 


After collision with surface 2, a particle stream has energy 
E,’. The same particle stream does not, however, approach 


the outer surface with energy Fy’, for only a fraction F of the © 


20 22 24 


— ERE ae OUR OVI od Oe 


an = Sion SSR af (8] 
26 28 4«4OFr 
ie = F(E,' — E;") (9] 
' FE, + 1—F)E)'— EB, 
From Equations [7] and [9] 
BE, — £ 
BY! — Ee = - , Pele. <. [10] 


The symbol F is defined as the fraction of particles leaving the 
outer surface taking direct paths to the inner. Its value can be 
calculated on the basis of Knudsen’s cosine law (7), which may 
be stated as follows: The probability that a gas molecule leaves a 
surface on a path making an angle 6 with the normal to the sur- 
face is proportional to cos 6.8 

For a given path between the surfaces if 6 and 6 represent the 
angles made by the path with the normals to the surfaces at the 
points of contact, and if p is the path length, then 


1 5 
p a2 f eaten go a2 eta [11] 
w SA p re 


Returning to Equation [10], since #, — EH, represents maxi- 
mum energy transfer as given by Equation [6], the general equa- 
tion for heat conduction between two infinite concentric cylinders 
separated by a rarefied gas is ‘ 


1 p(T, — T 1 
a: = 3600n 2 ; ar =F = 1121 
| Pee 
T2 a2 ay) 
Using M = 29.0and y = 1.404 for air 
T, — T. 1 
(*). = 3980 r, 2 a 2 - : 5+ [18] 
air Th Tr (2 st 1) beats 
T2 \ Oe ay 


In the experiments, 7 varied from 0.1040 to 0.1042, m = 
0.1094, and L = 0.667. For this case 


p(T; — T2) 1 


6 This law is analogous to Lambert’s cosine law for diffuse reflec- 
tion of radiation. 
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Although Equation [14] should be used only with small values 
of T,; — Ts, Knudsen’s experiments show (1) that the maximum 
value of 60 F in the present experiments is well within the range 
of application. 


EXPERIMENTAL PROCEDURE 
The apparatus generally used for the measurement of accom- 
modation coefficients (1, 2, 10) is not well adapted to an investi- 
gation of various types of surfaces. This objection is overcome 
in the new design, shown in Fig. 6, which has the additional ad- 
vantage of a center, or test cylinder of simple construction which 
can be easily changed without duplication of test equipment. 
With the exception of the disk washers at the ends of the test 
cylinder, which were designed for a press fit, all other parts were 

used interchangeably for all ten inner cylinders. 
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WIAA SSSA 


RS 
== 
AA 


Sg 
ASSEMBLY OF TEST UNIT SCALE :HALF SIZE 


Fic. 6 AsspmBiy or Test Unit 


Heat supplied by the central heating coil flows outward be- 
tween the two test surfaces approximately 2.50 and 2.62 in. 
diam. The test length is 8.00 in. Energy input to the auxiliary 
coil at each end is regulated to maintain a very small temperature 
differential between the disk washer and the inner face of the end 
plug. Therefore, substantially all the heat input to the central 
coil flows through the test surface. 

The outer test surface was used as a reference and was given a 
double coat of flat black paint in order to make it nearly “gray” 
in response to radiation. Its properties were determined by test- 
ing a similarly prepared surface on an inner cylinder. The emis- 
sivity was 0.932 and the accommodation coefficient 0.89. The 
five radiation tests on this surface for temperature differences 
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from 40 to 80 F checked the Stefan-Boltzmann law with an av- 
erage deviation of 0.25 per cent. To prevent damage to the ref- 
erence surface while changing test. cylinders, a special fixture, 


not shown in the illustration, was designed. 


Energy to the heating coils was supplied by a direct-current 
generator driven by a synchronous motor. Heat loss from the 
central coil leads was experimentally measured and a correction 
established. ; 

The test section was evacuated by a Cenco Megavac pump, and 
the pressure was measured by a McLeod gage with a dry-ice trap. 
The system included a phosphorus pentoxide drier. The pressure 
was corrected for thermal diffusion between air at water tempera- 
ture entering the test casing and air between the test cylinders. 
No correction was made for pressure drop in the piping from gage 
to casing because the correction is less than 4 per cent at the low- 
est conduction test pressures. The diffusion correction is actu- 
ally applied by substituting 7» for gas temperature 7 in Equa- 
tion [3]. 

Temperatures inside the test unit were measured with cali- 
brated copper-constantan thermocouples of 24 B&S gage wire, 
with are-welded junctions. A Leeds and Northrup Type K poten- 
tiometer served to measure the temperature difference between 
the test cylinders, and a galvanometer sensitive to 0.05 deg F was 
used for the end-loss couples. : 

The thermocouple installations should be noted. Contact 
between the couple junction and the inner test cylinder was es- 
tablished mechanically by a set screw, Fig. 9(c). The junction 
on the reference cylinder, being electrically insulated from the 
casing, provides poor thermal contact and therefore conduction 
through the leads should be held to a minimum. The plug con- 
taining the element is shown in Fig. 7, its assembly in the cylin- 
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drical shell extending from the casing in Fig. 9(a), and the posi- 
tion of the shell on the casing in Fig. 10. Installation of the 
end-loss couples is shown in Fig. 9(d), and Fig. 8 shows them 
connected to the end plugs but detached from the disk washers. - 

To draw leads out of the test section, stainless-steel bolts, 
Figs. 8 and 9(b), were made part of each thermocouple circuit. 
An emf correction for the copper-to-steel junctions in each cir- 
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THERMOCOUPLE CONNECTIONS 


Casine or Test UNIT 


(C) TEST CYLINDER COUPLE 


cuit was determined by joining the bolt heads with 
copper wire and operating the equipment under 
previously established conditions. The correction 
in no instance exceeded 2 deg. 


INNER CYLINDER 
SURFACE NEAR 
HEATING COIL 


TC.WIRES Preparation of Test Surfaces. The test cylin- 


COLLAR . 3 
ders were machined from castings made by com- 


mercial foundries in Connecticut. The bronze was 
an alloy of 85 per cent copper, 5 per cent each of 
tin, lead, and zine by weight; the cast iron was 
actually a semisteel, since steel scrap had been 
added; the aluminum was alloyed with traces of 


INNER SURFACE 


CERAMIC TUBE 


ENO OF TEST 
CYLINDER 


(0) END LOSS COUPLES 


other elements to provide free machining and cor- 
rosion resistance. All castings of one metal came 
from the same lot. Only cylinders without flaws 
on the finished test surface were used. 

Test cylinders were machined without lubrication. 
Before the final lathe cuts were made the cylinder 
and the plugs which supported it between centers 
were removed from the lathe and cleaned with 
soap solution, dioxane, and distilled water. After 
this cleaning the cylinder was not touched by 
hand. No further cleaning of any kind was given 
the machined or polished surface. 

Polishing, except for cast iron, was done on a lathe with five 
grades of polishing paper ending with 4/0. Because of metal por- 
osity cast iron was finished only to 1/0 paper. On aluminum the 
finer papers produced a dull finish on which the many fine surface 
pits were entirely covered up. Polishing at various speeds and 
pressures was tried without obtaining a brighter finish. The 
machined aluminum, on the other hand, was bright but clearly 
showed the porosity of the casting. The fact that the emissivity 
of the machined surface was found to be less than that of the 
polished was not entirely unexpected. 

Etching was done with nitric acid, reducing the over-all diame- 
ter of a cylinder 0.004 to 0.005 in. Upon removal from the acid 
bath, the cylinders were washed and scrubbed in tap water, a 
dilute ammonium hydroxide solution, more tap water, and finally 
distilled water. 


INSULATOR 


TRANSACTIONS OF THE A.8.M.E. JANUARY, 1946 


62 


SrianGe Gl’ ee 929°? T6e* °6 40° TT 928° 4° TBS L°9€S $°06 9-9 

ne zB °E ot0’s €98°2 T6E° OTe*S t's Lz9° T° 185 9°SES 1S 6-4 

(EN CHinG 66L°7 €oL*e Tée* 601° 6€°L 955° 1°78S G°9ES yO Yat 

gIe ALOT 190°¢ Siu hard T6E° zS0°€ 9€°9 ost7* 9° BLS 9°9€S as €-a 

€0° SEC°S L8€°* Stet S7°S Cie o°SLS T°9ES 60°0 Cr 

TS0° LBE*E N6E* ZONE IL°9 L0S* 6°26S 9° LES 10) 00S2°T Ta 

O€T $8°TT_ AEN. O€7*T céT® 89°38 4L°OT 808° 1°85 O° Les ba) a 

102 60°% 8c6°6 T69°T c6T° 009°7 T8°L 685° 9° 265 1° 9ES 0°02 $-d 

9ez €9°% B26°S $9S°T one ZOE °E 65°9 T0S* L°68S T° LES 9°IT q-d 

NEe Ov7°T 086°% 6S7°T do) ted Bee °s €9°S ote 0°98S 8°9€S 0°9 €-d 

‘ 0zo° 89e°T Tél’ upg Tey OTE’ T°08S e°9ES $0°0 e-d 

TiOe 61S °Z T6T° 295 °% T8°s TINIE 8°19 B°SES 90°0 O82 *T Tee 

8eT ei 68° TI9°T ec 68°38 68°OT 918° 9° 08S 6° TES 0°26 Ne) 

S6T 69°47 9€°OT OL Th tee BEL°t T6°L 665° 9°85 9°2Es 0°te 99 

0% B7°t 689° Lever eee 108°z 60°9 09t7* z°9LS 9°ZES Sage $-0 

g0z Ge°L /Bint 2 6L2°T Tee" 806°T c0°s OBe° €° TLS 9° TES 0°9 Vint) 

7Z0° g0z°T 612° STe°T To*" €0€° 8°69S 8° TES L0°0 €-9 

TS0° Clad (Gage LST°% YES tov 9° 76S O°€eEes 60°0O ot) 

AGUA Be°S 00%7° 47° €6S G° ES 90°0 06%e°T itd 

SS0° © 919°T TST’ cE9°T 69°77 ste° 0°76 O° LES 0 ya 

26T 92°7 T19°6 81S °T 18T° NEES 7 €9°L 895° 1° T6S 0° 6S (ANKA o-d 

€lz 0¢°€ 918°9 OeT°T 48T° LUTE BL°9 470° T°33S 9°BES O°ST Hd 

T1z Cec 109° SEA 1E 181° 9LL°% OT’9 Sou" $°L8S S° LES LORE €-d 

0% ST°T €IL°2 889 °T 413T° €8E°% 99°S 9° 26S g° Les os ea 

= 190° o%70°e LgT° L50°% 92°S 6° 709 T°BEs 1530) 0672 °T Tai 

Z0e OL°€ TLL°3 4L8°0 €0T° THE 23°9 S0S° 9° 06S 9°9€S $°8T Vv 

612 €9°% 618°S $z8°0 €0T° 0€S*°e €3°S UKSUIS e° 88S 6°9ES 0°eT 6-¥ 

0S2 0°? ZEG°E 80L°0 €0T° 6S8°T TO°S TLE 1° 8S B° LES 2s Y 

Te? 20°T e6°T 989°0 €0T° 642° Tr? t0€° 6°6L5 z°9ES Ty €-V 

190° 99T°T r0) the 48T°T 00°47 962° 8°909 6°9€S 60°0 z-¥ 

8€0° T9L°T 4OT* SLT°T LO°E 962° %°909 T°S€s 90°0 S6%e°T Tv 

-T6°0O 61° €% TO°Ze Woeey/ cE6° LS°€T TL Er 066° 4°69 T° LES 005 ‘Sz OT-x 

TAAL TE*2T 83°72 Bee "9 Acton €9°ET LL°€T 066° 8°7%as 0°9€S “TEE 6-X 

BLT Suhel 9€8°6 6L0°L Zeb" 196°6 9L°TT L18° e°885 8°9ES 8° 9-X 

98T €9°€ 699°L 9eL°9 c€6° €L6°3 OT°TT 403° 1°985 LES S°6T LX 

Clie 65°% onT’s TELA9) 7E6° S0€°8 €L°OT yee 0°L8S €°SES OuAL 9-X 

9¢e c9°T CUES 087°9 e€6° Ost? L ST°OT HELL o°S8S 6° LES Capt, $-X 

zLO* T6T°CT Te6° Stee GO°ET 9€6° T°819 €°9ES 60°0 1-X 

T70° 619° Th° 199°S 6L°8 TI9° L°sls S°9ES T°O c=x 

8t0° 99L°9 BE6° 08l°9 €9°6 One 8°S8S 6°9ES T°0 e-X 

190° 816°8 626° 96°83 90° TT 808° - 8°86S €°SES T°0 eT-X 

T90°0 T9L°S Z€6°0 6LL°8 96°OT T08°0 6° L6S B°SES T°0 S0S2°T Tx 

(¢1 = T)d Cp Ty eee she Te s4yem S}TOA «dure Yde Vie OT x Soteey On 
b) eon aD ceo hi q H wm aT 

ag Lo a Ae g : Seog Ty AN ainssetg aa 


VIVGC LSHL ¢ WTaVL 


63 


WIEDMANN, TRUMPLER—THERMAL ACCOMMODATION COEFFICIENTS 


I uy YE Jo uopzenzonTd aay 
gz se yorum se G UF uoTyenqonTd L 4 TH 
4804 sod ¥5 qnoge ‘pearesqo doap einssetg 966 *1°E-D 
moy sed ZT yRoge espr eunssetg G66 fd 
% UT $6°T JO uotzenzonTd a 
I ut $I go uotzenzonTd g-0 
. pe}Tuo seyseMm BOTT T-0 
*sauout 00’S = 1 fseuouT Sz7TE°T = ex :syreueYy 
86T 96°) witt 29L"°0 6060° g0T’? enh €56° 41°06S T° LES 0°Se Bg-t 
g6T G6°7 ge*tt §=6$SL°0 6060° €S0°" Leek 055° S°06S g° Les 0°se or 
gTe 68° STe°3 S69°0 6060° ZOT°E 97°9 os* 8°98S g° Les O°sTt c-f 
Hoc 89° LEetT’s 1S9°0 6060° ~ 0S2°2 0S °S 6017" €°€3s T°9ES S$°Tt ele 
gee Clee TOE *S 9TS°0O 6060° T6T°T WO) L6c* 1°SLS 6° LES 0°9 SJe 
T70° CAME TE T160° ante © 68°€ 682° 0° 609 B° es 90°0 Cat 
890° 850°T 8060" gL0°T Ts°€ €8e° L°309 B°LEes eo) S6tie°T He 
sto" $99°S €SL° 619°S SL°s 619° 4°98S L°S€sS T°0 BL-H 
gt10° 9LL°S .8SL° 06L°S 18°38 $s9° 0°98S 4° 1ES T°O L-H 
OTe Te°s $8°OT S67°S SLL° ¢L9°3 €3°Ot Tos” 0° 78S B°SEes 3°e 9-H 
ote Bl°e €9¢°8 c6L°S SLL° eles €S°OC osl’ 6°S8S €°SES g°LT S-H 
OT OL’? 66L°S L9°S SLL* g7e* lL 16°6 6€L* T7186 etES IUCAE 4-H 
Tee 0S°T 605 °€ 6€2°9 GLL* L9e°L 06°6 nEL* 0° 68S 6°NES L°9 €-H 
a 190° €Lo°L e6L° T60°L gl°6 Gel 8°S6S 2° Les TT°0 e-H 
iOS 629°7 ces" T19°7 138°L 765° 1°SLS L°9€S EEO S6Ne°T T-H 
TiO’ 9¢8°7 L69° ges" To’s 409° 8°T8s GES 60°0 [K) 
O6T ee°s 29°66 meet 169° O€T*L gL°6 6eL° Z°BLS 9°SES S° Le ° [ea) 
ote €6°% L10°9 gle°" L69° T19°9 Albeo) SOL* e°€8S L°S€s Oyae Sr) 
12 St°? 989°" S9T°S L69° ges °9 €°6 ooL* 18S B° Ces 9°6 (f=) 
6€2% OTT 0S °% Bers L69° e@Lt’9 60°6 619° 6°98S €°7ES 9° €-) 
TS0° 079°9 €69° $59°9 1°6 SoL* 6° 96S O° Tes 60°0 (ed) 
+t0° L50°S ZOL* oL0°S 6T°S 619° €°€9s €°v7es T’0 oste*T Gr?) 
€LT 65°9 TLE oLO°T Ope? 6L1°9 90°6 739° 7°96 Z°SEes O° 8€ Fd 
fT 98°¢ $l0*6 685 °O OTT’ BSE S3°9 gTS° 6° L8S O°9€S 0°0¢ S-a 
Liz 09°% 9L9°S €S8°0O OTL’ 915°? LL’°S gen 0°S8S 9° TES o°et ; ted 
€€e GE°T 9*18° 9¢8°0 Ome 099°T 69°7 qSE° S° 18S B°SES B°S €-d 
8t0° T18°0O s0T’" S78°0 Oe °€ 0Sz* €°€8S €°res 0) c-d 
850° TSe°T ager g9e°T 60° OTE*O 6°€09 L°ses 60°O 00S¢°T Td 
(C= tn) Ged 
z z i) a) 
AO) LA oO Uy Asaieun amssetg 


(panuzuop) 2 ATAVL 


64 TRANSACTIONS OF THE A.S.M.E. 


Method of Testing. All metallic surfaces subjected to low pres- 
sures were cleaned so as to be free of volatile impurities. Degas- 
sing was done at test-cylinder temperatures not exceeding 160 F. 
During radiation tests the vacuum pump was run continuously, 
but in conduction tests the system was closed at the desired pres- 
sure. The pressure rise for an 8- to 12-hr test was usually negli- 
gible. The one exception was cylinder Ff, machined bronze, for 
which the pressure increased as much as 1 per cent per hr. 

For each surface at least two radiation tests (air pressure 
about 10-4 mm of mercury) and four conduction tests were made. 
The conduction tests were in or near the rarefied-gas region, the 
upper limit of which is about 0.01 mm of mercury for the appara- 
tus. Each test was continued until the temperature change was 
within 0.1 deg F per hr. The temperature of the inner cylinder 
ranged from 109.5 to 158 F, and the water temperature, assumed 
to be the same as that of the reference surface, from 71.6 to 79 F. 
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Certain Aspects of High-Pressure 


Centrifugal Pumping Cycles 


By I. J. KARASSIK,! HARRISON, N. J. 


The author discusses the hydraulic and mechanical 
limitations imposed upon centrifugal boiler feed pumps by 
certain features of the modern regenerative feedwater 
heating cycles, as well as the influence of the practical 
limitations of pump design on the cycles themselves. 


posed of a boiler feed pump, a boiler, and a heat engine 

such as a turbine. All other equipment represents re- 
finements which supplement the functions of these three elements 
and which are directed at improving the over-all efficiency of the 
power plant. 

As the size and number of steam-generating stations grew, the 
desire for improvements in operating economy dictated many 
refinements in the steam cycle. These refinements lay in two 
separate directions. The first of these refinements involved a 
steady increase in operating pressures and temperatures, in order 
to provide for a greater heat drop between boiler and condenser. 

Operating pressures grew until at the present time 1250-psi 
boilers are quite common, and a few plants are operating at even 
higher pressures. Temperatures rose from a conservative 600 
and 750 F, until today 900 to 950 F is not considered radical prac- 
tice. Wartime experience with superchargers and gas turbines 
promises to lead to even higher steam temperatures. The direct 
effect of this rise in steam-plant operating conditions will be dis- 
cussed later. 

The second refinement was directed toward greater utilization 
of heat through increased feed heating. The problem presented 
by the introduction of multiple heaters centered in the choice be- 
tween direct-contact heaters and closed heat exchangers. From 
the thermodynamic point of view, the former have certain definite 
advantages. It must be considered, however, that a separate 
boiler feed pump is required after each direct-contact heater, as 
shown in Fig. 1, while the use of a group of closed heaters permits 
the application of a single boiler feed pump discharging through 
these heaters and into the boiler, Fig. 2. 


| Daeeeters to its essential elements, a power plant is com- 


Moutrrete Direct-Contact Haters 


The first of the problems introduced by a feed cycle with 
multiple direct-contact heaters is the requirement of separate 
feed pumps for each of the heaters. This of course may be 
partially circumvented by the use of a multistage boiler feed 
pump with multiple inlets and discharge connections, each stage 
or group of stages serving an individual heater, and each such 
group operating at a successively elevated temperature. 

The most difficult problem, however, resides in the control of 
the boiler feed pump or pumps... A centrifuga] pump will always 
operate at those conditions of head and capacity which corre- 
spond to the intersection between its head-capacity characteristic 
at its operating speed and the system-head curve. Thus in 


1 Application Engineer, Worthington Pump & Machinery Corpora- 
tion. Mem. A.S.M.E. 

Contributed by the Hydraulic Division and presented at the Spring 
Meeting, Birmingham, Ala., April 3-5, 1944, of Tue AMERICAN 
Society or MecHANICAL ENGINEERS. 

Nore: Statements and opinions advanced in papers are to be 
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of the Society. 
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order to alter the operating conditions, it becomes necessary to 
alter the shape either of the head characteristic by varying the 
operating speed, or of the system-head curve by varying one of 
its components. In a boiler-feed cycle this variation is provided 
by the feedwater-level regulator, as illustrated in Fig. 3, where 
curves (a) to (d), inclusive, correspond to various positions of the 
feedwater-regulator valve. Thus, depending on the boiler 
load and consequently on the position of the feedwater-regulator 
throttling valve, the pump-discharge capacity will correspond 
to the intersection of the H-Q curve with curves (a) to (d), or with 
any intermediate curve. 

In the case of a feed cycle with multiple direct-contact heaters, 
as in Fig. 1, it becomes necessary to interpose a feedwater regu- 
lator in the discharge of each individual boiler feed pump. The 
controlling impulse would be the boiler-drum level in the case of 
the last discharge and the level in the storage space of the direct- 
contact heaters in the intermediate units. This multiplication 
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Fig. 3 OPERATION OF FEEDWATER REGULATOR 


of controls introduces additional expense as well as an increase in 
the source of potential maintenance troubles. 

Of course, any centrifugal pump is, to a certain extent, self- 
regulating as regards delivery through the effect of suction condi- 
tions on the pump head-capacity curve. A typical example of a 
self-regulating installation is that of a hot well or condensate- 
removal pump. In this particular case, the delivery of the pump 
is generally controlled by the submergence over the pump center 
line, this submergence representing the net positive suction head 
available over and above the vapor pressure. If the inflow to the 
condenser hot well is reduced, the pump will withdraw more water 
than is returned to the condenser with a resultant change in the 
level of the hot well, a reduction in the available net positive 
suction head, and a consequent reduction in the delivery of the 
condensate pump. This operation is illustrated in Fig. 4. 

In theory the same principle of operation could be applied to a 
feed cycle with multiple direct-contact. heaters. The delivery of 
the last boiler feed pump would be controlled by a feedwater 
regulator while the intermediate units would be subjected solely 
to the control of the available net positive suction head, that is, 
of the submergence from the storage space of the direct-contact 
heater. As the boiler load would be reduced, each pump in turn 
would continue for a short period of time to handle more feed- 
water than was being returned into the direct-contact heater at 
the head of the system, with the result that levels in the storage 
space would fall and reduce the pump capacity through the re- 
duction of the net positive suction head. ’ 

In practice, however, this would prove to be a difficult under- 
taking. It must be considered, in the first place, that in order to 
obtain satisfactory capacity variations through level control, 
the available net positive suction head itself must undergo vari- 
ations of as much as 80 per cent or even more of its maximum 
value. This is entirely feasible in the case of a condenser hot 
well, where the available net positive suction head may vary 
between 4 ft and 1 ft. Condensate pumps are properly propor- 
tioned and designed to avoid excessive inlet and rotative speeds 
and are provided with special condensate-type impellers. The 
application of such designs to high-head boiler feed pumps 
would not be practical, as the latter operate at high speeds and 
require net positive suction heads of 15 to 20 ft or even more. 
Even if it were possible to provide suction-level control to the 
operation of boiler feed pumps; the storage space of the direct- 
contact heaters would have to extend to within 2 or 3 ft of the 
pump center line, increasing the cost of the heaters quite appre- 
ciably. 

There is, finally, one more factor which presents difficulties in 
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the application of a feed cycle with,multiple direct-contact heat- 
ers. Operation of boiler feed pumps at extremely reduced capaci- 
ties is conducive to a variety of undesirable effects due to the 
temperature rise experienced within the boiler feed pump at light 
flows. In the normal case, the boiler feed pump is provided with 
a by-pass in its discharge, leading to some region of lower pres- 
sure in the feed cycle where the excess heat in the by-passed feed- 
water may be dissipated without harmful effects. This by-pass 
is operated either manually or automatically whenever the flow 
through the pump falls to some predetermined minimum or when- 
ever such a reduction in capacity is expected. A boiler-feed 
cycle such as illustrated in Fig. 1 would require the installation 
of an individual by-pass beyond each separate boiler feed pump, 
and each by-pass would require its own piping and its own 
individual control. An automatic-control system would be 
quite feasible but, again, would add to the initial investment and 
probable maintenance expense. f 
We have not touched upon the problem of heater construction 
because this analysis is concerned strictly with the effect of 
various feed cycles on the centrifugal pumping equipment in-- 
volved. There is little doubt, however, that certain problems 
must arise in the design and operation of direct-contact heaters 
at the elevated pressures and temperatures imposed by the 
application of a feed cycle with multiple direct-contact heaters. 


MuttieLe Ciosep HEatTErs 


It cannot be said that the application of a feed cycle with 
multiple closed heaters is entirely devoid of problems. As a 
matter of fact, the average power-plant cycle is based on a 
modified system; one direct-contact heater fulfills the function 
of feedwater deaeration, while several additional heaters of the 
closed type are located upstream as well as downstream of the di- 
rect-contact heater and its boiler feed pump. 

The condensate drains from the closed heaters must be returned 
to the cycle in order to avoid a waste of treated feedwater and of 
its heat and pressure energy. Two basic methods are available 
for the reintroduction of the heater drains into the feed cycle: 

1 Pumping the drains into the feedwater cycle at some higher 
pressure point, as in Fig. 5. 

2 Flashing the drains into the steam space of a lower-pressure 
portion of the feed cycle, as in Fig. 6. (When this is done suc- 
cessively with several closed heaters, the process is referred to as - 
“cascading.’’) ; 

Thermal-balance calculations indicate that the first method 
shows a gain in over-all plant economy over the straight cascading 
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system. However, the economy to be gained by the inclusion of 
individual heater drain pumps at each closed heater is hardly 
justifiable in view of the added initial expense and greater mainte- 
nance required. As a matter of fact, pumping from the closed 
heaters on the discharge side of the boiler feed pump in a high- 
pressure installation is quite impractical. The capacity of the 
heater drain pumps is but a fraction of the design capacity of the 
feed pump itself, while the discharge pressures of the two are of 


an equivalent order or magnitude. As a result, it may become 
impossible to apply a centrifugal pump for the service strictly 
because of hydraulic-design considerations. 

A compromise between these two basic methods is illustrated in 
Fig. 7. It consists of cascading the heater drains of all closed 
heaters located on the discharge side of the boiler feed pump 
back into the direct-contact heater and of pumping the drains of 
the low-pressure heater into the discharge from the condensate 
pump. The over-all losses in plant economy brought about by 
this modified cascading system, when compared to the straight 
pumping system, are negligible. 


Errects or RIsE IN OPERATING PRESSURES 


It has been mentioned that the rise in power-plant operating 
pressures had a very definite effect on centrifugal boiler-feed- 
pump construction and operation. The pump of some 20 or 30 
years ago seldom operated at speeds in excess of 1800 rpm. 
While boiler pressures hovered around 200 to 250 psi, there was 
relatively little advantage to be gained from the introduction of 
higher operating speeds for any but the smallest-capacity boiler 
feed pumps. However, as soon as boiler pressures started on 
their upward climb, the lower operating speed of 1800 rpm be- 
came no longer practical and the modern practice of 3600 rpm 
for boiler feed pumps was launched. 

It is interesting to note that high-pressure installations of the 
central stations have distinctly showed the way for the smaller 
units which started on a parallel growth toward 300, 600, and 
now even 1200 psi. What 3600 rpm did for the medium- and 
large-capacity boiler feed pump, it could no longer do for the 
smaller-capacity unit, until today there is a definite growing 
trend toward operating speeds in excess of 3600 rpm. 

This might serve, incidentally, as a logical occasion to discuss 
the oft-quoted technical superiority of our European counterpart 
in the field of centrifugal-pump design. The more frequent use 
of high operating speeds in Europe stems not from the greater 
talent and ingenuity of their engineers but rather from a series of 
economic and geographical circumstances. 

The first and probably least important of these circumstances 
is the fact that European electric-distribution practice is based 
on 50-cycle current, permitting a maximum synchronous speed 
of only 3000 rpm, compared to a possible 3600 rpm in the United 
States. This lower limit handicapped European designers suffi- 
ciently to make it desirable to exceed the 3000-rpm speed, either 
through greater use of turbine drive or through the application of 
stepup gears. Once the natural speed limit of motor-driven 
pumps was exceeded and abandoned, there was little reason to 
carry the increase to only 3600 rpm, and higher-speed pumps 
were widely introduced. 

The more important reason, however, is the fact that, while 
raw-material costs are generally higher in Europe, labor costs 
were always abysmally low. The European designer has 
through necessity learned to spare materials even at the cost of 
increasing labor man-hours, both in production and in mainte- 
nance. Where an American engineer will use a safety factor of 5, 
European training dictates the choice of lighter structures, close 
figuring, and safety factors of 2 or even less. 

Even after the advent of 3600 rpm as standard practice in this 
country, the rise in boiler pressures kept pace. Soon, discharge 
pressures in excess of 1200 psi became a common requirement and 
even the use of cast steel for the pump casings proved an in- 
sufficient barrier to internal leakage of the then common axially 
split casing boiler feed pump. The next step was the introduc- 
tion of radially split casing designs, culminating in the well- 
known forged-steel-barrel double-casing construction of today for 
the higher-pressure range. It is interesting to note, incidentally, 
that the dividing line of approximately 1250 psi between axially 
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and radially split casings appears due for a slight upward revision, 
probably to 1350 or 1400 psi. 

Damage to pump stationary parts was found to be the most 
troublesome item standing in the way of such upward increases 
in operating pressures, and it is being overcome now by the use 
of stainless-steel alloys of similar grade to those already in use 
for the rotating pump parts, such as impellers, impeller rings, 
and shaft sleeves. Thus the old dividing line will be greatly 
affected by a more general adaptation of erosion-resisting ma- 
terial in all the parts of the pump. Part of this revision may be 
said to derive from greater accumulated experience with axially 
split casings operating near the old dividing line; part may be 
ascribed to the availability of better materials than heretofore 
used, 

Once the general principles of design for this higher pressure 
range had become more or less stabilized, there still remained the 
question of the limiting stage pressures. From a strictly hy- 
draulic point of view, there was much to be gained from main- 
taining low- or medium-stage pressures, at the cost of multiplying 
the required number of stages. The lower heads per stage kept 
disk-horsepower losses to a minimum and avoided possible reduc- 
tion of pump efficiencies. It also became necessary to reconsider 
the effects of the higher velocities brought by higher stage pres- 
sures on the erosion problems in the pump casings. The use of 
extremely low heads per stage, however, leads to a greatly in- 
creased number of stages, greater shaft spans and shaft deflec- 
tions, and greater stresses on the pump casings. The increased 
number of running joints in a pump designed with too great 
a number of stages added seriously to its disadvantages. A satis- 
fying compromise had to be reached, taking into account available 
materials, the optimum number of stages, and the required oper- 
ating speed, the main consideration being to obtain a pump de- 
sign which can be maintained in service with the least expenditure 
of cost, effort, and expert knowledge. 

As a result of the availability of better materials, it has been 
the practice lately to increase stage pressures and keep the num- 
ber of stages to an absolute minimum. Six stages are considered 
today the best possible combination for most pressures in excess 
of 1000 psi, although in the case of relatively low-capacity 
pumps it becomes necessary once in a while to provide seven or 
even eight stages. Assuming a net pressure of 1800 psi, a six- 
stage pump would develop 300 psi per stage (or about 725 ft). 
Once thought excessive and even impossible, 725 ft per stage 
can be considered as a very safe and reliable figure today. 


Hicu Sucrion TEMPERATURES 


Along with rising discharge pressures, boiler feed pump de- 
signers have had to contend with a trend toward higher suction 
temperatures and higher suction pressures. This trend can be 
traced to two separate power-plant developments. In certain 
cases, the desire to improve the station heat balance dictates the 
choice of a higher pressure stage in the main turbine for feed heat- 
ing at the direct-contact heater. As a result, instead of 3 to 15 
psig, direct-contact heaters are being more frequently designed 
for 65 psia, and even higher. In other cases, the pumping tem- 
peratures are increased through the application of two boiler 
feed pumps operating in series with one or more closed heaters 
being located between the primary and main feed pumps. This 
frequently occurs in topping installations, where the pumps 
which fed the low-pressure boilers are retained and serve as 
booster pumps to the higher-pressure addition. 

When the direct-contact heater is fed with steam from a rela- 
tively high pressure stage of the main turbine, special attention 
must be given to the boiler feed pump suction conditions if the 
heater pressure is not held at a constant value regardless of bleed- 
pressure fluctuations, because there exists a danger of flashing in 
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the pump suction whenever a severe and sudden drop is experi- 
enced in the station load. 

As the generator dumps its electric load, the first tendency of 
the main turbine is to overspeed, because it is receiving an excess 
of steam over the requirements of the lighter load. However, 
since turbogenerators are so arranged in their controls as to pre- 
vent asynchronous operation, the turbine governor immediately 
acts to throttle the steam supply to the turbine and holds the 
turbine speed at its rated value. 

The consequence of this reduction in pressure at the steam inlet 
is a proportionate reduction at all the successive stages of the 
turbine and of course at the bleed stage supplying the direct- 
contact heater. In the meantime, the water in the storage space 
of the heater, which was at a temperature corresponding to the 
pressure in the heater, is suddenly subjected to a severe decrease 
in pressure without a corresponding decrease in temperature. 
The net result of this situation is severe flashing which occurs in - 
the heater as well as in the feedwater in the suction piping to the 
boiler:feed pump. 

The time element involved in again bringing the pressure and 
temperature within the heater to a correlated value will of course 
vary with each individual case. The variation in the time com- 
ponent is due to the great variation in the size and shape of the 
storage space, since these two elements are important factors in 
the rate of flashing at the free surface of the feedwater. 

It must also be considered that the downward velocity of the 
feedwater in the suction piping is too high to permit the vapor 
bubbles generated by flashing, to rise up to the heater and thus 
to clear the boiler feed pump suction of the flashing steam. A 
safety factor must be provided to take care of the time lag which 
exists between a certain reduction of pressure in the heater and 
the ultimate reduction of temperature at the pump suction. Itis 
therefore necessary either to provide an excessive factor of 
safety in the installation of the boiler feed pump through the 
addition of an empirical percentage to the calculated net positive 
suction head or else some means whereby any appreciable pres- 
sure reduction in the heater would be accompanied immediately 
by a temperature reduction at the pump suction. 

The latter solution is obviously predicated upon the injection 
of colder water into the pump suction. While the design of a 
control for such injection is not an insurmountable problem, the 
practicability of this solution is governed by circumstances _ 
peculiar to each individual installation. The cooler injection 
water must be taken from the downstream side of the cycle, that 
is, either from the discharge of the condensate pump or from the 
drains of a low-pressure closed heater ahead of the direct-contact 
heater. To introduce nondeaerated condensate into the boiler 
feed pump suction is not recommended, especially in the case of 
high-pressure boilers where the possible effects of contamination 
are severe. It must be considered, however, that a high-pressure 
installation, including a direct-contact heater operating at high 
temperature, will probably also include a low-pressure closed 
heater between the condensate discharge and the direct-contact 
heater itself. It must be further considered that the drains from 
this heater, generally pumped by a separate pump into the direct- 
contact heater, are to all intents and purposes pure deaerated 
condensate, and that their injection into the feedwater suction 
stream will cause no harmful effects. 

A logical precaution, therefore, might consist in the installation 
of a by-pass line from the low-pressure heater drain pump dis- 
charge into the boiler feed pump suction. This by-pass line 
would be properly valved and controlled by some means such as 
the difference between the suction pressure and the vapor pressure 
at the suction, any decrease of this difference below a predeter- 
mined minimum causing the by-pass valve to operate and inject - 
colder water into the suction of the feed pump. Two charts 
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have been prepared to illustrate the effect of colder water injec- 
tion on the available net positive suction head. Fig. 8 shows 
the temperature depression resulting from the injection, plotted 
against the difference in temperature between the feedwater and 
the injection stream and given for varying ratios of injection. 
Fig. 9 illustrates the effect of temperature depression on the 
available net positive suction head at various feedwater tem- 
peratures. 
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Thus, for instance, if it is desired to provide 20 ft surplus net 
positive suction head at the boiler feed pump suction and if the 
feedwater temperature is 325 F, the required temperature depres- 
sion is 6 deg F. Assuming that the temperature of the heater 
drains is 227 F, allowing a 98-deg differertce between feedwater 
and injection water, the by-pass line must be proportioned to pass 
0.062 of the maximum boiler feed pump capacity, or 6.2 per cent. 

Another serious problem introduced by high feedwater tem- 
peratures, especially when coupled to high suction pressures, re- 
sides in the operation of the stuffing boxes. The design of the 
pump interior parts can be accommodated to these conditions 
without undue difficulties, now that high-grade alloy stainless 
steels are readily available. The practice of shrinking impellers 
on the shaft in high-pressure applications has further helped 
solve some of the earlier problems. Finally, the design of the 
pump casing itself, through the application of materials such as 
forged steel and in some cases even stainless alloys, as well as 
through certain mechanical improvements, has likewise elimi- 
nated some of the earlier temperature limitations. 

There remains, however, the problem of providing a reliable 


and economical stuffing-box design. Certain objectives may be 
said to have been attained even in this field. For instance, within 
the last few years a more thorough understanding of the inter- 
relation of stuffing-box pressures, rubbing speeds, and leakage 
temperatures has led to an improved style of water-cooled stuffing- 
box design. Such stuffing boxes have been built suitable to with- 
stand feedwater temperatures up to 400 F and stuffing-box pres- 
sures up to 500 psi, without recourse to pressure-reducing laby- 
rinths. The construction of such a stuffing box is illustrated in 
Fig. 10. It will be noted that one of the outstanding features of 
this design is that the stuffing-box leakage is precooled before it 
reaches the packing proper. Operating records in excess of 6 
months have been obtained with this construction without re- 
newing the packing. 

While this combination of pressures and temperatures does not 
necessarily represent the ultimate limit, further increases will 
probably be reached more slowly and after much careful labora- 
tory and field experimentation. Furthermore, it may be noted 
that, at the present time and with very few exceptions, this 
range of pressures and temperatures represents the maximum 
conditions which might have to be imposed on boiler feed pump 
stuffing boxes. 

Whenever boiler feed cycle requirements are such that these 
limits have to be exceeded, the only solution resides in the 
application of pressure-reducing devices ahead of the stuffing box 
proper. Since the design of such devices has now reached a high 
degree of proficiency, the only remaining drawback to their 
application is the necessity of bleeding a part of the effective. 
capacity back to a lower pressure level and the attendant reduc- 
tion in efficiency. 


PRESSURE AND TEMPERATURE CONTROL IN CLOSED FEED CYCLES 


There is a certain interesting although little known phase in 
the control problem of the so-called closed feed cycle, where 
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several boiler feed pumps are operated in series with closed 
heaters interposed between successive pumps, and where one or 
more of these pumps is operated at variable speed, as illustrated 
in Fig. 11. 

Normally, the variable-speed operation is restricted to the 
low-pressure unit of the combination, the controlling mechanism 
being either an excess-pressure or a constant-pressure regulator. 
This regulator is used to vary the speed of the centrifugal pump in 
question, which is driven by a steam turbine, a wound-rotor 
motor, or even through a variable-speed coupling. 

This control is obviously unaffected by the relationship between 
the temperature and the pressure at the suction of any one of the 
pumps so that conditions might arise where the speed regulation 
on one pump lowers the suction pressure of the following pump 
dangerously close to or even below the vapor pressure of the 
feedwater at the suction. 

While at full load the temperature of the feedwater leaving the 
individual pumps is for practical purposes equivalent to the tem- 
perature at the pump suction, this does not hold true at light 
loads, and for that reason the illustration indicates this by using 
different subscripts for the temperatures at the suction and dis- 
charge of each pump. It is obvious that the margin between the 
pressure at the suction of the secondary feed pump (P4) and the 
vapor pressure corresponding to the temperature 7’, must never 
be allowed to fall below a certain predetermined value which de- 
pends on the pump design. 

The selection of the individual pump-discharge pressures is a 
relatively complex problem which is affected by the amount of 
required pressure regulation, that is, by the system head curve 
of the boiler feed cycle, by the head-capacity characteristics of 
the individual pumps, the maximum temperatures of the feed- 
water at the outlet of the various heaters, the permissible speed 
variation of the variable-speed drivers, and finally by the desired 
maximum hydraulic pressure to be imposed on the inter- 
mediate closed heaters. 

Fig. 12 illustrates qualitatively the characteristics of a typical 
closed system such as in Fig. 11. For purposes of simplification, 
frictional losses between the pumps have been neglected. It will 
be seen that the pressure at the discharge of the secondary pump 
considerably exceeds the required system head pressure at light 
loads if all the pumps are operated at constant speed, and that 
even if the primary feed pump is by-passed altogether, the 
secondary feed pump can develop sufficient pressure from shut- 
off to capacity ‘“‘A.” In other words, the speed variation re- 
quired of the driver of the primary pump ranges all the way from 
0 to 100 per cent. This, obviously, is impractical for a great 
many reasons, especially if the pump is driven by a wound-rotor 
motor where reduction of speed below 50 per cent of design speed 
is seldom practiced. Asa result, the operation of such a cycle at 
variable speed generally takes place along the lines illustrated in 
Fig. 13. Whatever excess pressure remains is throttled by the 
feedwater regulator. 

The important feature of a proper selection of primary feed 
pump conditions and of a suitable variable-speed control is that 
means must be made available to prevent the discharge pressure 
of the primary feed pump (P; in Fig. 11), from falling below a 
certain minimum determined by the relationship between P, (P3 
less losses through the closed heater), and the vapor pressure 
corresponding to the temperature 7. Should certain operating 
conditions arise which might lead to a dangerous reduction of the 
difference between P, and the vapor pressure at 7's, some control 
must be instituted, actuated by this difference and arranged to 
set a definite limit to the speed reduction called for by the excess- 
pressure regulator, 

Fig. 13 indicates the introduction of such a speed-reduction 
limiting control, whereby the discharge pressure of the primary 
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feed pump is not permitted to fall below a predetermined mini- 
mum. For the sake vf simplification, it has been assumed that 
the value of 7, is constant at all loads, and therefore that the 
minimum permissible value of P; is likewise constant. Actually, 
the closed heaters are fed with steam bled from some given stage 
of the main turbine and the pressure and temperature of this 
steam are reduced as the load is reduced. Thus in an actual 
installation, the process of checking the minimum permissible 
primary feed pump discharge pressures may be a little more 
tedious than indicated here. 


CoNCLUSIONS 


The foregoing discussion of boiler feed cycles has been pre- 
sented not as an attempt to catalogue in an exhaustive manner 
all the possible variations and ramifications which can be found 
in these cycles, but rather with the desire to illustrate the close 
relationship between cycle and pump design. There is no ques- 
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tion that the degree of integration of the pump into the pumping 
cycle itself is a direct criterion of the ultimate success of the 
entire installation. It is toward such complete integration that 
the pump designer, the application engineer, and the technical 
operating personnel have had to apply their efforts, and it. is 
through this integration that modern centrifugal-pump designs 
have achieved their relative success. 

It is hoped that the present analysis will give rise to more 
studies of all the factors connected with high-pressure cycles, and 
thus the seed will be sown for further improvements based upon 
mutual understanding and co-operation between designer and 
operating personnel. The problems which high-pressure pump- 
ing cycles have presented have not grown easier with time. But 
it may be said, with a certain pardonable pride, that the centrifu- 
gal-pump designer has so far fully met the challenge. 


Discussion 


R. P. Moors.? The author describes a clever means of pre- 
venting loss of the suction pressure on boiler feed pumps drawing 
water from direct-contact deaerating heaters when the load on the 
turbine is reduced. This scheme includes some automatic valy- 
ing. Experience has proved, in certain cases at least, that it is 
possible to accomplish the same result without any automatic 
valves or other moving equipment. Proper design of the de- 
aerating heater, as is described in an article by the writer,® ac- 
complishes the desired result. The scheme consists of placing 
baffles in the heater so as to guide the freshly deaerated water 
directly to the boiler feed pump suction and to allow any flash 
steam which forms in the water in the heater to rise to the steam 
space without interference with this flow. 

Fig. 14 of this discussion shows the baffles and vents added to 
the bleeder deaerating heater of standard design in use in Huntley 
Station No. 2 of the Buffalo Niagara Electric Corporation. The 
results accomplished by this addition are shown in curve form in 


2 Mechanical Engineering Department, Buffalo Niagara Electric 
Corp., Buffalo, N.Y. Mem. A.S.M.E, 
3 Electrical World, October 14, 1944. 
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Figs. 16 and 17, taken from the article referred to, compared to 
the conditions prior to the change, Fig. 15. In Fig. 16, the rate 
of pressure drop was such that deaerator and stage pressures re- 
mained the same. In Fig. 17 the stage pressure dropped so 
rapidly that the bleeder nonreturn valve closed, and the deaerator 
pressure could not follow. Not only has loss of pump suction 
been eliminated when the turbine load drops, even if instantane- 
ously, but the water continues to pass through the deaerator and 
is deaerated under all conditions. 


R. M. Sraco.! 


In connection with higher boiler feed pump 


4 Efficiency: Engineer, New Orleans Public Service, Inc., New 
Orleans, La. Mem. A.S.M.E. 
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suction temperatures, the author states that in the event of a 
severe and sudden drop in station load, turbine-stage pressures 
will quickly drop resulting in severe flashing in the direct-contact 
* heater serving as boiler feed storage tank and also in suction pip- 
ing to the boiler feed pump. While drop in pressure of bleed 
steam to the boiler feed heater undoubtedly contributes to the 
condition of severe flashing, it is not, in the writer’s opinion, the 
only factor involved. 

To illustrate, let us follow the train of events resulting from a 
severe and sudden drop in station load: 


1 Due to sudden decrease in load on the electric generator, 
the turbine tends to overspeed, resulting in a rapid throttling of 
control valves by the governor. 

2 Asa result of control valves throttling, stage pressures drop 
rapidly and throttle pressure rises rapidly. 

3 Sudden rise in throttle pressure may result in one of two 
effects, in so far as the boiler is concerned, depending on the 
severity of the load drop, that is, depending on whether it is 
within or without the capability of combustion control. 


If the drop is too severe for combustion control to handle, one 
or more safety valves will lift, causing an instantaneous rise in 
boiler-drum water level. The feedwater regulator might close 
completely, and if automatic recirculation is provided the boiler 
feed pump will be subjected to minimum flow as determined by 
its recirculating orifice capacity. If manual recirculation control 
is provided, under these conditions, the pump will probably be 
subjected to zero flow, resulting in severe flashing and probable 
damage to its internal parts. 

On the other hand, if the load drop is within the capability of 
combustion control to prevent popping of safety valves, another 
train of events, detrimental to boiler feed pump operation, will 
be started. Rise in throttle steam pressure will result in an 
initial rise in boiler-drum pressure causing a drop in drum water 
level, followed by a more severe drop in water level as the firing 
rate is rapidly decreased. If the feedwater regulator is of the 
single-element type, it will open rapidly, placing a sudden demand 
on the boiler feed pump and also on the boiler-feed tank. 
If the regulator is of the multiple-element type, the water-level 
effect will probably override the effect of diminished steam flow 
to result in a similar though less pronounced increase in feeding 
rate. This will result in a rapid decrease in boiler feed pump 
npsh over and above the decrease imposed by drop in turbine 
stage pressure. 

Further aggravation will also be imposed in some installations 
by quenching of vapor in the boiler feed tank by make-up water 
supplied from an auxiliary storage or surge tank. 

If return of turbine loading is rapid, as in the case of a short- 
duration electrical interruption, similar hardships will be placed 
on the boiler feed pump. Rapid increase in boiler firing rate 
will result in almost instantaneous swelling of water within the 
boiler, often causing the feed regulator to close completely. 
Again the boiler feed pump will be subjected to minimum flow 
conditions followed by a gradual return to normal flow. 

The foregoing is certainly not intended as a criticism of the 
author’s analysis of causes of erratic boiler feed pump suction 
conditions, but rather is intended as an amplification of this 
analysis in line with a personal opinion that correction of such 
conditions should be, in the majority of cases, thrown back in 
the lap of engineers responsible for the heater, piping, and pump 
layout. Obviously, sufficient submergence of boiler feed pumps 
to provide ample static suction head in excess of any reductions 
imposed by normal surges would remedy the situation without 
recourse to injection of colder water into the boiler feed pump 
suction. For the majority of cases where space or operating 
pressure limitations will not allow provision of sufficient static 
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suction head, the injection of colder water offers definite possi- 
bilities. ; 

The use of low-pressure heater condensate as injection water 
presents another problem in connection with many existing sys- 
tems. Ina typical installation, for example, at full turbine load, 
low-pressure heater condensate represents 8 per cent of throttle 
steam flow (or boiler feed pump flow). However, at one half 
load, low-pressure heater condensate is only 5 per cent of throttle 
flow. This may or may not provide sufficient injection water to 
insure adequate temperature depression at the boiler feed pump 
suction. 

In line with thesg possible restrictions, it might be well to re- 
examine the possibility of using condensate direct from main 
condensate pumps. It has been our experience, based on 
numerous tests on two modern condensers designed for de- 
aeration of hot-well water, that condensate leaving condensate 
pumps is oxygen-free within limits of detection, by the modified 
Winkler test, and consequently, from this standpoint, is suitable 
for injection directly into the boiler. 

Papers of this type, together with the discussions which they 
evoke, are of tangible benefit to the engineer faced with similar 
problems. Encouragement should be given to preparation of more 
papers in which the equipment designer freely discusses the 
limitations as well as the capabilities of the apparatus. 


A. J. Srepanorr.’ The writer will call attention to a few 
problems in high-pressure boiler feed practice not mentioned in 
the paper, which are not caused directly by the high pressure, 
high temperatures, or pump high operating speed. One of them 
is the wear of internal parts resulting from corrosion and erosion 
commonly referred to as “cutting.” The pH value of the 
pumped water does not indicate the corrosiveness of water when 
applied to pure water. A definite trend'to purer feedwater in 
central-station application will demand wider use of stainless 
steels for pump parts than has been the case in the past. 

Mechanical problems connected with high-pressure high- 
temperature boiler feed pumps were not different from those 
encountered with the hot-oil charging pumps. A number of 
these were in service some 10 to 15 years ago, pumping oils at 800 
F, maximum pressures from 1000 to 1500 psi (2300 psi was used 
at one refinery) and 4000 rpm pump speed. Even with the most 
corrosive (sour) oils, pump-casing cutting was practically un- | 
known with hot-oil pumps. 

Pressure “‘surges,’’ accompanied by the capacity, speed, and 
power fluctuations, presented another problem in boiler feed 
pump operation not previously encountered in hot-oil or any 
other pumping service. Such surges are of the same nature as 
those observed in hydroelectric power plants frankly and com- 
pletely discussed in a paper by Rheingans.* Both systems have a 
mass of water suspended between the free surfaces at their ends. 
This is free to oscillate if some impulses appear at a frequency 
synchronous with that of the natural frequency of the system. 
As a result of such fluctuations, in the case of boiler feed pumps, 
the piping may require additional bracing, and some pumps can- 
not be operated at low capacities. R. Dziallas? has studied the 
surges and was able to reproduce those in a laboratory setup. 

It should be noted that condensate-pump cycle presents 
another swinging system with two free water surfaces at the ends, 
but no surges are observed on this system as long as the con- 
densate pumps operate beyond the breakdown point on the head- 


5 Engineer, Cameron Engineering Department, Ingersoll-Rand 
Company, Phillipsburg, N. J. Mem. A.S.M.E. 

6 ‘*Power Swings in Hydroelectric Power Plants,” by W. J. 
Rheingans, Trans. A.S.M.E., vol. 62, 1940, pp. 171-184. ; 

7 “Schwingungen bei Kreiselpumpen mit labiler Kennlinie,” by R. 
Dziallas, Archiv fiir Warmewirtschaft und Dampfkesselwesen, March, 
1941, pp. 63-66. : 
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capacity curve, Fig. 4 of the paper, because cavitation destroys 
the regularity of impulses in the system. Similarly, surges can 
be stopped on any boiler feed pump by throttling the pump dis- 
charge slightly (sometimes as little as 20 psi), thus destroying the 
regularity of the water mass oscillation. 

The writer does not believe that the author’s criticism of 
European engineering or training is justified. The high-pressure- 
boiler development is only one example of European engineering 
achievement. 

In so far as pump speeds are concerned, the general trend 
toward higher speeds is very definite in this country too. While 
3600 rpm was a top speed only recently, 6000 to 7000 rpm is 
quite common now. — Cost is only secondary to weight and space 
saving, and hydraulic advantages in the advance in the operating 
speeds of boiler feed pumps. 


G. H. Van Henceu.’ The paper gives a good detailed descrip- 
tion of the difficulties encountered when operating centrifugal 
pumps with net positive suction heads close to their guarantee. 
These pitfalls are generally not realized when the layout of a 
plant is made and are difficult to correct afterward because they 
are dependent on the elevations which have been fixed by the 
structure of the building. 

In large power stations where bleeding is used, a change of load 
will change the pressures in all the turbine stages practically pro- 
portionally to the load. Thus in case 20 per cent of the turbine 
load is dropped, the stage pressures will also drop 20 per cent, and 
with them the heater pressures, All the pumps extracting from 
these bleeding heaters will therefore be suddenly robbed of 20 
per cent of their suction pressure. 

The stage temperatures, however, because of the heat storage 
in the condensate present in the heater, will not drop as rapidly. 
It is this phenomenon of the lag in temperature with respect to 
pressure that causes a difficult condition for the heater pumps. 
Because this lag depends on various conditions of the layout 
and the rate change of load, it differs for each individual case. 
The whole matter reduces to the question of how fast the load on 
the turbine drops in comparison to the ensuing rate of tempera- 
ture changes in the heater. 

This comment is intended to emphasize the severity of the 
condition of pump operation caused by the persistence of this lag, 
and to show a ready method of calculation to provide for any con- 
dition that may occur. 

Let us now consider different cases of installation and required 
operation. For all cases, we assume: That no undercooling of 
the condensate leaving the heaters exists in regard to the steam 
entering the heaters or that the pump is running with a net posi- 
tive suction head at the least value required for its operation. 


The following cases are discussed: 


1 Instantaneous load drop. 

2 Very gradual load drop. 

3 Decrease of load throughout a given period. 
(a) With no storage capacity. 
(b) With storage capacity. 


Case 1. Instantaneous Load Drop. As a basis of calculation, let 
us assume that the load change on the turbine is instantaneous 
and that the temperature in the heater does not change, a case 
which corresponds to infinite storage capacity. For safe opera- 
tion, an excess suction head has to be available to provide a 
corresponding pressure drop. To calculate this excess suction 
head, we can readily determine the head equivalent of the 
assumed sudden pressure drop by using the conservative formula, 
1 psi pressure = 2.6 ft head. Thus a sudden drop of 10 psi ona 
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heater would demand an excess net positive suction head of 26 ft 
on the pump to prevent unsafe operation. 

Ezample: A condensing 80,000-kw turbine with 17 stages has 
a Main steam pressure and temperature of 800 psia and 900 F. 
It has four direct-contact heaters supplied with steam bled from 
the 4th, 7th, 10th, and 13th stages at full-load stage pressures of 
250, 100, 40, and 8 psia, respectively. If there is a sudden load 
drop of 20 per cent, the conditions which result are as listed in 
Table 1 of this discussion. 


TABLE 1 CONDITIONS aa Yep FROM 20 PER CENT LOAD 


——Heater stages——. Con- 

4th 7th 10th 13th denser 

1 Full-load pressure, psia.............. 250 100 40 8 0.50 
2 80 per cent load pressure, 0.80 X (1), 

TT ale ARO GaNncoraee op Omeaber 200 80 32 6.4 0.47 

3 Sat.temp corr. to (1), deg F.......... 401 327 267 183 80 

4 Sat. temp corr. to (2),degF.......... 382 312 254 173 28 

5 cemp. diff. due to load, (3-4), deg F... 19° 915 13 10 2 
6 Head increment needed on _ heater 

pumps [2.6 (1-2)], ft water........ 130 52 21 4 0.08 


Table 1 shows that with an instantaneous drop of load of 20 
per cent and with the greatest possible temperature lag in the 
heaters, we would have to provide 130, 52, 21, 4 ft additional head 
on the pumps if this instantaneous drop in load is to be provided 
for. 

For a hot-well pump on a condenser with automatic level con- 
trol, however, the table shows that the sudden 20 per cent drop of 
turbine load demands only 0.08 ft excess suction head. There- 
fore, for the condenser hot-well pump, sudden drops of turbine 
load can never be a serious factor to be considered in design. 

For the heaters this is somewhat different. It must be noted 
here that for the heaters the additional net positive suction head 
required is proportional to the absolute heater pressure. The 
higher the heater pressure, the more probability there is for losing 
the suction of the heater pump and the more margin in suction 
head we should have on the pump. 

-Case 2. Gradual Load Drop. At the other extreme is the 
case of dropping the load over a very long period which will 
cause the heater condensate temperature to follow the bled- 
steam temperature with no lag. 

The suction head which must be provided on the pumps in 
excess of the manufacturer’s guarantee would then be zero feet 
for this case. 

The actual provision of suction head to be made must lie be- 
tween the two extreme Cases 1 and 2. 

Case 3(a). Decrease in Load Throughout a Given Period of 
Time—No Storage Capacity. Now assume that the load drops 
not instantaneously, but 20 per cent in 5 min; furthermore, that 
the condensate level is not in the heater but in the suction pipe; 
a case with practically no storage capacity. Then for the 4th- 
stage heater in the example, the condensate should increase its 
pressure by flowing downward at a rate of at least 130 ft in 5 min, 
or 0.43 fps if we want to prevent flashing. This 0.43 fps will thus 
determine the maximum size of our pump suction line. 

If we have an existing installation, then our suction velocity 
should be greater than 0.43 fps for successful operation. 

Case. 3(b). Decrease in Load Throughout a Given Period of 
Time—With Storage Capacity. Assume again 20 per cent load 
drop in 5 min. The case can be considered first as Case 3(a), but 
the hot-well cross section replaces the suction-line area. In 
addition, the colder condensate flow will mix with the hot-well 
storage, reducing the required suction head. This problem has 
to be solved mathematically and is dependent on heater hot- 
well capacity and rate of steam going to the heater. Space 
limitations prevent treating this case. 


All cases so far treated assumed no undercooling or no addi- 


74 TRANSACTIONS OF THE A.S.M.E. 


tional net positive suction head. In direct-contact heaters, 
the undercooling will be practically zero. In closed heaters some 
undercooling always exists. Assume a 5-deg F undercooling for a 
closed heater on the 4th stage of the example given. The net 
positive suction head needed would be reduced by 


5 
= 130 = 34 ft 
19 y 


Also, if regulation of the pump discharge is taking place, it 
may be possible that additional suction head is available at the 
time the load is dropped. 

For safe operation of the heater pump, with time for load drop 
asin Cases 3(a) and (0), the undercooling and the additional suction 
head available during operation reduce the maximum head re- 
quired under Case 1. 

Conclusion. Power stations should be designed not only to 
meet maximum load conditions, but also certain rate changes in 
load. For instance, the requirements for the net positive suction 
head on heater extraction pumps (boiler feed or heater drain 
pumps), due to sudden load drops, are not sufficiently recognized. 
Oversight of this fact often leads to dangerous operation which is 
difficult to correct because of the limitations of the building 
structure. The problem of variability in load during operation 
must therefore be thoroughly investigated during the pre- 
liminary layout of the station. 


AvuTHOR’s CLOSURE 


The author wishes to acknowledge with appreciation the con- 
tributions of the discussers. 

Mr. Moore’s solution of the problems introduced by sudden 
load fluctuations is an extremely valuable contribution to satis- 
' factory heater operation and, the author believes, should be well 
studied by feedwater heater designers. While the possibility 
remains that this solution may not meet all cases, it is gratifying 
to see that it attacks the problem exactly where it arises, that is 
in the feedwater heater. One of the features of this solution 
which remain to be investigated resides in the effect of the baffling 
arrangement described on installations embodying very large 
feedwater storage volumes, where such baffling may lead to ex- 
cessive subcooling of that portion of the stored feedwater which 
is farthest removed from the heater outlet. The author is con- 
fident, however, that even this problem is not unsurmountable. 

The additional effect of boiler controls on the phenomena 
which take place at a boiler feed-pump suction on a sudden drop 
of main turbine load are very well brought out by Mr, Seago. 
While a sudden increase in boiler feedwater demand may or may 
not impose additional difficulties, depending upon the margin 
available over the minimum npsh at maximum pump flow, the 
reverse case—that of complete closure of the feedwater regu- 
lator—introduces very serious possibilities. 

If we were to assume that the sudden load drop occurs when 
operation was taking place between 50 and 100 per cent full load 
and that the by-pass is manually controlled, it is logical to expect 
that all by-passes are closed at the time of the load drop and the 
pumps may be severely damaged. This, incidentally, is quite 
independent of the suction pressure and temperature and. could 
occur even if the pressure in the heater was pegged to a definite 
value. Much as it may be desired to avoid the introduction of 
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additional controls, the sequence of events outlined by Mr. 
Seago shows the advantages and the added protection afforded 
by an automatically controlled by-pass. 

If, as Mr. Seago states, the low pressure heater condensate at 
one-half load corresponds to but 5 per cent of the throttle flow, 
its injection into the feed pump suction may or may not be suf- 
ficient to afford the required protection against flashing. Assum- 
ing the conditions outlined in the paper, that is, 325 F feedwater, 
227 F heater drains, and a required margin of 20 feet—the in- 
jected condensate must be 6.2 per cent of the boiler feed pump 
capacity and in accordance with Mr. Seago’s statement, the 
available condensate is insufficient. On the other hand, assum- 
ing the same temperature conditions, a 5 per cent injection will 
provide about 17 ft surplus npsh (see Figs. 8 and 9) and may be 
sufficient in some cases. Obviously, where it is possible to use 
the discharge from the condensate pumps, the available flow is 
fully ample and because of the greater temperature differential, 
the percentage of injected feedwater appreciably reduced. It 
must be remembered that the injection is an emergency measure 
and that the slight danger of feedwater contamination is justifia- 
ble, especially since the emergency is not a frequent one. 

The problems of corrosion and erosion of internal boiler feed- 
pump parts mentioned by Mr. Stepanoff have been purposely 
omitted from this paper as their importance merits a complete 
and separate study, while their inclusion would have lengthened 
the present paper unnecessarily without exhausting the subject 
matter. It may be of some interest to mention, however, that 
very. serious studies have been carried out on these problems 
within the past twelve months. The preliminary results seem to 
justify the growing trend toward stainless steels mentioned in 
Mr. Stepanoff’s discussion as well as in the author’s paper. It 
is already becoming apparent that stainless steel casings in a 
range from a minimum of 4 per cent chromium content to the 
25 per cent chromium-20 per cent nickel type, will be utilized 
to a much greater extent than heretofore. Coupled to the al- 
ready prevalent use of stainless steel fittings, this trend gives 
definite indications of solving the twin problems of corrosion and 
erosion, even at the higher stage pressures predicted. 

No criticism of European engineering or training was voiced 
or even implied by the author. The comparison of practices was 
brought out to illustrate the possible effect of economic and geo- 
graphical circumstances on engineering designs and no issue can 
be taken with the fact that European labor is cheaper than in the 
U.S. A. and raw materials less available. 

Mr. Van Hengel’s calculations illustrate very forcefully the 
theoretical head increment needed on pumps taking suction from 
direct contract heaters operating at high pressures, both with 
instantaneous and gradualload drop. Fortunately for the power- 
plant designer, the usual case corresponds to Mr. Van Hengel’s 
Case No. 3 and the entire theoretical increment in npsh need 
not be provided for. The determination of the minimum per- 
missible safety margin in npsh is rather complex but can be 
carried out fairly accurately by the feedwater heater designer, 
once all the contributory factors and component data have been 
made available. It is highly important, however, that power 
stations be designed, as Mr. Van Hengel states, not only to meet 
maximum load conditions, but also the expected range of load 
changes. ; 


The Behavior of a Hot-Wire Anemometer 
Subjected to a Periodic Velocity 


By R. C. MARTINELLI! anp R. D. RANDALL,? BERKELEY, CAL. 


A method is presented which allows a rapid graphical 
prediction of the response of the constant-current hot- 
wire anemometer to a periodic velocity. From such a 
prediction the suitability of the hot-wire anemometer for 
any particular application can be readily determined. 
The method is demonstrated by several illustrative ex- 
amples. 


NOMENCLATURE 


The following nomenclature is used in the paper: 


A = surface area of wire, sq ft 
a = temperature coefficient of resistance, deg F~! 
c = constant of integration 


C, = heat capacity of wire, Btu/(Ib) (deg F) 

= base of Napierian logarithms 

= voltage across hot-wire anemometer at any temperature, 

t, volts 

E., = Voltage across hot-wire anemometer when wire is at 
room temperature, To, volts 

E — Eo, voltage across hot wire referred to room-tem- 
perature voltage, volts 

€m = voltage across hot wire, referred to room-temperature 

voltage, corresponding to mean velocity, volts 


é) = voltage across hot wire, referred to room-temperature 
voltage, corresponding to conditions of zero thermal 
lag, volts . 

f = unit thermal conductance between hot wire and air, 
Btu/(hr) (sq ft) (deg F) (actually f includes convection, 
radiation, and conduction to anemometer supports) 

2 = current through hot wire, amperes 

P(@) = function of time defined in text 
P(v) = function of velocity defined in text 

Q = “characteristic’’ of hot wire, defined in text, volts per hr 
or volts per sec 

R = resistance of wire at any temperature, t, ohms 

Re = resistance of wire at room temperature, ohms 

Sa “ = rate of change of voltage with time, volts per hr or 

volts per sec 
t = temperature of hot wire, deg F 


7. = temperature of air around hot wire, deg F 
-W = weight of hot wire, lb 
6 = time, hr 
3.413 = number of Btu per hr in 1 watt 


ll 


INTRODUCTION 


The hot-wire anemometer, as normally operated, consists of a 
small wire (usually of platinum) suspended between supports in a 
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fluid stream and heated by an electric current. Variations of 
the velocity of the fluid flowing past the hot wire produce corre- 
sponding variations in the heat lost from the wire and thus vary 
the wire temperature, which in turn alters its electrical resistance. 
By proper calibration, the wire resistance can be correlated with 
the velocity of the fluid flowing past the wire, and therefore the 
hot wire can be used for the measurement of fluid velocities. 

The application of the hot wire to the measurement of steady 
velocities or to small velocity fluctuations about a large mean has 
been the subject of many studies (1 to 11). Some work has also 
been done on the behavior of the hot wire when subjected to small 
velocity fluctuations about a zero mean, mainly in conjunction 
with the hot-wire ‘“‘microphone”’ (12 to 15). The hot wire has 
also been utilized to measure fluctuating wind velocities, which 
are usually relatively large variations in velocity of rather low 
frequency (16, 17). 

In all these applications the lag of the response of the hot wire 
to velocity variations, due to the heat capacity of the wire, has 
been given considerable attention. However, no complete 
analysis of the problem has been presented which would allow 
the prediction of the hot-wire behavior when subjected to large 
velocity variations of relatively high frequency. It is the purpose 
of this paper to present such an analysis. 


Fie. 1 Hort-Wire ANEMOMETER 


DEscrRIPTION OF IDEAL SYSTEM 


In the following discussion a hot wire of length Z and diameter 
D is considered suspended between two supports, Fig. 1. The 
actual system will be idealized as follows: 


1 The temperature of the wire is considered uniform at all 


3 Numbers in parentheses refer to the Bibliography at the end of 
the paper. 
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times. This is not strictly true because of the cooling of the ends 
of the wire by the supports, and the radial temperature distribu- 
tion existing in the wire during transient operation. 

2 The wire operates at constant current. This requirement 
is easily met in the actual system. ° 

3 All heat losses from the wire are proportional to the differ- 
ence in temperature between the wire and the ambient air. 

4 The heat loss from the wire at any fluid velocity during 
transient operation is exactly equal to the heat loss at the same 
velocity during steady-state operation. 


Basic Equations 


Since the operation of the hot wire depends upon its heat loss 
to the surrounding fluid, the basic equation describing the opera- 
tion of the anemometer is a heat balance on the wire. During 
transient operation, the energy supplied to the wire by the elec- 
tric current appears in two places; part of it is lost to the sur- 
roundings, and the remainder is used to heat the wire from one 
temperature to another. ‘Thus the heat balance becomes 


we, + fAt—7ro) =34137R.......... {1] 


The first term represents the heat required to change the tempera- 
ture of the wire, by an amount dt in a time d0; the second term 
represents the heat lost to the surroundings; and the last term 
the heat supplied by the electric current. All term sare in Btu 
per hour. Equation [1] is the basic equation for the idealized 
hot-wire anemometer. 

In practice it is more convenient to measure the voltage drop 
across the wire than to measure the wire temperature. The tem- 
perature and voltage are simply related, as shown below. 

The resistance of the hot wire R, at any temperature t, may be 
written as 


R= Ro[l + alt—To)].......-...- SORIA 


where FR is the resistance of the wire at the temperature of the 
ambient airt#. The voltage across the wire at constant current 
is 


E = ik = iRa[l + alt —T.)] = Ho[l + alt—72)]).. [3] 


dE 
dt = TS EY sel Eiki [4] 
and 
(t—To) = joke [5] 


Substituting in Equation [1] yields 


dE fA 
do We we Ba) =( 


413 7Ro 
3 a *) E. 6] 


WC, 
For simplicity, a new variable e = H — He may be defined. 
Then Equation [6] becomes 


ae (4 — 3.413 7Rea\ | _ (Caettesto (7) 
do WC, - WC, ares 


The problem of determining the behavior of the hot-wire ane- 
mometer subjected to a periodic (or fluctuating) velocity is then to 
solve Equation [7], obtaining e as a function of time, corre- 
sponding to a given variation of the velocity as a function of time. 
Although the velocity, per se, does not appear in Equation [7], 
variations in velocity cause definite variations of the unit thermal 
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conductance f, and thus the coefficient of ¢ in Equation [7] is a 
function of velocity, and finally of time, since the velocity de- 
pends upon time. An idea of the variation of f with time can be 
obtained as follows: 

If the frequency of velocity variation is not too great, it may be 
postulated that, during a cycle of velocity variation, the value of f 
at any instantaneous velocity v is equal to the value of f which 
would have been obtained at an equal steady velocity.’ With a 
given current flowing through the wire, values of f versus » may 
be obtained at various steady velocities. A typical plot of f ver- 
sus » is shown in Fig. 2, The variation of f with v is complex be- 
cause at zero velocity, free convection, radiation, and conduction 
are important; while at high velocities forced convection con- 
trols. Now, assuming the data for Fig. 2 are applicable at each 
instant of a cycle of variable velocity, the construction of Fig. 3 
would yield the corresponding variations of f with 0. 
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It is noted that although the velocity v becomes negative, the 
value of f is, of course, always positive since heat is transferred re- 
gardless of the direction of air velocity. At the point v = 0 the 
curve of f versus v may not be exactly symmetrical due to free 
convection effects (2, 9, 12, 14). 

Since f varies with time, the coefficient of ¢ in Equation [7] is 
a function of time, and Hquation [8] may be formally written as 


de 
Br PoP (6) evan Qi. laste, Pe eee {8] 
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where 
== 52 
P() = (4 3.4137 ae) 
wc, 
and : 
mys Ge PReak« 
WC, 
The formal solution of Equation [9] is (18) 


@ = eS Pal [QS Pa dg +. ¢| Woe neers [9] 


Unfortunately, this formal solution is of little use, since the 
analytical expression for P(@) is of such a nature that direct in- 
tegration is impossible. Maxwell (15) has solved Equation [8] 
for small periodic-velocity variations about a zero mean by ex- 
panding the voltage ein a Fourier series. This method, however, 
is laborious and not practical for engineering purposes. Fortu- 
nately, another method of solution of Equation [8] is available, 
namely, the method of isoclines (19), which in this particular case 
leads quickly to an engineering answer. 

In the isocline method, Equation [8] is plotted on the e-6 
plane, using constant values of (de/dé) asa parameter. Having 
established loci along which de/d@ is constant, short lines are 
drawn with a slope of de/d@ on these loci. The plane therefore 
becomes covered with these ‘‘direction’’ lines. Having the initial 
condition given locates a point in the plane, and the desired 
solution is found by drawing a curve from this initial point, 
always tangent to the direction lines. This curve will, of course, 
satisfy Equation [8] and also the initial condition and therefore 
is a solution of the equation. 

A difficulty still remains, however, since the function P(@) is 
not known. Inspection of Equation [8] reveals a method of 
obtaining this function by a simple graphical method. 

When the steady-velocity calibration curve of the anemometer 


d 
is obtained, since os = 0, Equation [8] reduces to 


where P is now a function of v rather than 6, since f depends on v. 
The subscript 0 denotes a condition of zero thermal lag. A plot 
of é versus v obtained experimentally then yields 


A typical calibration curve, which can be readily obtained ex- 
perimentally, is shown in Fig. 4. It is noted that Fig. 4 is closely 


akin to Fig. 2. 
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Now if a periodic velocity of very low frequency, Fe — 0, is 


applied to the hot wire, the construction shown in Fig. 5 can be 
made. Projecting the instantaneous velocities on the steady- 
velocity calibration yields the voltage eo, existing across the wire 
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d 
a — 0, Wire No THEerMat Lag 


at each instant of time. Thus on the right of Fig. 5, the curve 
plotted represents the voltage-time variation across the anemome- 
ter if there is no thermal lag. The curve therefore represents 


d 
and is the locus of the isoclines = = 0. Equation [12] allows 


the evaluation of the function P(@). Returning to Equation 
[8] and substituting the foregoing relation 
de Q 
eée= 
dé €0(8) 


d 
Thus curves of constant S = S become 


= *) €(0) = e0(0) (: as ;) pe eee (14] 


Equation [14] represents the locus of isoclines of value S and is 
immediately known once é(6) and Q are known. 

It should be remembered that this method is based on the 
postulate that the heat loss at any velocity v, during a cycle of 
periodic velocity variation, is exactly equal to the heat loss at the 
same velocity during steady-flow operation. This postulate may 
not be exactly true because of the difference in the hydrodynamic 
conditions about the wire in the two cases. Experimental work 
needs to be performed to check this point. 

If the method is experimentally verified, it may be possible to 
correct an experimental response curve by means of Equation 
[14] to obtain e(@), from which the true velocity existing at any 
instant may be determined. 


NUMERICAL EXAMPLES 


Two numerical examples, illustrating the application of the iso- 
cline method of solution, will be presented. The first example 
shows the prediction of the response of the hot-wire anemometer 
subjected to a periodic velocity of 13 fps about a zero mean, 
with a frequency of 5 cycles per sec. The second example shows 
the response of the hot wire to a periodic velocity with an ampli- 
tude of 8 fps about a mean of 20 fps, with a frequency of 5 cycles 
per sec, and also 50 cycles per sec. 
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The following data are available for both examples: 


1 Steady-velocity calibration curve of the wire operating with 
a constant current of 0.140 amp. 
2 The properties of the hot wire as follows: 
(a) Material = platinum 
(b) Diameter = 0.000532 in. (may be calculated from re- 
sistance and resistivity of wire) 
(c) Length of wire = 0.285 in. 
(d) Resistance at room temperature = 4.60 ohms 


The following properties of platinum are used in the calcula- 
tions: 


(a) Density = 1340 lb per cu ft 
(b) Heat capacity = 0.032 Btu/(Ib) (deg F) 


From these data the quantity Q can be calculated 
Q = 2.12 X 10® volts per hr = 59 volts per sec 


Example 1. The velocity cycle shown in Fig. 6, which has an 
amplitude of 13 fps and a frequency of 5 cycles per sec, is pro- 
jected on the steady-state calibration curve, and the zero lag 
response determined. This curve is also the locus of the isocline 
S = 0. To determine the loci of other isoclines, the values of eo 


S 
at each instant are multiplied by ( 59) where S is measured 


in volts persecond. These loci, for several values of S, are shown 
in Fig. 6. Having established the isocline field, the response of 
the anemometer can be determined by starting at any point 
and following the isoclines until the resulting curve becomes peri- 
odic. The periodic curve represents the steady-state behavior of 
the hot wire. The initial point at which the curve is started is 
of no consequence, since the steady-state solution is independent 
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of the initial conditions. The response of the hot wire in the first 
cycle or so represents the transient behavior of the wire before 
reaching steady state. Two initial conditions are shown in Fig. 6, 
both, of course, leading to the same steady-state solution. Fig. 
6 reveals clearly the effect of the heat capacity on the hot-wire 
response. It is interesting to note that due to the nonlinearity of _ 
the steady-state calibration curve, the lag differs at various parts 
of the cycle. Also the wire never heats up to its “still-air”’ tem- 
perature even when the velocity cycle passes through zero. Thus 
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the voltage across the wire never exceeds 1.0 volts, although the 
still-air voltage is 1.65 volts. 

Haample 2. In example 2, shown in Fig. 7, the wire is sub- 
jected to a periodic velocity of 8 fps. The isocline method is 
again used to predict the behavior of the hot wire for frequencies 
of 5 and 50 cycles per sec. For simplicity, the construction lines 
are not shown. The figure reveals that for a frequency of 5 
cycles per sec, the lag and the reduction in amplitude are very 
small, but for 50 cycles per sec, they have become appreciable. 
It is also noted that at higher frequencies the higher harmonics 
in the response curve are damped out and the response curve 
approaches a pure sine wave. A comparison of the lags and 
reductions in amplitude obtained by the isocline method, with 
the approximate equation of Dryden and Keuthe (6), reveals a 
good check for the low-frequency response, and considerable 
divergence for the high-frequency curve. 


CoNCLUSIONS 


A method has been presented which allows a rapid graphical 
prediction of the response of a hot-wire anemometer to periodic 
velocities. From such a prediction, the suitability of the hot 
wire for any particular application can be readily determined. 
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(Owing to travel emergency conditions existing when this paper was presented, discussion will be accepted until February 10, 1946) 
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